Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_pgi	464,358	0	T	G	25.4%	118.3 / 55.0	67	intergenic (+50/‑44)	ybaV/fadM	putative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base G (6/11);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGG  >  NC_000913_3_pae_pgi/464211‑464497                                                                                                                                                    |                                                                                                                                            aaGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                                                           <  1:520117/150‑1 (MQ=255)  agagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGACGAGGTCATACCAGTTAt                                                                                                                                          <  2:917475/150‑1 (MQ=255)    agTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGa                                                                                                                                        <  1:3536/150‑1 (MQ=255)       aCGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     <  1:983816/150‑1 (MQ=255)       aCGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     <  3:580329/150‑1 (MQ=255)        cGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctc                                                                                                                                    <  2:216602/150‑1 (MQ=255)         ggTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACAATGTAATCAACTTGCATAGCCGCAATTTTTTGACAGACAGAAGCAGTCATACCAGTTATGACatct                                                                                                                                   >  1:722202/1‑150 (MQ=255)            ccGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttca                                                                                                                                <  3:465216/150‑4 (MQ=255)            ccGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttca                                                                                                                                <  4:120206/150‑4 (MQ=255)              ggttAAAGCCGTGGAGGATCTAAAGCCGGTGCCGGGGGTGGGCAATTTGCTGGTGGAACGTAATCTGGCGGTATTAACACTCGAATTAATTTGCATAGTGGCAATCTTTTGCCAGACTGCAGAGGTCATCACTGGTATGATctctccagt                                                                                                                              <  2:542585/148‑6 (MQ=255)               tttAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACtt                                                                                                                             <  2:554388/150‑1 (MQ=255)                         tgGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaacaa                                                                                                                   <  2:33568/150‑1 (MQ=255)                           gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaa                                                                                                                         >  1:179414/1‑142 (MQ=255)                           gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaa                                                                                                                         <  2:179414/142‑1 (MQ=255)                           gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACg                                                                                                                 >  1:549866/1‑150 (MQ=255)                           gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACg                                                                                                                 >  1:400717/1‑150 (MQ=255)                                  tAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTGATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGtt                                                                                                          <  2:885543/150‑1 (MQ=255)                                   aaaGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTa                                                                                                         <  4:15708/150‑1 (MQ=255)                                   aaaGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTa                                                                                                         >  3:15708/1‑150 (MQ=255)                                         ggTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     <  4:42577/116‑1 (MQ=255)                                         ggTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     >  3:42612/1‑116 (MQ=255)                                         ggTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     >  3:42577/1‑116 (MQ=255)                                         ggTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     <  4:42612/116‑1 (MQ=255)                                               ggggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                                                             <  3:144132/150‑1 (MQ=255)                                                gggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                                                                    <  3:434384/126‑1 (MQ=255)                                                gggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                                                                    >  4:434384/1‑126 (MQ=255)                                                gggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATGACAACAACGTAAGGTTATTGCGCTATGCaa                                                                                            <  1:518588/150‑1 (MQ=255)                                                 ggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                           <  2:531923/150‑1 (MQ=255)                                                       gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                                                           >  2:704700/1‑96 (MQ=255)                                                       gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                                                           <  1:704700/96‑1 (MQ=255)                                                           ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaa                                                                                 >  1:21076/1‑150 (MQ=255)                                                                  gtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACct                                                                                                                                     >  3:146828/1‑91 (MQ=255)                                                                  gtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACct                                                                                                                                     <  4:146828/91‑1 (MQ=255)                                                                      aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGat                                                                      <  4:514950/150‑1 (MQ=255)                                                                          tgaTCTGGCGGTCTTAACCCTGTAATTAATTTGCATAGTGGCAATTCTTTGCCAGACTGAAGAGGTCATACTAGTTATGGGCTCGGTAATTATATCAACAACGTAAGGTTATTGCGCTATGCAAACACAAAACACAGTTCGTGGTTatca                                                                  <  2:243916/148‑1 (MQ=255)                                                                             tCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtct                                                               <  2:408660/150‑1 (MQ=255)                                                                                ggCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGa                                                            <  2:400717/150‑1 (MQ=255)                                                                                 gCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGAc                                                           <  3:428019/150‑1 (MQ=255)                                                                                      aTTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                                                          >  1:329975/1‑114 (MQ=255)                                                                                      aTTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                                                          <  2:329975/114‑1 (MQ=255)                                                                                        tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTAc                                                                                                                               <  1:374331/75‑1 (MQ=255)                                                                                        tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTAc                                                                                                                               >  2:374331/1‑75 (MQ=255)                                                                                          aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAg                                                  >  2:654222/1‑150 (MQ=255)                                                                                            ccTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa                                                >  1:849764/1‑150 (MQ=255)                                                                                                                  gCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCtt                          >  1:1021501/1‑150 (MQ=255)                                                                                                                    aaTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGa                        >  4:304007/1‑150 (MQ=255)                                                                                                                    aaTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGa                        >  3:539428/1‑150 (MQ=255)                                                                                                                          ttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTtct                  <  1:654222/150‑1 (MQ=255)                                                                                                                              cAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa                                          <  3:238319/122‑1 (MQ=255)                                                                                                                              cAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa                                          >  4:238319/1‑122 (MQ=255)                                                                                                                               aGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaag             <  2:549866/150‑1 (MQ=255)                                                                                                                               aGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaag             >  3:433819/1‑150 (MQ=255)                                                                                                                                  cTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa                                          <  2:752958/118‑1 (MQ=255)                                                                                                                                  cTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa                                          >  1:752958/1‑118 (MQ=255)                                                                                                                                    acagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                           >  2:827059/3‑67 (MQ=25)                                                                                                                                    acagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                           <  1:827059/65‑1 (MQ=25)                                                                                                                                      agagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCg      >  4:463709/1‑150 (MQ=255)                                                                                                                                       gagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaagaa           <  4:104828/144‑1 (MQ=255)                                                                                                                                       gagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaagaa           >  3:104828/1‑144 (MQ=255)                                                                                                                                        agGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCt    <  3:438774/150‑1 (MQ=255)                                                                                                                                        agGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                                                                       >  2:850160/1‑67 (MQ=37)                                                                                                                                        agGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                                                                       <  1:850160/67‑1 (MQ=37)                                                                                                                                          gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaa                                                                                            >  1:474472/1‑60 (MQ=255)                                                                                                                                          gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaa                                                                                            <  2:474472/60‑1 (MQ=255)                                                                                                                                          gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTgg  >  2:925072/1‑150 (MQ=255)                                                                                                                                             aTACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTc                                                                            >  1:662293/1‑73 (MQ=255)                                                                                                                                             aTACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTc                                                                            <  2:662293/73‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                            AAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGG  >  NC_000913_3_pae_pgi/464211‑464497

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑