Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2870GC37.9% 114.6 / 34.6 89intergenic (+93/+136)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (17/37);  new base C (33/0);  total (52/37)
Fisher's exact test for biased strand distribution p-value = 1.80e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCC  >  NC_000913/3423145‑3423437
                                                                                                                                               |                                                                                                                                                          
acccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAAACCATTAATTGCATGCTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCCCGCCCCCCCCCCT‑c                                                                                                                                                   >  2:1169389/3‑151 (MQ=255)
     aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGc                                                                                                                                              <  2:1333715/151‑1 (MQ=255)
     aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCACCCC‑CCcgg                                                                                                                                               >  1:725151/1‑150 (MQ=255)
      aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCCCCCC‑CCGGGcc                                                                                                                                             >  1:266711/1‑151 (MQ=255)
       aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCACCCC‑CCCGGcc                                                                                                                                             >  1:849782/1‑150 (MQ=255)
        tCAGAGCGTACGAGGGCATCTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCACTTTATGGCGGTC‑GGCCTGCccccccccacccccg                                                                                                                                                >  2:1599010/1‑140 (MQ=255)
          agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGCATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCCCGCCCGCCACCCT‑CCggg                                                                                                                                               >  2:709585/1‑145 (MQ=255)
          agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑Gtt                                                                                                                                         <  2:1020349/150‑1 (MQ=255)
           gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GCCCCGGCCGCCCCCCC‑CCCGGcccgtg                                                                                                                                         >  2:904631/1‑146 (MQ=255)
              cGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCCGCCCGCCACCCC‑CCcgg                                                                                                                                               >  2:742935/1‑141 (MQ=255)
               gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGGTCGCCCCCCAACCC‑CCTGGCC‑Gt                                                                                                                                          >  2:1299687/1‑144 (MQ=37)
                  cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑aa                                                                                                                               <  2:78504/151‑1 (MQ=255)
                    aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACg                                                                                                                             <  1:104597/151‑1 (MQ=255)
                     gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGt                                                                                                                            <  2:250902/151‑1 (MQ=255)
                        cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCGCCCCCCCCCCC‑CCCGGCC‑GGTGCTTCGC‑CACAtt                                                                                                                           >  1:1035337/1‑149 (MQ=25)
                        cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTTCCCCCCCCCCC‑CCCGGCC‑GCTGGtt                                                                                                                                     >  2:797947/1‑140 (MQ=37)
                         aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑aa                                                                                                                       <  1:1047441/151‑1 (MQ=255)
                         aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTATGGCAGTTTATGGCGGGCGGTCCT‑CCCTCCACCCC‑CCCGTCC‑GGGTCTTCGcccagctta                                                                                                                          >  2:1288614/1‑142 (MQ=255)
                             ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCTCCCCCCCCCCC‑CCCGGCC‑GGTGCTTCGC‑AACGTTC‑AAATcc                                                                                                                   >  1:1346088/1‑151 (MQ=25)
                              tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCTGCCCGCCCCCCC‑CCCGGCC‑GGTGCTTCGC‑AACGTTC‑AAATccc                                                                                                                  >  1:416553/1‑150 (MQ=255)
                                  gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCCTGCCCCCCCCCCC‑CCCGGCC‑CTTTCTTCGC‑AACATTC‑AAAACCGCTcc                                                                                                              >  1:1135709/1‑151 (MQ=255)
                                    aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTcccc                                                                                                            >  1:1388850/1‑150 (MQ=255)
                                     ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAAAcccccccc                                                                                                             >  1:766768/1‑149 (MQ=21)
                                           ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGAt                                                                                                      <  1:1466829/150‑1 (MQ=255)
                                            cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAAACCGCTCCccgcggggtt                                                                                                    >  1:1272993/1‑147 (MQ=18)
                                            cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCC‑CCCGGCC‑GGTGCTTCGC‑AACGTTT‑AAATCCCCTCCccgcggggtt                                                                                                    >  1:1311088/1‑147 (MQ=255)
                                             aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTGCCGCCCGCCCCCCC‑CCGCGCC‑GTGGGTTCCC‑AACGCTT‑CCATCCGCCCCCCGCCGGttt                                                                                                    >  1:397191/1‑150 (MQ=255)
                                             aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTg                                                                                                   <  1:323719/151‑1 (MQ=255)
                                               tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTGCCGTCCGCCACCCC‑CCCGGCC‑GTTTCTTTGC‑TACGTTCAAAACCCGCCCCCGCCGGttt                                                                                                     >  2:539827/1‑148 (MQ=18)
                                                  ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCACCCC‑CCCGGCC‑GTTGCTTCGC‑AACGGTC‑AAAACCGCCCCCCGCGGAttttttc                                                                                                >  1:1455465/1‑145 (MQ=18)
                                                   tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTGCCGCCCCCCACCCC‑CCCGGCC‑GGTGCTGCGC‑AACGGTC‑ACAACCCCTCCCCCCGCATTTGTTCttc                                                                                             >  2:148815/1‑149 (MQ=18)
                                                    gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCCTCCCCCCACCCC‑CCCGGCC‑GTTTCTTCCC‑AAAGCTC‑AAAATCGCCCCCCGCGGGTTTGTTCTAc                                                                                             >  2:1246775/1‑150 (MQ=255)
                                                        gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCgg                                                                                         <  1:175471/150‑1 (MQ=255)
                                                         ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCGTCCCCCGCCCT‑CCCGGGC‑GTTGCTTTGA‑AACCTTC‑AAAACACCTCCCCGGGGGATTGTCCTACTc                                                                                           >  2:863213/1‑147 (MQ=17)
                                                          tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCACCCC‑CCCGGCC‑GGTAGTTCGC‑AACCATC‑CAATTCCCACCCCGCCGCTTTGTCCCACTCcgg                                                                                        >  2:1512243/1‑149 (MQ=17)
                                                               cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtgt                                                                                 <  1:477220/151‑1 (MQ=255)
                                                               cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                  <  2:1324881/150‑1 (MQ=255)
                                                                  aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGGCTGCCCGCCCCCCC‑CCCGGCC‑GTTGCTTTGC‑AACCTTT‑AAATCCCCTCCCGGCGGAGTTGTTCTACTCcgggg                                                                                      >  1:732653/1‑141 (MQ=17)
                                                                       aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCCTCCCCCCCCCCCC‑CCCGGCC‑GGTTTGTTGC‑AACCTTC‑ACATTCGCCCCCCGGGGGTTTGTCCCAccctgg                                                                                        >  2:845058/1‑131 (MQ=255)
                                                                       aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGA‑GGCCTGCCCGCCACCCC‑CCCGGCC‑GGTGCTTCCC‑CAAGTTT‑ACAACTGCCCCCCGCGGGTTTTTTCTACTCTCGAGGGGGTTCTCCGac                                                                         >  2:1194788/1‑151 (MQ=255)
                                                                        aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGGCCTCCCCCCCCCCT‑CCCGGCC‑GGGTCTTTGC‑AACCATC‑AAATACGCCCCCCGCGGGGTTGTTCTTCTCGGGAGGGGGttttcc                                                                            >  2:1106600/1‑143 (MQ=255)
                                                                        aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCC‑GGCCTGCCCCCCCCCCTCCCGGGCC‑GCT‑CTTCGC‑AAAATTT‑AAATACCCCCCCCGCGGGGTTGTCCTACTCCGGGGAgggtgctccc                                                                           >  2:1259320/1‑142 (MQ=255)
                                                                          aatgCCATCCTTTCGGCTGTCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTCCCCGCCACCCT‑CCGGGTC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGTATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                      <  2:237985/148‑1 (MQ=17)
                                                                          aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAAACCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                      <  1:1173829/151‑1 (MQ=255)
                                                                            ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                    <  2:884933/151‑1 (MQ=255)
                                                                            ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                    <  2:1424295/151‑1 (MQ=255)
                                                                            ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                    <  2:109735/151‑1 (MQ=255)
                                                                              ccATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTTCCCCCCCCCCT‑CCCGGCC‑GTTGCTTCCC‑AACGTTC‑AAATCCCCTCCCCGCGGGTTTGTCCTACTCCGGGGAGTGGTCACCGccaaacaa                                                                   >  1:1573775/1‑150 (MQ=17)
                                                                               cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                 <  1:323475/151‑1 (MQ=255)
                                                                               cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCC‑CCGGGCC‑GTTGCTTCGC‑CAAGATC‑AAATCCCCCCCCCGCGGGTTTGTCCTACTCCGGGGAGTGTTTACCCacaaacaa                                                                   >  1:954444/1‑149 (MQ=17)
                                                                               cATCCTTTCGGATGGCCTTTCGCTCGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                 <  1:1283453/151‑1 (MQ=255)
                                                                                  ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGGCCACCAACCC‑CCCAGCC‑GTTGCCTACC‑AAACGTC‑CAATCCCCTCCCCGCGGGTTTTTTCCACTCCGGGCAGAGGTCACCCacaaagaaaaaa                                                               >  1:963655/1‑146 (MQ=9)
                                                                                   cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGGTTTGC‑CACCTTC‑AAATACCCCCCCCGCGGGTTTTTTCTTCTCCGGGGAGGGTTTACCCacaaaaaa                                                                   >  1:468817/1‑145 (MQ=255)
                                                                                    tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCACCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAAACCCCCCCCCGCGGGTTTTTTCTACTCGGGGGAGTGTTCACCCAAAAAAAACAGATa                                                             >  1:6686/1‑150 (MQ=14)
                                                                                     ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCACCCC‑CCCGGCC‑GTTTCTTCCC‑AACGGTC‑AAATCCCCTCCCCGCGGGTTTTTCCCCCTCGGGGGAGAGGTCACCCACAAACAAACGATaa                                                            >  1:917828/1‑150 (MQ=14)
                                                                                      tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCT‑CCCGGCCGGTTGTTTC‑C‑AAAGTTT‑AAAACCCCCCCCCGCGGATTTGTCCTACTCCGGAGAGGGTTTACCCACCAACAAAAGAaaaa                                                           >  1:1456553/1‑150 (MQ=9)
                                                                                       cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCCCTCCCCGCGGGTTTTTCCTACTCGGGGGAGTGGTTACCGACAAACAACAGATaaaa                                                          >  1:28093/1‑150 (MQ=14)
                                                                                        ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACa                                                        <  1:1146701/151‑1 (MQ=255)
                                                                                        ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCCGCGGGTTTTTCCTACTCGGGAGAGGGTTCACCCACAAACAACAGATaaaa                                                          >  1:603749/1‑149 (MQ=14)
                                                                                          aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaa                                                      <  2:1153615/151‑1 (MQ=255)
                                                                                          aTGGCCTTTCCCTTGATTTGATGTCTGGCAGTTTATGGCGGTC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaa                                                      <  2:814750/151‑1 (MQ=255)
                                                                                               cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTT‑AAAACCCCTCCCCGCGGGTTTTTCCTACTCCGGGGAGGGTTTACCCACAAACAAAAGATAaaaaaaaa                                                     >  1:1148344/1‑147 (MQ=14)
                                                                                                tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTTCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCCCCCCGCGGGTTTTTCCTACTCCGGAGAGTGTTCACCGACAAACAACAGATAAAAAAAAAAGccc                                                >  1:1282376/1‑151 (MQ=14)
                                                                                                tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTTCCCCCCCCCCC‑CCCCGCC‑GTTGCTTCCC‑AACGTTT‑CAAACCCCCCCCCGCGGGTTTTTTCTCCTCCGGGGGGGGTTTGCCCCCCAACAACAGATAaaaaaaaaagac                                                 >  2:1454432/1‑146 (MQ=9)
                                                                                                 ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCCCCCT‑CCGGGCC‑GGTGCTTCTC‑AACCTTC‑AAAAACGCCCCCCGCGGATCTGTCCTACTCGGGGGGGGGGTCACCGACAAACAACCGATAaaaaaaaa                                                     >  2:1575074/1‑145 (MQ=14)
                                                                                                 ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACGGATAAAACAAAAGGcc                                                 <  1:802482/149‑1 (MQ=255)
                                                                                                  tCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCa                                               <  1:1030182/150‑1 (MQ=255)
                                                                                                   cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                              <  2:101309/150‑1 (MQ=255)
                                                                                                   cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGt                                             <  1:1300633/151‑1 (MQ=40)
                                                                                                   cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTCCCCCCCACCCC‑CCCGGGC‑GGTGCTTCGC‑AAAGGTC‑AAATTCGCCCCCCGCGGAGTTTTCCCACCCGGGAGAGTGTTGACCCGCAAACAACAGATAAAAAAAAAAGcc                                                 >  1:118277/1‑147 (MQ=255)
                                                                                                      ttgatttgatGTCTGGCAGTTTATGGCGGGC‑GGCCTCCCCGCCCCCCT‑CCCGGGC‑GTTGCTTCGC‑CACGTTC‑TAAACCGCCCCCCGCGGGTTTTTCCTCCTCCGGGGGGGGGTTTCCCTCAAACAACAGAGAAAACaaaaaacc                                                 >  2:747326/1‑140 (MQ=2)
                                                                                                       tgatttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCt                                           <  2:1270203/149‑1 (MQ=40)
                                                                                                         atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTTCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAAACCCCTCCCCGCGGGTTTTTCCTTCTCCGGGGAGAGTTTACCGACAAACAAAAGATAAAAAAAAAAGGCCCGTCTTcc                                        >  1:458888/1‑150 (MQ=2)
                                                                                                         atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                       <  1:516253/151‑1 (MQ=38)
                                                                                                         atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCC‑CCGGGCC‑GTTTCTTCGCAAACGTTC‑‑AATCCCCCCCCGGCGGGATTTTCCTCCTCAGGAGAGTGTTTACCCACAAACAAACAATAAAAAAAAAAGACCAGTCTTcc                                        >  1:1618754/1‑150 (MQ=2)
                                                                                                            tgatGTCTGGCAGTTTATGGCGGGC‑GTCCTTCCCCCCACCCC‑CCCGGGC‑GGGGCTTCGC‑AACGTTC‑AAAACCGCCCCCGCCGGGGTTTTCCTATTCGGGGGAGGGTTCACCCACCAACAAAAGCTAAAAAAAAAAGccccagcttccccc                                     >  2:1517153/1‑139 (MQ=2)
                                                                                                                   tGGCAGTTTATGGCGGGC‑GTGCTGCCCCCCCCCCC‑CCCGGCC‑GGTTCTTTTC‑CACGTTC‑AAATTCGCCCCCCGCCGGATTGTTCTTCTTGGGGGAGGGGTTACCGACCAACAACAACAAAAAAAAAAAGCACCGCCTTCCGCCTGCGCCtt                             >  2:1090978/1‑151 (MQ=2)
                                                                                                                      cAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                          <  2:36922/151‑1 (MQ=37)
                                                                                                                           tATGGCGGGC‑GTCCTGCCCGCCCCCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGTATTTGTCCTACTCGGGAGAGTGTGCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTTATCCTTTTGtttt                      <  2:557555/150‑1 (MQ=255)
                                                                                                                           tATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                      <  1:1388401/150‑1 (MQ=38)
                                                                                                                               gcgggc‑gTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                 <  2:1400080/151‑1 (MQ=38)
                                                                                                                                tgggc‑gTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                 >  2:950476/2‑150 (MQ=38)
                                                                                                                                 gggc‑gTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt               <  1:803324/151‑1 (MQ=38)
                                                                                                                                      gTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc            <  1:1346141/150‑1 (MQ=40)
                                                                                                                                      gTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt           >  2:323276/1‑151 (MQ=40)
                                                                                                                                       tCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg          >  2:846037/1‑151 (MQ=40)
                                                                                                                                        ccTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg         <  2:624464/151‑1 (MQ=40)
                                                                                                                                           gcacGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGACCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTGGATGTCggg         <  2:334136/145‑1 (MQ=255)
                                                                                                                                           gCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc        <  1:1483683/149‑1 (MQ=40)
                                                                                                                                           gCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa       <  2:724860/150‑1 (MQ=40)
                                                                                                                                               gCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc  <  1:1155050/151‑1 (MQ=40)
                                                                                                                                               gCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc  >  1:1117816/1‑151 (MQ=40)
                                                                                                                                               |                                                                                                                                                          
ATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCC  >  NC_000913/3423145‑3423437

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