Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2970GC39.9% 263.9 / 72.5 143intergenic (+103/+126)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (22/64);  new base C (57/0);  total (79/64)
Fisher's exact test for biased strand distribution p-value = 5.29e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCG  >  NC_000913/3423147‑3423446
                                                                                                                                                         |                                                                                                                                                          
ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCg                                                                                                                                                            <  1:1271412/151‑1 (MQ=255)
   aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCCCCC‑CCCggg                                                                                                                                                          >  1:253211/1‑150 (MQ=255)
     aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCCCCACCC‑CCCcgg                                                                                                                                                          >  3:487955/1‑148 (MQ=255)
      tCAGAGCGTACGCGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGAC‑GCCCT‑CCCGCCACCC‑CCCCGGcc                                                                                                                                                        >  3:271929/1‑149 (MQ=255)
       cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTGCTGCCCCCCCCCC‑CCCCGGcc                                                                                                                                                        >  2:1219963/1‑148 (MQ=255)
        agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCCGCCCGCCCCCC‑CCCGGGc                                                                                                                                                         >  1:1362909/1‑146 (MQ=255)
               aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCCCCCCCC‑CCCCCGCCGGTGCtt                                                                                                                                                 >  2:1381000/1‑147 (MQ=37)
                cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCaa                                                                                                                                            <  3:576764/151‑1 (MQ=255)
                 gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑c                                                                                                                                          <  2:879721/151‑1 (MQ=255)
                 gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCCCCCCCC‑CCCCGGCCGTTGCTTCGCaa                                                                                                                                            >  3:172826/1‑150 (MQ=255)
                  aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑G‑GGCCTTCCCCCCCCCC‑CCCccgccggccctt                                                                                                                                                 >  1:258983/1‑138 (MQ=255)
                   gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑G‑GGCCCGCCCGCCCCCC‑CCCCGGGCGGTGCTTCTCaa                                                                                                                                            >  1:1249997/1‑148 (MQ=255)
                     gCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGTCTGCCCCCCCCCC‑CCCCGGCCGCTGTTTCGCCA‑A‑Gtt                                                                                                                                      >  1:1553684/1‑150 (MQ=255)
                     gCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCTGCCCGCCCCCC‑CCCCGGCCGTTGCCTCGCaa                                                                                                                                            >  2:1026256/1‑146 (MQ=255)
                      cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCCCCC‑CCCCCGGCGGTTCGTCGCCA‑C‑Gtt                                                                                                                                      >  1:1695446/1‑149 (MQ=25)
                       aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCaa                                                                                                                                   <  2:1531992/151‑1 (MQ=255)
                       aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCa                                                                                                                                    <  3:108898/150‑1 (MQ=255)
                           ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCTGCCCGCCACCC‑TCCCGGCCGGTGCTTCTCCA‑C‑GGTTaaa                                                                                                                                  >  1:996255/1‑148 (MQ=255)
                              tCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTCCCCCCCCCCC‑CCCCGGCCGCTGCTTCGCAA‑C‑GGTCAAAAcccc                                                                                                                             >  1:54219/1‑148 (MQ=21)
                              tCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGTCTGCCCCCCCCCC‑CCCCGGCCCGTGGTTCTCCA‑C‑ATTCaaa                                                                                                                                  >  1:30580/1‑145 (MQ=255)
                                gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCTTCCCCCCACCC‑CCCCCGCCGTTGCTTCGCAA‑C‑GTTCAAAACCCCTcc                                                                                                                          >  2:568724/1‑151 (MQ=21)
                                 cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTcca                                                                                                                         <  1:20986/151‑2 (MQ=255)
                                     tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCTTCCCCCCCCCC‑CCCCGGCCGGTGCTTCTCAA‑A‑GGTCaaaaacccccccc                                                                                                                        >  2:1670754/1‑138 (MQ=255)
                                       aTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCCCCC‑TCCCGGCCGCTGCTTCGCAA‑A‑GTTTAAATCCGCTCCccgc                                                                                                                      >  3:283629/1‑148 (MQ=21)
                                       aTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGa                                                                                                                   <  2:1380672/151‑1 (MQ=255)
                                          cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑CCCCGGCCGGTGCTTCGCCA‑C‑GTTCAAATCCCCTCCCCGCGGGttt                                                                                                                >  2:1191502/1‑151 (MQ=18)
                                             tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCCCCACCC‑CCCCGGCCCGTGCTTTGCCA‑A‑GTTTAAATCCCCCCCccgccgggttgt                                                                                                              >  2:927746/1‑142 (MQ=18)
                                                ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCTTCCCCCCACCC‑CCCCGGCCGATGCTTCTCAA‑C‑ATTCAAATCCGCCCCCGGCGGGTTTGt                                                                                                              >  2:1436368/1‑147 (MQ=18)
                                                  gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTCCCCCCCCCCC‑CCCCGGCCGTTGCTTCGCAA‑C‑GTTCAAAACCCCTCCCCGCGGGTTTTTCCTAc                                                                                                         >  2:95563/1‑150 (MQ=18)
                                                   aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                       <  3:360167/151‑1 (MQ=255)
                                                   aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GCCCTGCCCGCCACCCGTCCGGG‑CGTTGCCTCGCAA‑C‑AGTGCAATACGCTCCCGGCGCGCTTGTCCttc                                                                                                         >  1:231968/1‑147 (MQ=255)
                                                    aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg                                                                                                      <  2:189834/151‑1 (MQ=255)
                                                      gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCggg                                                                                                    <  4:434165/151‑1 (MQ=255)
                                                        tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGga                                                                                                   <  1:857028/150‑1 (MQ=255)
                                                         aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgag                                                                                                  <  2:383777/150‑1 (MQ=255)
                                                          gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑G‑GGCCTGCCCCCCCCCC‑CCCCGGCCGTTGCTTTGCCA‑C‑CGTCAAAACCCCTCCCCGCGGGTTTGTTCTACTCCGga                                                                                                   >  2:838825/1‑148 (MQ=255)
                                                          gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCCCCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGaga‑g                                                                                               <  4:579029/151‑1 (MQ=255)
                                                          gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCTGCCCCCCCCCC‑CCCCGGCCGGTGCTTTGCAA‑C‑CATTCAATTCGCTCCCCCCGGAttttttccact                                                                                                        >  2:1532450/1‑135 (MQ=17)
                                                             cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTTCTGCCCGCCACCC‑CCCGGGGCGTTGCTTTGCAA‑C‑GTTTTAAACCGCCCCCCGCGGGTTTGTTCTACTCGGGAGA‑gtg                                                                                             >  3:519425/1‑150 (MQ=17)
                                                             cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑gtg‑t                                                                                           <  2:1147528/151‑1 (MQ=255)
                                                             cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑gtg‑t                                                                                           <  3:363982/151‑1 (MQ=255)
                                                             cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑gtg                                                                                             <  1:1660441/150‑1 (MQ=255)
                                                              gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTCCCCCCCCCCC‑CCCCGGGCGTTGGCTTTCAA‑C‑CTTCAAAACCGCCCCCCGGGGAgtttttccattccgga                                                                                                   >  2:1169611/1‑131 (MQ=12)
                                                               aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCCGCCCCCCCCCC‑CCCCGGCCGTTTGTTCTCAA‑A‑GTTTAAAACCGACCCCCCCGGGTTTTGCCTACCCGGGcga‑g                                                                                               >  1:1549928/1‑146 (MQ=255)
                                                                aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCa                                                                                        <  2:1485018/151‑1 (MQ=255)
                                                                   gCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCCCCC‑CCCCGGCCCTTGCTTCGCAA‑C‑GTTCAAATCCCCTCCCCGCGGGTTTTTTCTACTCGGGGGG‑GGG‑TTCAcc                                                                                      >  1:759660/1‑150 (MQ=17)
                                                                    cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGa                                                                                    <  2:1250142/151‑1 (MQ=255)
                                                                     aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCTGCCCCCCACCC‑TCCCGGCCGTTGCTTCGCAA‑C‑ATTCAAATTCGCTCCCCGCGGGGTTGTCCCACTCCGGGGA‑GGG‑TTTAcc                                                                                      >  2:1453685/1‑148 (MQ=255)
                                                                       aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGacaa                                                                                 <  1:1235882/151‑1 (MQ=255)
                                                                        aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCCGCCCCCCCCCC‑CCCCGGGCGGTGCTTTGCCA‑C‑GTTCAAATCCCCCCCCCGCGGGTTTTTCCTTCTCcgggggggggttccc                                                                                        >  4:24169/1‑133 (MQ=9)
                                                                        aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCCGCCCCCCCCCC‑CCCCCGCCGTTGCCTCGCAA‑C‑ATTCAAATCCCCCCCCCGCGGGGTTTTCCTTCTCCGGgga‑ggggtttcccacaa                                                                                  >  1:527984/1‑137 (MQ=9)
                                                                         aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGT‑T‑CTCCCGCCCGCCACCC‑TCCGGGCCTATTCTTTGCAA‑A‑GGTCAAAACCCCCCCCCGCCGGATTGTCCTACTCGGGCGA‑GAG‑TTcgac                                                                                      >  2:565926/1‑141 (MQ=255)
                                                                             cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAacaaca                                                                           <  2:1016496/151‑1 (MQ=255)
                                                                                ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCCCCC‑CCCCGGCCGTTGCTTCGCAA‑C‑GTTCAAATTCACCCCCCGCGGGTTTGTCCTACTCCGGAGA‑GTG‑TTCACCCacaaaaaaaa                                                                           >  2:496903/1‑146 (MQ=14)
                                                                                 cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTGCTGCCCCCCCCCC‑CCCCCGCCGTTTCTTCGCAA‑C‑GTTCAAATCCCCTCCCCGCCGGTTTGTCTTCCCCGGGGGA‑GGG‑TTCACCCCCAaaccaccaac                                                                        >  1:989358/1‑146 (MQ=9)
                                                                                 cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCCCCCCCC‑CCCCGGCCCTTGCTTTTCAA‑C‑GTTTCAATCCCCCCCCCGCGGGGTTGTTCTAATCCGGGGG‑GGG‑GTTTCCGCCAAACAAAAGac                                                                        >  4:572811/1‑149 (MQ=9)
                                                                                    tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑G‑GTCCTGGCCGCCCCCC‑CCCCGGCCGTTTTTTCGCCA‑C‑GCTCAAAATCGCCCCCCGCGGGTTTGTTCCTCTCCGGGGG‑GTG‑GTTACCTCCCAACAACACAa‑aaa                                                                    >  1:288369/1‑150 (MQ=255)
                                                                                    tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑aaa                                                                    <  1:1477148/150‑1 (MQ=255)
                                                                                      ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACa                                                                 <  2:1215022/151‑1 (MQ=255)
                                                                                        aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACaaa                                                               <  2:466164/151‑1 (MQ=255)
                                                                                         tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACaaa                                                               <  1:927922/150‑1 (MQ=255)
                                                                                         tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCCCCCCCC‑CCCCGGCCGGTGCTTTGCAA‑A‑GTTCAAATCCGCCCCCCGCGGGTTTGTCCTACTCCGGAGA‑GGG‑GTTACCCACAAAAAAAAGAT‑AAAACaaa                                                               >  3:501071/1‑150 (MQ=14)
                                                                                          ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAg                                                             <  1:731846/151‑1 (MQ=255)
                                                                                            ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCTGCCCGCCACCC‑TCCCGGCCGGTGGTTCGCCC‑C‑GCTTAAAACCGCCCCCCCCGGGGTTGTCCTACTCCGGGGAGGGG‑TTACCCCAAAAACAA‑AGAT‑AAAACaaaaa                                                             >  1:1454207/1‑148 (MQ=9)
                                                                                              tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCCCCCCCC‑CCCCGGCCGGTGCTTTGCAA‑A‑GTTCAAAACCGCCCCCCGCGGGTTTGTCCTACTCGGGAGA‑GGG‑TTCACCGACAAACAAAAGAT‑AAAACAAAAAGccc                                                         >  2:1439711/1‑151 (MQ=14)
                                                                                                tCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGCCTCCCCCCCACCC‑ACCCGGCCGCTTCGTCCCCA‑CGGTCCAAACCCCCTCCC‑GCGGGTTTGTGCTTCTCTGTAGA‑GAGATTCA‑CGACAAACAAAAAAT‑Aaaaaaaaaacacc                                                         >  2:1322396/1‑144 (MQ=2)
                                                                                                 cGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GGTCCGCCCCCCCCCC‑TCCCGGCCGCTGCTTTTCAACC‑GTTC‑AATCCGCCCCCCGCGGGGTTGTCCCACTCCGGGGG‑GTG‑TTCTCCGCCAAAAAACAGAT‑Aaaaaaaaaagacc                                                         >  1:361850/1‑143 (MQ=255)
                                                                                                  gCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTc                                                     <  2:253254/151‑1 (MQ=40)
                                                                                                  gCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑CCCCGGCCGGTTCTTTCCAA‑A‑GGTTCAATTCGCCCCCCGCGGGCTTGTCCTACTCTGGGGG‑GGG‑GTTACCCACACACAACAAAT‑Aaaaaaaaaaaccc                                                         >  1:627156/1‑142 (MQ=2)
                                                                                                   cTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCCGCCCCCCCCCC‑TCCCGGCCGGTGCTTCGCAA‑A‑GTTCAAATCCCCTCCCCGCGGGGTTGTCCTACTCCGGAGA‑GGG‑TTTACCCCCAAAAAACAGAT‑AAAACAAAAAGCCCAGTCt                                                    >  3:430016/1‑151 (MQ=11)
                                                                                                      gatttgatGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTc                                                  <  1:1208384/150‑1 (MQ=39)
                                                                                                        tttgatGTCTGGCAGTTTATGGCGGG‑G‑GGGCTGCCCCCCACCC‑TCCCGGCCGGTGCTTTGCAA‑C‑CTTTAAATCCGCTCCCCGCGGAGTTTTCCTACTCGGGGGA‑GTG‑GTTACCCACAAACAACAGAA‑AAAACAAAAAGGCCACGCTTcc                                                 >  4:220797/1‑149 (MQ=255)
                                                                                                          tgatGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGAc                                              <  2:92743/150‑1 (MQ=38)
                                                                                                           gatGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTg                                            <  3:230405/151‑1 (MQ=38)
                                                                                                             tGTCTGGCAGTTTATGGCGGG‑CGGCCCTCCCCCCCCCCC‑CCCCGGCCGTGGTTTGCCAA‑C‑CTCCAAAACC‑CCCCCCGCGGGCTTGTTCCTCCCCGGAGA‑GGG‑GTCTCCCCCCAACAAAAGATAAAAAAAAAAGGCCCAGCCTTc                                                  >  4:405043/1‑144 (MQ=255)
                                                                                                              gTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAg                                          <  2:28296/150‑1 (MQ=38)
                                                                                                              gTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGc                                         <  3:468045/151‑1 (MQ=37)
                                                                                                                  ggCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt                                     <  4:76005/151‑1 (MQ=35)
                                                                                                                  ggCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCA‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                       <  2:232006/149‑1 (MQ=25)
                                                                                                                          aTGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                              <  1:300732/150‑1 (MQ=38)
                                                                                                                           tGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                            >  2:665928/1‑151 (MQ=38)
                                                                                                                            ggcggg‑c‑gTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                           <  1:1174620/151‑1 (MQ=38)
                                                                                                                            ggcggg‑c‑gTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                           >  1:10243/1‑151 (MQ=38)
                                                                                                                            ggcggg‑c‑gTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                           <  2:1733514/151‑1 (MQ=38)
                                                                                                                            ggcggg‑c‑gTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                            <  4:70483/150‑1 (MQ=38)
                                                                                                                              cggg‑c‑gTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                          >  2:1629692/1‑150 (MQ=38)
                                                                                                                                 g‑c‑gTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt                      >  2:1425402/1‑151 (MQ=39)
                                                                                                                                       ccTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                  <  2:682236/151‑1 (MQ=40)
                                                                                                                                       ccTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                  >  3:132042/1‑151 (MQ=40)
                                                                                                                                          gCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                 <  2:817323/149‑1 (MQ=40)
                                                                                                                                          gCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                 >  2:1292737/1‑149 (MQ=40)
                                                                                                                                          gCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                 >  1:255601/1‑149 (MQ=40)
                                                                                                                                          gCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa                <  2:1347233/150‑1 (MQ=40)
                                                                                                                                             cGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc            <  1:343093/151‑1 (MQ=40)
                                                                                                                                               ccACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc           <  1:565843/150‑1 (MQ=40)
                                                                                                                                                cACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTTCTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt         <  2:1114745/151‑1 (MQ=255)
                                                                                                                                                    c‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTATTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActc     >  1:1337015/1‑151 (MQ=255)
                                                                                                                                                      tCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct    >  2:871947/1‑151 (MQ=40)
                                                                                                                                                      tCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActc     >  1:650939/1‑150 (MQ=40)
                                                                                                                                                        cGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCg  >  1:1262986/1‑151 (MQ=40)
                                                                                                                                                         |                                                                                                                                                          
CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGG‑C‑GTCCTGCCCGCCACCC‑TCCGGGCCGTTGCTTCGCAA‑C‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTG‑TTCACCGACAAACAACAGAT‑AAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCG  >  NC_000913/3423147‑3423446

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