Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,691,169 A→T 5.7% noncoding (172/184 nt) glmY ← sRNA activator of glmS mRNA, glmZ processing antagonist

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,691,1690AT 5.7% 794.3 / 28.1 228noncoding (172/184 nt)glmYsRNA activator of glmS mRNA, glmZ processing antagonist
Reads supporting (aligned to +/- strand):  ref base A (103/112);  new base T (9/4);  total (112/116)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.19e-01

Only 100 of 232 total aligned reads displayed.

TCTCCAGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTC  >  NC_000913/2690897‑2691429
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                    
tctcCAGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTACCCGAGGGGCAGTTCAGCTACAACAAAGCCGGGAAt                                                                                                                                                                                                                                                                    >  2:212290/1‑275 (MQ=255)
                                                 cGTAGTGCGGTTCACCCGCGCCGCCTGGTCGCCAAGCTCATTCAGGCTTTGCCAGGCTTGCCAGGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGcc                                                                                                                                                                                                                   <  4:57113/275‑1 (MQ=255)
                                                     gTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTAGACGTTTCGCCCg                                                                                                                                                                                                                 >  2:57113/1‑273 (MQ=255)
                                                                 agggcgaccTGAACGCGAAGGGGAGTCAGGCTTTGCCAGGCGTGCAAGGGCAGGACCAACAGGGGAAGCAGGATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGGCAGACGCGAATAGCCTGATGCTAACCGGGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAg                                                                                                                                                                                                    <  4:119217/267‑1 (MQ=255)
                                                                      cgctggATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCAAGCACCAAAAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACg                                                                                                                                                                                               <  4:8780/269‑1 (MQ=255)
                                                                        ccTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGAAATTACCCGGATTTGTTATGGAAGAAGGCGGTGCCTAACTCGACGATTCACCCGAAGGGTGATATAGGTaagatg                                                                                                                                                                                            >  1:29182/1‑270 (MQ=255)
                                                                                                            gccatgccaGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGc                                                                                                                                                                                                   >  4:64397/1‑232 (MQ=255)
                                                                                                                tgccaGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACg                                                                                                                                                                                               <  2:64397/232‑1 (MQ=255)
                                                                                                                                    gcagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCg                                                                                                                                                                                                                                                              >  5:1205043/1‑149 (MQ=255)
                                                                                                                                    gcagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCg                                                                                                                                                                                                                                                              >  7:275116/1‑149 (MQ=255)
                                                                                                                                     cagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCgg                                                                                                                                                                                                                                                             <  6:127009/149‑1 (MQ=255)
                                                                                                                                       gAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCt                                                                                                                                                                                                                                                           >  8:926621/1‑149 (MQ=255)
                                                                                                                                       gAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCt                                                                                                                                                                                                                                                           >  5:1070069/1‑149 (MQ=255)
                                                                                                                                         aTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCttt                                                                                                                                                                                                                                                         <  6:1031951/149‑1 (MQ=255)
                                                                                                                                          tCAGCAAAAATGCCAGCATTATCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTg                                                                                                                                                                                                                                                        <  7:1338310/149‑1 (MQ=255)
                                                                                                                                           cAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCTTTGt                                                                                                                                                                                                                                                       <  8:669592/149‑1 (MQ=255)
                                                                                                                                            aGCAAACATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGtt                                                                                                                                                                                                                                                      <  7:136139/149‑1 (MQ=255)
                                                                                                                                              cAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCCCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTAt                                                                                                                                                                                                                                                    <  8:690540/149‑1 (MQ=255)
                                                                                                                                               aaaaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATg                                                                                                                                                                                                                                                   >  5:536393/1‑149 (MQ=255)
                                                                                                                                                aaaaTGCCAGGATATACCAGTGTCGTAATGGGGGGGGAAAAACGGGCCGGCGTTTCAAGGGGTTACTCGCGGAAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTAAGATACAACAAAGCCTGGAAGTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGc                                                                                                                     <  4:212290/274‑1 (MQ=255)
                                                                                                                                                  aaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAGGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGaa                                                                                                                                                                                                                                                <  6:943102/149‑1 (MQ=255)
                                                                                                                                                  aaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCTTTGttgtatct                                                                                                                                                                                                                                                <  8:570945/149‑7 (MQ=39)
                                                                                                                                                       cAGCATTACCAGTTGACGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTAc                                                                                                                                                                                                                                                                 <  7:611825/127‑1 (MQ=255)
                                                                                                                                                        aGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCttt                                                                                                                                                                                                                                                         >  2:160653/1‑134 (MQ=255)
                                                                                                                                                         gCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCg                                                                                                                                                                                                                                         <  6:219627/149‑1 (MQ=255)
                                                                                                                                                           aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCAGt                                                                                                                                                                                                                                       <  7:722649/149‑1 (MQ=255)
                                                                                                                                                            ttACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTa                                                                                                                                                                                                                                                     <  4:160653/134‑1 (MQ=255)
                                                                                                                                                                ccAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                                                                                                                                   >  6:344592/1‑149 (MQ=255)
                                                                                                                                                                  agatGTCGTAATGACCGGGGAAAAACGGGCTAGCATTTCAAGGTGTTACTCTCGTCAGACGTGAATAGCCTGATGCTAACCGAGGGGAAGGTCAGATACCACAAAACCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTaa                                                                                                                                                                                                                                 <  8:899578/146‑1 (MQ=255)
                                                                                                                                                                  aGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTaa                                                                                                                                                                                                                                 >  8:154420/1‑149 (MQ=255)
                                                                                                                                                                  aGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTaa                                                                                                                                                                                                                                 <  7:189096/149‑1 (MQ=255)
                                                                                                                                                                   gTTGTCGTAATGAGCGGGGAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACt                                                                                                                                                                                                                               >  5:427485/1‑149 (MQ=255)
                                                                                                                                                                        gggtagggggggggAAAAACGGGCCAGCGTTTCAAGGTGTGACTCTCGTCAGACGCGAAGAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCa                                                                                             <  4:301467/265‑1 (MQ=255)
                                                                                                                                                                          tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTGAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCa                                                                                             <  2:299150/273‑1 (MQ=255)
                                                                                                                                                                          tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGGATAAGGCGGTGCCTAACTCGACGt                                                                                                                                                                                                                         >  7:243405/1‑149 (MQ=255)
                                                                                                                                                                          tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGt                                                                                                                                                                                                                         >  5:716295/1‑149 (MQ=255)
                                                                                                                                                                             tGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                                                                                                                                      <  6:777813/149‑1 (MQ=255)
                                                                                                                                                                             tGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                                                                                                                                      <  8:1013065/149‑1 (MQ=255)
                                                                                                                                                                                    ggAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGAt                                                                                                                                                                                                               <  7:926621/149‑1 (MQ=255)
                                                                                                                                                                                     gAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATg                                                                                                                                                                                                              >  6:413262/1‑149 (MQ=255)
                                                                                                                                                                                        aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                                           <  6:589891/149‑1 (MQ=255)
                                                                                                                                                                                        aaaCGGGCCAGCGTTTCAAGGTGTTACTCCCGTCAGACTCGAATAGCATGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTGATGGTATAAGGCGGTGCCTAACTAGACGTTTCGCGCGATGGtt                                                                                                                                                                                                           <  6:759488/149‑1 (MQ=255)
                                                                                                                                                                                         aaCGGGCCAGTGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTg                                                                                                                                                                                                          <  7:310784/149‑1 (MQ=255)
                                                                                                                                                                                         aaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCg                                                                                                                                                                                                                                         >  7:336063/1‑118 (MQ=255)
                                                                                                                                                                                         aaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAAATCGACGTTTCGCCCGATGGTTg                                                                                                                                                                                                          >  7:826452/1‑149 (MQ=255)
                                                                                                                                                                                          aCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                                           <  6:198015/147‑1 (MQ=255)
                                                                                                                                                                                          aCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                                           >  5:198015/1‑147 (MQ=255)
                                                                                                                                                                                          aCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGa                                                                                                                                                                                                         <  7:154420/149‑1 (MQ=255)
                                                                                                                                                                                                 aGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                                           <  8:310243/140‑1 (MQ=255)
                                                                                                                                                                                                 aGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCt                                                                                                                                                                                                  >  7:876399/1‑149 (MQ=255)
                                                                                                                                                                                                 aGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCt                                                                                                                                                                                                  <  6:687223/149‑1 (MQ=255)
                                                                                                                                                                                                    gTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACg                                                                                                                                                                                               >  7:1103973/1‑149 (MQ=255)
                                                                                                                                                                                                    gTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACg                                                                                                                                                                                               >  7:412533/1‑149 (MQ=255)
                                                                                                                                                                                                            gtgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATc                                                                                                                                                                                       <  6:960126/149‑1 (MQ=255)
                                                                                                                                                                                                            gtgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATc                                                                                                                                                                                       >  8:724750/1‑149 (MQ=255)
                                                                                                                                                                                                             tgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGAtata                                                                                                                                                                                                     <  7:856743/134‑1 (MQ=255)
                                                                                                                                                                                                             tgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCa                                                                                                                                                                                      >  5:295558/1‑149 (MQ=255)
                                                                                                                                                                                                             tgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCa                                                                                                                                                                                      <  8:5966/149‑1 (MQ=255)
                                                                                                                                                                                                              gtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGttt                                                                                                                                                                                                                       >  8:1128404/1‑115 (MQ=255)
                                                                                                                                                                                                                  ctctCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGt                                                                                                                                                                                 >  8:435891/1‑149 (MQ=255)
                                                                                                                                                                                                                   tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAAc                                                                                                                                                                                                                                <  7:66998/101‑1 (MQ=255)
                                                                                                                                                                                                                   tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGatat                                                                                                                                                                                                      >  5:1429407/1‑127 (MQ=255)
                                                                                                                                                                                                                   tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGtt                                                                                                                                                                                >  6:669358/1‑149 (MQ=255)
                                                                                                                                                                                                                   tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGtt                                                                                                                                                                                >  5:790924/1‑149 (MQ=255)
                                                                                                                                                                                                                       gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGAc                                                                                                                                                                            <  5:669358/149‑1 (MQ=255)
                                                                                                                                                                                                                       gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATa                                                             >  2:249738/1‑260 (MQ=255)
                                                                                                                                                                                                                       gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATa                                                             <  4:249738/260‑1 (MQ=255)
                                                                                                                                                                                                                       gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTgg                                              >  2:141514/1‑275 (MQ=255)
                                                                                                                                                                                                                           gACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                                                                                                                                        <  8:526617/149‑1 (MQ=255)
                                                                                                                                                                                                                            aCGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggcag                                                                                                                                                                       >  7:484136/1‑149 (MQ=255)
                                                                                                                                                                                                                                      ccTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                                                                                                                                             <  5:964140/149‑1 (MQ=255)
                                                                                                                                                                                                                                        tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                                                                                                                                           <  8:103301/149‑1 (MQ=255)
                                                                                                                                                                                                                                              tAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCa                                                                                                                                                     >  8:1388865/1‑149 (MQ=255)
                                                                                                                                                                                                                                              tAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCa                                                                                                                                                     >  6:863931/1‑149 (MQ=255)
                                                                                                                                                                                                                                                aCCGAGGGGAAGTTCAGATACATCAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCTACGTTTCTCCCGTTGGTTTTTATTTCTACGCTTATATCAGAAGtttg                                                                                                                                                                              >  2:132502/1‑120 (MQ=255)
                                                                                                                                                                                                                                                aCCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCt                       >  2:330680/1‑273 (MQ=255)
                                                                                                                                                                                                                                                 ccGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTc                                                                                                                                                  <  5:863931/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTg                    <  4:141514/275‑1 (MQ=255)
                                                                                                                                                                                                                                                   gAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                                           <  5:411614/90‑1 (MQ=255)
                                                                                                                                                                                                                                                     ggggAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAt                                                                                                                                              <  8:166688/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      gggAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAtt                                                                                                                                             >  7:1282764/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       ggAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAttt                                                                                                                                            <  5:888781/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       ggAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACGCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAttt                                                                                                                                            <  8:826452/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            ttCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGGTATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGt                                                                                                                                       >  7:1175928/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             tCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTa                                                                                                                                      >  5:1288305/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              cAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAg                                                                                                                                     <  8:1346482/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aacaaAGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGATACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                                                                                                                              >  8:25016/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     aacaaAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                                                                                                                              >  8:1280666/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       caaAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCg                                                                                                                            <  6:437576/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aaaGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGa                                                                                                                           >  5:46112/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         aaGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa                                                                                                                          >  7:210891/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         aaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa                                                                                                                          <  8:192761/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa                                                                                                                          >  8:355396/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          aGCCGGGTAATTCCCGGCTTTGTTGTGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                                         >  5:937083/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                                         >  6:237905/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          aGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                                         <  8:280680/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                                         <  8:393508/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAgg                                                                                                                        <  6:975684/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAgggg                                                                                                                      <  7:435891/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aaTTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGa                                                                                                                  <  7:355396/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                    
TCTCCAGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTC  >  NC_000913/2690897‑2691429

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: