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

GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAC  >  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)
                     gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAAc                                                                                                                               <  2:757266/150‑1 (MQ=255)
                        atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTCCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCCCCCCCTGCACCCCCCGGCGCTGGCggaaaa                                                                                                                              >  1:1271802/1‑144 (MQ=255)
                        ataaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAATCGCCCCCGCCCCCCGCGCCCCCCGGctcttgggcctaaa                                                                                                                             >  1:885148/4‑140 (MQ=25)
                         tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCCGCACCCCCGGGCTCTGGCGAACaaaa                                                                                                                           >  2:108870/1‑150 (MQ=255)
                          caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATCGCCGCCCCCCCCCGCACCCCCCGGCTCCGGCGGAAAAAATc                                                                                                                         >  1:836156/1‑151 (MQ=25)
                          caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATc                                                                                                                         <  1:70363/151‑1 (MQ=255)
                            ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCGCCCCCCGCCCCCCCCGGCCCTGGCGAACAAAATcac                                                                                                                       >  2:650003/1‑149 (MQ=25)
                             gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgcc                                                                                                                      <  2:355153/151‑1 (MQ=255)
                                cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGTGGCCAACCCCCGCCGCCCCCCCCACCCCCCCGGTCTGGGGGACAAAAAcgc                                                                                                                       >  2:304968/1‑147 (MQ=255)
                                cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGCGGCCAAACGGCCCCGCCCCCCGCACCCCCCGGCTCTGGCGAAACaaaaagc                                                                                                                        >  1:918576/1‑144 (MQ=255)
                                 gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTg                                                                                                                  <  2:393516/151‑1 (MQ=255)
                                  aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCCCCCCCCCCTCCACCCCCGGGGTCTGGCGAACAAAATctc                                                                                                                       >  1:524208/1‑143 (MQ=255)
                                     aTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGcgcg                                                                                                              <  2:371959/151‑1 (MQ=255)
                                       aaCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCAATCGCCCCCCCCCCCCCCCCCCCCCCGCCCTTGCGCACaaaaaaccccccgcgcgg                                                                                                             >  2:770770/1‑136 (MQ=21)
                                       aaCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCGGGCTCTGGGGAAAAAAAACCCCCCGGCGCggg                                                                                                            >  1:157579/1‑150 (MQ=21)
                                        aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTg                                                                                                           <  2:920589/151‑1 (MQ=255)
                                         cGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCCGGCTCTTGCGAACAAAAACGCCGCCGCGCgggg                                                                                                           >  1:1082568/1‑148 (MQ=21)
                                           agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCCCCCCCCGGCTCTGGCGAACAAAATCGCCCCCCCGCgggcgcc                                                                                                        >  1:837466/1‑146 (MQ=21)
                                             aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                                                      <  2:84114/151‑1 (MQ=255)
                                                   tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATCGCCGCCCCCCCCCCCCCCCCCCGGCTCTGGCGGACAAAAACTCCGCCCCGCGGGGccc                                                                                                        >  1:1237696/1‑143 (MQ=18)
                                                   tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAATAGCCGCCGCCCCCCCCACCCCCCCGGTCTGGCGAAAAAAAACCACGCCGCGCGGGGccccccgc                                                                                                   >  1:812835/1‑143 (MQ=18)
                                                   tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATAGCCGCCGCGCGGTGCCCTCGGcgt                                                                                                 >  1:607063/1‑148 (MQ=255)
                                                    cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCCCCGCCGCCCCCCGCCCCCCCCGGCTCCCGCGGACAAAAACGCCCCTGCGCCGGGccccccgcctct                                                                                               >  2:367880/1‑142 (MQ=18)
                                                    cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCCGGCCCTGGCGAACAAAAACGCCGCTGCGCGGGGCCCCCCGCTTAt                                                                                               >  1:8589/1‑151 (MQ=18)
                                                      gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                                             <  2:1267619/151‑1 (MQ=255)
                                                        ctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCAAACGCCCCCGCCCCCCCCACCCCCCGGCCCTGGCGAACAAAAACGCCGCCGCGCGGGGCCCCCGGCCTATccc                                                                                            >  2:337512/1‑150 (MQ=18)
                                                         tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGGGGCCAATCACCGACCCCCCCCCCCCCCCCCCGCCTTCGGGAACAAAATAGCCCCTGCGCGGGGcccccccctttatccct                                                                                          >  1:848337/1‑137 (MQ=16)
                                                          ctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCCCCACCCCCCGGCTCTGGCGGAAAAAAAAGCCCCTCCGCGGGGCCCCCCGCCTATccc                                                                                            >  2:1110479/1‑148 (MQ=18)
                                                            cATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACCCCGCCGCCCCCCGCACCCCCCGGCTCTGGCGAACAAAAACGCCGCTGCGCGGGGCCCCCCGCTTATccc                                                                                            >  1:604741/1‑146 (MQ=18)
                                                                cTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCAAACCCCCCCGCCCCCCCCACCCCCCGGGTCTGGCGAACAAAAACACCGCCGCGCGGGGCCCCCGGCCTCTCCCCTTCGGccc                                                                                   >  1:1121666/1‑149 (MQ=17)
                                                                    tGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAATCGCCGCCGCCCCCCGCACCCCCCGGCTCCGGCGGACAAAAACGCCGCCGCGCGGGGCCCTCCGCTTATCCCTTACGGCTa                                                                                   >  1:611636/1‑147 (MQ=17)
                                                                             gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                                                       <  2:313174/150‑1 (MQ=255)
                                                                                  gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                                                  <  2:579921/150‑1 (MQ=255)
                                                                                    gCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCCAACGCCCCCCCCCCCCCCCCCCCCCCGCCCCCCCGGACAAAAACGGCGCCGCCCGGGGCCCTCGGCCTATTCCCTACTGCCAcccccgc                                                                            >  2:443080/1‑133 (MQ=9)
                                                                                    gCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCCCCCCCCGCACCCCCCGGCCCTGGCGAACAAAAACGCCGCCGCGCGGTGCCCCCGGCTTATTCCTTACCGCTACCCGGgcgggggcgggg                                                                  >  2:1176563/1‑148 (MQ=14)
                                                                                     cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCCCCCCCTGCACCCCCCGGCTCCGGCGAACAAAAACCCCGCTGCGCGGGGCCCCCCGCTTCTCCCTTACGGCTACCCGGgcgggggcggggg                                                                 >  1:335462/1‑147 (MQ=14)
                                                                                       ttGCTTCATATAATCACTGGCGCAAGGAGCGCGCAAGGGGGGGCCCATCGCCTCCCCCCCCCGCGCCCCCGCGCTCTCGCGGACAAAATCGCCGCCGCCGGGGGCCCCCCGCCTATCCCCTACGGCCTCCCGGgcgggggcgggggga                                                               >  1:604267/1‑145 (MQ=255)
                                                                                        tGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                                            <  1:106253/150‑1 (MQ=255)
                                                                                        tGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCCCCCCCCGCCCCCCCGGGCTCTGGGGGACAAAAACACCGCCCCGCGGCGCCCTCCGCCTATCTCTTCCCGCTATCcgggcgggg                                                                        >  1:607767/1‑137 (MQ=9)
                                                                                         gCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGGCCAACGCCGCCCCCCCCCCCCCCCCCCCGCTCCGGCGGACAAAAAAACCGCCGCGCGGGGCCCCCCGCCTTTCCCCCACCGCTACCCCGgcgggggcggggggaaaacc                                                          >  2:100451/1‑143 (MQ=9)
                                                                                            tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                                                        <  1:836914/150‑1 (MQ=38)
                                                                                            tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAACCGCCGCCGCCCCCTCCCCCCCCCGGCTCTGGCGAACAAAATCGCCCCGGCGCGGTGCCCTCGGCCTATCCCTTACGGCCAACCGGTCcggggcggggggggcttcc                                                         >  2:857754/1‑140 (MQ=14)
                                                                                             cATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                                                      <  1:135393/151‑1 (MQ=37)
                                                                                                   aTCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAATCGCCGCCGCCCCCCGCACCCCCCGGCTCTGGCGGACAAAAACGCCGCTGCGCGGGGCCCTCCGCTTATCCCCTACGGCCAACCGGGCGGGGGCGGGGGAACATCCCCGTaaaa                                                  >  2:559210/1‑149 (MQ=11)
                                                                                                         ggCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                                           <  1:366474/150‑1 (MQ=31)
                                                                                                            gcAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCa                                        <  1:1050018/150‑1 (MQ=31)
                                                                                                            gcAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg                                       <  1:1032552/151‑1 (MQ=31)
                                                                                                              aaGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc                                     <  2:297633/151‑1 (MQ=31)
                                                                                                                ggAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGTTGCCCTCGGCTTATCCCTTACGTCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca                                   <  1:1208033/151‑1 (MQ=14)
                                                                                                                 gAGCGCGCAGGGGGCGGCCAATCGGCGCCGCCCCCTGCACCCCCGGTCTCTGTCGAACAAAATCGCCGGTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac                                  <  1:1008458/151‑1 (MQ=9)
                                                                                                                     gcgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc                              <  2:885184/151‑1 (MQ=31)
                                                                                                                      cgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa                             <  2:882791/151‑1 (MQ=31)
                                                                                                                        cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                            >  1:145619/1‑150 (MQ=31)
                                                                                                                         aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                          <  1:297184/151‑1 (MQ=31)
                                                                                                                         aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGCCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg                           >  2:699056/1‑150 (MQ=14)
                                                                                                                          gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                          <  2:1100525/150‑1 (MQ=31)
                                                                                                                            gggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg                       >  2:1278738/1‑151 (MQ=31)
                                                                                                                                    aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCATTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                >  2:530024/1‑150 (MQ=14)
                                                                                                                                     aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCgg              <  1:94720/151‑1 (MQ=31)
                                                                                                                                       cgccgccgcCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTGCGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCTCTCAGCCCACATCCATGTGGGCTGCCCCGGcc            <  1:395691/151‑1 (MQ=14)
                                                                                                                                             cgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAAGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCCgg      >  1:685310/1‑151 (MQ=14)
                                                                                                                                             cgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAAACGCCCTCACCCCACATCCATGTGGGCTGCCCCCGCCTTc         >  1:625581/1‑148 (MQ=9)
                                                                                                                                              gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCTGGTCGGGCGCGAGGGAACATCCCTGTAAAACGCGCCCTCAGCCCCCATCCATGTGGGCTTCCCCGGCCTTCAg       >  1:349874/1‑149 (MQ=9)
                                                                                                                                              gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg     >  2:151594/1‑151 (MQ=31)
                                                                                                                                               cccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGa    <  2:540026/151‑1 (MQ=31)
                                                                                                                                                  ccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc  >  1:120912/1‑150 (MQ=31)
                                                                                                                                                    |                                                                                                                                                    
GATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAC  >  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: