Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,229,2040TC39.0% 94.9 / 19.6 77intergenic (+139/+234)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base T (18/29);  new base C (30/0);  total (48/29)
Fisher's exact test for biased strand distribution p-value = 4.16e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAC  >  NC_000913/2229056‑2229351
                                                                                                                                                    |                                                                                                                                                    
gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCCCCGCCGCCCCCCGc                                                                                                                                                    >  1:786781/1‑151 (MQ=255)
        aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTGATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCg                                                                                                                                             <  1:241786/150‑1 (MQ=255)
        aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCgg                                                                                                                                            >  1:487124/1‑151 (MQ=255)
         aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCCATCTcccccccccccccccccccc                                                                                                                                              >  2:1060236/1‑148 (MQ=38)
         aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATCGCCGCCGCCCCCCGCACCCCCcgg                                                                                                                                           >  2:307446/1‑151 (MQ=255)
          aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATCGCCGCCCCCCCCCCCAccccc                                                                                                                                              >  2:1070966/1‑147 (MQ=255)
           attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTTGCGCAAGGAGCGCGGAAGGGGGGGCCAATCGCCCCCCCCCCCCTCACCCCCcgg                                                                                                                                           >  1:1100465/1‑149 (MQ=255)
           attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGTACCCCCGGGct                                                                                                                                         <  1:203805/151‑1 (MQ=255)
              aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGctct                                                                                                                                       <  1:1024636/150‑1 (MQ=255)
                 taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCTCACGGGGCGGCCAATGGCCGCCCCCCCCTGCACCCCCCGGCTCCGGc                                                                                                                                    >  1:36886/1‑150 (MQ=255)
                  aatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCCCCCCCGGGCTCTGGCg                                                                                                                                   >  1:394128/1‑150 (MQ=255)
                    tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCAAGCGCCGCCCCCCCCCGCACCCCCCGGcccttg                                                                                                                                     >  1:673316/1‑144 (MQ=255)
                     gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AAc                                                                                                                               <  2:757266/150‑1 (MQ=255)
                        atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTCCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCCCCCCCTGCACCCCCCGGCGCTGGCggaaaa                                                                                                                              >  1:1271802/1‑144 (MQ=255)
                        ataaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAATCGCCCCCGCCCCCCGCGCCCCCCGGctcttgggcctaaa                                                                                                                             >  1:885148/4‑140 (MQ=25)
                         tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCCGCACCCCCGGGCTCTGGCG‑AACaaaa                                                                                                                           >  2:108870/1‑150 (MQ=255)
                          caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATCGCCGCCCCCCCCCGCACCCCCCGGCTCCGGCG‑GAAAAAATc                                                                                                                         >  1:836156/1‑151 (MQ=25)
                          caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATc                                                                                                                         <  1:70363/151‑1 (MQ=255)
                            ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCGCCCCCCGCCCCCCCCGGCCCTGGCG‑AACAAAATcac                                                                                                                       >  2:650003/1‑149 (MQ=25)
                             gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATcgcc                                                                                                                      <  2:355153/151‑1 (MQ=255)
                                cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGTGGCCAACCCCCGCCGCCCCCCCCACCCCCCCGGTCTGGGG‑GACAAAAAcgc                                                                                                                       >  2:304968/1‑147 (MQ=255)
                                cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGCGGCCAAACGGCCCCGCCCCCCGCACCCCCCGGCTCTGGCGAAACaaaaagc                                                                                                                        >  1:918576/1‑144 (MQ=255)
                                 gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTg                                                                                                                  <  2:393516/151‑1 (MQ=255)
                                  aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCCCCCCCCCCTCCACCCCCGGGGTCTGGCG‑AACAAAATctc                                                                                                                       >  1:524208/1‑143 (MQ=255)
                                     aTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGcgcg                                                                                                              <  2:371959/151‑1 (MQ=255)
                                       aaCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCAATCGCCCCCCCCCCCCCCCCCCCCCCGCCCTTGCG‑CACaaaaaaccccccgcgcgg                                                                                                             >  2:770770/1‑136 (MQ=21)
                                       aaCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCGGGCTCTGGGG‑AAAAAAAACCCCCCGGCGCggg                                                                                                            >  1:157579/1‑150 (MQ=21)
                                        aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTg                                                                                                           <  2:920589/151‑1 (MQ=255)
                                         cGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCCGGCTCTTGCG‑AACAAAAACGCCGCCGCGCgggg                                                                                                           >  1:1082568/1‑148 (MQ=21)
                                           agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCCCCCCCCGGCTCTGGCG‑AACAAAATCGCCCCCCCGCgggcgcc                                                                                                        >  1:837466/1‑146 (MQ=21)
                                             aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTc                                                                                                      <  2:84114/151‑1 (MQ=255)
                                                   tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATCGCCGCCCCCCCCCCCCCCCCCCGGCTCTGGCG‑GACAAAAACTCCGCCCCGCGGGGccc                                                                                                        >  1:1237696/1‑143 (MQ=18)
                                                   tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATAGCCGCCGCCCCCCCCACCCCCCCGGTCTGGCG‑AAAAAAAACCACGCCGCGCGGGGccccccgc                                                                                                   >  1:812835/1‑143 (MQ=18)
                                                   tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATAGCCGCCGCGCGGTGCCCTCGGcgt                                                                                                 >  1:607063/1‑148 (MQ=255)
                                                    cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCCCCGCCGCCCCCCGCCCCCCCCGGCTCCCGCG‑GACAAAAACGCCCCTGCGCCGGGccccccgcctct                                                                                               >  2:367880/1‑142 (MQ=18)
                                                    cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCCGGCCCTGGCG‑AACAAAAACGCCGCTGCGCGGGGCCCCCCGCTTAt                                                                                               >  1:8589/1‑151 (MQ=18)
                                                      gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                                             <  2:1267619/151‑1 (MQ=255)
                                                        ctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCAAACGCCCCCGCCCCCCCCACCCCCCGGCCCTGGCG‑AACAAAAACGCCGCCGCGCGGGGCCCCCGGCCTATccc                                                                                            >  2:337512/1‑150 (MQ=18)
                                                         tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGGGGCCAATCACCGACCCCCCCCCCCCCCCCCCGCCTTCGGG‑AACAAAATAGCCCCTGCGCGGGGcccccccctttatccct                                                                                          >  1:848337/1‑137 (MQ=16)
                                                          ctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCCCCACCCCCCGGCTCTGGCG‑GAAAAAAAAGCCCCTCCGCGGGGCCCCCCGCCTATccc                                                                                            >  2:1110479/1‑148 (MQ=18)
                                                            cATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACCCCGCCGCCCCCCGCACCCCCCGGCTCTGGCG‑AACAAAAACGCCGCTGCGCGGGGCCCCCCGCTTATccc                                                                                            >  1:604741/1‑146 (MQ=18)
                                                                cTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCAAACCCCCCCGCCCCCCCCACCCCCCGGGTCTGGCG‑AACAAAAACACCGCCGCGCGGGGCCCCCGGCCTCTCCCCTTCGGccc                                                                                   >  1:1121666/1‑149 (MQ=17)
                                                                    tGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAATCGCCGCCGCCCCCCGCACCCCCCGGCTCCGGCG‑GACAAAAACGCCGCCGCGCGGGGCCCTCCGCTTATCCCTTACGGCTa                                                                                   >  1:611636/1‑147 (MQ=17)
                                                                             gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                                                       <  2:313174/150‑1 (MQ=255)
                                                                                  gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                                                  <  2:579921/150‑1 (MQ=255)
                                                                                    gCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCCAACGCCCCCCCCCCCCCCCCCCCCCCGCCCCCCCG‑GACAAAAACGGCGCCGCCCGGGGCCCTCGGCCTATTCCCTACTGCCAcccccgc                                                                            >  2:443080/1‑133 (MQ=9)
                                                                                    gCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCCCCCCCCGCACCCCCCGGCCCTGGCG‑AACAAAAACGCCGCCGCGCGGTGCCCCCGGCTTATTCCTTACCGCTACCCGGgcgggggcgggg                                                                  >  2:1176563/1‑148 (MQ=14)
                                                                                     cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCCCCCCCTGCACCCCCCGGCTCCGGCG‑AACAAAAACCCCGCTGCGCGGGGCCCCCCGCTTCTCCCTTACGGCTACCCGGgcgggggcggggg                                                                 >  1:335462/1‑147 (MQ=14)
                                                                                       ttGCTTCATATAATCACTGGCGCAAGGAGCGCGCAAGGGGGGGCCCATCGCCTCCCCCCCCCGCGCCCCCGCGCTCTCGCG‑GACAAAATCGCCGCCGCCGGGGGCCCCCCGCCTATCCCCTACGGCCTCCCGGgcgggggcgggggga                                                               >  1:604267/1‑145 (MQ=255)
                                                                                        tGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                                            <  1:106253/150‑1 (MQ=255)
                                                                                        tGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCCCCCCCCGCCCCCCCGGGCTCTGGGG‑GACAAAAACACCGCCCCGCGGCGCCCTCCGCCTATCTCTTCCCGCTATCcgggcgggg                                                                        >  1:607767/1‑137 (MQ=9)
                                                                                         gCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGGCCAACGCCGCCCCCCCCCCCCCCCCCCCGCTCCGGCG‑GACAAAAAAACCGCCGCGCGGGGCCCCCCGCCTTTCCCCCACCGCTACCCCGgcgggggcggggggaaaacc                                                          >  2:100451/1‑143 (MQ=9)
                                                                                            tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                                                        <  1:836914/150‑1 (MQ=38)
                                                                                            tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAACCGCCGCCGCCCCCTCCCCCCCCCGGCTCTGGCG‑AACAAAATCGCCCCGGCGCGGTGCCCTCGGCCTATCCCTTACGGCCAACCGGTCcggggcggggggggcttcc                                                         >  2:857754/1‑140 (MQ=14)
                                                                                             cATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                                                      <  1:135393/151‑1 (MQ=37)
                                                                                                   aTCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAATCGCCGCCGCCCCCCGCACCCCCCGGCTCTGGCG‑GACAAAAACGCCGCTGCGCGGGGCCCTCCGCTTATCCCCTACGGCCAACCGGGCGGGGGCGGGGGAACATCCCCGTaaaa                                                  >  2:559210/1‑149 (MQ=11)
                                                                                                         ggCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                                           <  1:366474/150‑1 (MQ=31)
                                                                                                            gcAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCa                                        <  1:1050018/150‑1 (MQ=31)
                                                                                                            gcAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg                                       <  1:1032552/151‑1 (MQ=31)
                                                                                                              aaGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc                                     <  2:297633/151‑1 (MQ=31)
                                                                                                                ggAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGTTGCCCTCGGCTTATCCCTTACGTCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca                                   <  1:1208033/151‑1 (MQ=14)
                                                                                                                 gAGCGCGCAGGGGGCGGCCAATCGGCGCCGCCCCCTGCACCCCCGGTCTCTGTCG‑AACAAAATCGCCGGTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac                                  <  1:1008458/151‑1 (MQ=9)
                                                                                                                     gcgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc                              <  2:885184/151‑1 (MQ=31)
                                                                                                                      cgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa                             <  2:882791/151‑1 (MQ=31)
                                                                                                                        cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                            >  1:145619/1‑150 (MQ=31)
                                                                                                                         aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                          <  1:297184/151‑1 (MQ=31)
                                                                                                                         aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGCCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg                           >  2:699056/1‑150 (MQ=14)
                                                                                                                          gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                          <  2:1100525/150‑1 (MQ=31)
                                                                                                                            gggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg                       >  2:1278738/1‑151 (MQ=31)
                                                                                                                                    aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCATTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                >  2:530024/1‑150 (MQ=14)
                                                                                                                                     aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCgg              <  1:94720/151‑1 (MQ=31)
                                                                                                                                       cgccgccgcCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTGCGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCTCTCAGCCCACATCCATGTGGGCTGCCCCGGcc            <  1:395691/151‑1 (MQ=14)
                                                                                                                                             cgccCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAAGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCCgg      >  1:685310/1‑151 (MQ=14)
                                                                                                                                             cgccCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAAACGCCCTCACCCCACATCCATGTGGGCTGCCCCCGCCTTc         >  1:625581/1‑148 (MQ=9)
                                                                                                                                              gccCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCTGGTCGGGCGCGAGGGAACATCCCTGTAAAACGCGCCCTCAGCCCCCATCCATGTGGGCTTCCCCGGCCTTCAg       >  1:349874/1‑149 (MQ=9)
                                                                                                                                              gccCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg     >  2:151594/1‑151 (MQ=31)
                                                                                                                                               cccccTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGa    <  2:540026/151‑1 (MQ=31)
                                                                                                                                                  ccTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc  >  1:120912/1‑150 (MQ=31)
                                                                                                                                                    |                                                                                                                                                    
GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAC  >  NC_000913/2229056‑2229351

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