Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2940TC42.9% 153.6 / 42.0 112intergenic (+100/+129)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base T (17/47);  new base C (48/0);  total (65/47)
Fisher's exact test for biased strand distribution p-value = 1.59e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAC  >  NC_000913/3423147‑3423441
                                                                                                                                                    |                                                                                                                                                      
ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTcc                                                                                                                                                      <  1:132058/150‑1 (MQ=255)
ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTcc                                                                                                                                                      <  2:179952/150‑1 (MQ=255)
   aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCcgg                                                                                                                                                   >  1:446860/1‑150 (MQ=255)
   aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCA‑CCCCCCcgg                                                                                                                                                   >  1:357808/1‑150 (MQ=255)
    aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTGCTCCCCGCCA‑CCCTCCcgggc                                                                                                                                                 >  2:1049227/1‑149 (MQ=255)
    aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGcc                                                                                                                                                 <  2:997381/151‑1 (MQ=255)
     aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCA‑ccccccc                                                                                                                                                     >  3:165736/1‑145 (MQ=255)
        tgagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCA‑CCCCCCCGGcc                                                                                                                                                 >  3:57744/2‑147 (MQ=255)
         gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTGCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGt                                                                                                                                               <  2:120498/148‑1 (MQ=255)
                cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCC‑CCCCCCCGGCCGtt                                                                                                                                              >  3:408510/1‑142 (MQ=255)
                 gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCC‑CCCCCCCGGCCGGTGGtt                                                                                                                                          >  2:718725/1‑145 (MQ=37)
                      cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCA‑CCCCCCCCGCCCGTGCTTCGCaccaatgt                                                                                                                               >  2:855206/1‑144 (MQ=255)
                       aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCCGCCTGCCCGCCC‑CCCTCCCGGGCGGTTCTTTTCCAC‑GTT‑c                                                                                                                              >  3:233608/1‑149 (MQ=255)
                        tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGCCCGTGCTTCGCCAC‑GGT‑cca                                                                                                                            >  2:1072079/1‑148 (MQ=25)
                           ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCC‑CCCCCCCGGGCGGTGCTTTGCAAC‑GTT‑TCAATcc                                                                                                                        >  1:622190/1‑151 (MQ=25)
                            tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATcc                                                                                                                        <  2:750798/150‑1 (MQ=255)
                            tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCC‑CCCCCCGGGCCGTTGCTTCGCAAC‑GTT‑Caaa                                                                                                                           >  1:329611/1‑147 (MQ=255)
                            tATCGCAGGTAATCCATTAATGGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATcc                                                                                                                        <  2:366659/150‑1 (MQ=255)
                             aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGGCGTTGCTTCCCAAC‑GTT‑CAAATcccc                                                                                                                      >  1:256768/1‑149 (MQ=21)
                               cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCt                                                                                                                     <  1:1126530/150‑1 (MQ=255)
                                gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTcc                                                                                                                   <  3:199626/151‑1 (MQ=255)
                                gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTcc                                                                                                                   <  1:272451/151‑1 (MQ=255)
                                gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCGCCC‑CCCCCCCGGCCGCTGTTTCGCAAC‑GTT‑CAAAACCGCTcc                                                                                                                   >  4:146006/1‑151 (MQ=21)
                                 cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGCCGGTGCTTCGCCAC‑GTT‑CAAATCCCCTccc                                                                                                                  >  1:541116/1‑151 (MQ=21)
                                  aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCGCCCCCCCA‑CCCCCCCGGCCGCTGCTTTCCAAA‑GTT‑CAAAAccccccccc                                                                                                                 >  2:143438/1‑150 (MQ=21)
                                    gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGCCCCTGCTTTGCAAC‑CTT‑CCAATccccccccc                                                                                                                 >  2:260245/1‑148 (MQ=21)
                                    gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCA‑CCCTCCCGGCCCTTGCTTCGCAAC‑CAT‑CAAATCCGCTccc                                                                                                                  >  1:563374/1‑148 (MQ=21)
                                       aTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCCCCCGGCCGTTGCTTCGCAAC‑CTT‑CAAATCCCCTCCccgcgg                                                                                                             >  1:390335/1‑150 (MQ=255)
                                         ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCCCCC‑CCCCCCCGGCCGTTGCCTCGCCAC‑GTT‑CTAAACCGCTcccccgcgggt                                                                                                          >  2:497229/1‑143 (MQ=18)
                                           aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCCGGGCGTTTCTTTGCAAC‑CTT‑CAAAACCGCTCCcggcggggtt                                                                                                         >  4:454057/1‑146 (MQ=18)
                                           aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCC‑CCCCCCCGGCCGCTGCTTCTCCAC‑CTT‑CTAATTCTCTCCCGCCGGGttt                                                                                                         >  3:313078/1‑150 (MQ=255)
                                               aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                    <  4:412655/151‑1 (MQ=255)
                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCC‑CCCCCCCGGCCGGTGCTTTGCCAA‑GTT‑TAAAAACCCTCCccgcggggtttttc                                                                                                     >  1:188467/1‑140 (MQ=255)
                                                   aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCTCCCGGCCGTTTCTTTGCCAC‑CTT‑CAAATCCCCTCCCCGCGGGTTTTTCCTACTc                                                                                                >  3:168649/1‑151 (MQ=18)
                                                       ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGTCTGCCCCCCC‑CCCCCCCGGCCGTTGCTTTCCAAC‑GTT‑TAAATCCCCCCCCCGCGGGGTTGTcc                                                                                                     >  3:49694/1‑142 (MQ=255)
                                                        tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCA‑CCCCCCCGGCCCGTGCTTTGCCAA‑GTT‑TAAAAACGCCCCCCGCGGGATTGTcc                                                                                                     >  1:880983/1‑141 (MQ=17)
                                                        tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTCCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCGCCC‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGga                                                                                            <  3:21443/150‑1 (MQ=255)
                                                         aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCGCCCCCCC‑CCCCCCCGGCCGTCGCCTCTCCAC‑ATT‑TCAATACGCTCCCCGCGGGGTTGttcctccccgg                                                                                             >  2:1003514/1‑138 (MQ=255)
                                                         aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgaga                                                                                          <  2:1275401/151‑1 (MQ=255)
                                                           ttCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCC‑CCCCCCGCGCCCTTTCTTCTCAAA‑GCT‑TAAATCCCCCCCCCGCCGGTTTCTCCTtctctcg                                                                                             >  3:365919/1‑143 (MQ=255)
                                                             cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCC‑CCCCCCCGGCCGTTGCTTTGCAAC‑GTT‑CAAATCCGCCCCCCGCGGGTTTTTTCTACTCcgggg                                                                                           >  1:75640/1‑144 (MQ=17)
                                                              gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGtt                                                                                     <  2:218355/151‑1 (MQ=255)
                                                               aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCA‑CCCCCCCGGCCGTTGCTTCGCAACAGTT‑C‑AATCCCCCCCCCGCCGGGTTGTCCTCCTCCGGggaggg                                                                                       >  4:403436/1‑146 (MQ=17)
                                                                  aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCCCCCCCCC‑CCCCCCCGGCCGTTGCTTCGCAAA‑CTT‑TCAATCCGCTCCCCCCGGGTTTGTCCTACTCCGGCGGGGGTTCAcc                                                                                 >  2:229642/1‑151 (MQ=255)
                                                                   gCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCg                                                                                <  1:1224605/151‑1 (MQ=255)
                                                                    cAAAAAGGCCTTCCTTTCGGATGGCCTTTCGCTCGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCTCCGa                                                                               <  2:204466/151‑1 (MQ=255)
                                                                    cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGa                                                                               <  3:280501/151‑1 (MQ=255)
                                                                    cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCA‑CCCCCCCGGCCGTTGCTTCGCAAA‑GTT‑CAAAACCGCCCCCCGCGGGTTTTTCCTTCTCCGGGGAGTGTTTAcc                                                                                 >  4:108029/1‑149 (MQ=17)
                                                                     aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCCCCCGGCCGTTGCTTCGCCAA‑GTT‑CAAATCCCCTCCCCGCGGGTTTTTCCTACTCGGGAGAGGGTTTAcc                                                                                 >  4:290265/1‑148 (MQ=17)
                                                                      aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCA‑CCCCCCCGGGCGGTGCTTCGCAAC‑GCT‑CAAAACCCCCCCCCGCGGGTTTGTTCTTCTCcgggggggggt                                                                                     >  2:1171761/1‑135 (MQ=255)
                                                                        aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                           <  1:209576/151‑1 (MQ=255)
                                                                          ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac‑a                                                                        <  2:357748/151‑1 (MQ=255)
                                                                           gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCC‑CCCCCCCGGCCGTTGCTTTGCAAC‑CGT‑CAAAATCCCCCCCCGCGGGTTTTTCCTTCTCCGGGGAGGGTTCACCCacaaac‑a                                                                        >  1:840357/1‑150 (MQ=255)
                                                                                ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGa                                                                   <  2:975978/150‑1 (MQ=255)
                                                                                  tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGCCGTTTCTTCCCAAC‑GTT‑TAAATCCCCCCCCCGCGGGTTTTTCCTACTCGGGGGGGGGTTCTCCCACCAac‑aacaaaaa                                                                 >  2:1002501/1‑146 (MQ=9)
                                                                                    tCGGGTGGCCTGACGGTTGATTGGGGGTCTGGCAGTTTATGGTGGGGGTCCTGCCCGCCC‑CCCTCCGGGCCGTTGCTTCGCAGC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGAt                                                                  <  2:1259152/147‑1 (MQ=14)
                                                                                     cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATaaaa                                                              <  1:1266490/150‑1 (MQ=255)
                                                                                     cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATa                                                                 <  2:1233511/147‑1 (MQ=255)
                                                                                      ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACa                                                            <  3:423581/151‑1 (MQ=255)
                                                                                      ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCC‑CCCCCCGGGGCGGTGGTTCGCAAA‑GTT‑CAAAACCGCTCCCGGCGGGGTTTTCCTACTCGGGGGAGGGTTCACCGACAAAC‑AACAGATaaaa                                                              >  1:918582/1‑149 (MQ=14)
                                                                                      ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCC‑CCCCCCCGGGCCTTTCTTCTCAAA‑GTT‑CAAATCCCCCCCCCGCGGGGTTGTCCCACTCTGGGGAGGGTTTACCCACCaaa‑aaaaaaaaaaaa                                                             >  2:654333/1‑149 (MQ=255)
                                                                                       gATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACaa                                                           <  2:421768/151‑1 (MQ=255)
                                                                                        aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCCCCA‑CCCTCCCGGCCGTTGCTTCGCCAC‑CTT‑CAAATCCCCCCCCCGCGGGTTTGTTCCTCTCGGGAGAGGGTTTACCGCCaaaa‑aaaaaaaaaaaaaa                                                           >  1:1136823/1‑150 (MQ=9)
                                                                                           gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCT‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAg                                                        <  2:18207/150‑1 (MQ=255)
                                                                                            ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCACCCCCCCGGGCCG‑TGCTTCTCGAA‑CTT‑CAAATCCCCTCCCCGCGGATTTGTCCTTCTCTGGGGAGGGTTCACCCCCCAAA‑AAAAAATAAAACaaaa                                                         >  1:758783/1‑148 (MQ=255)
                                                                                            ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCC‑CCCCCCCGGCCGGTGCTTCGCAAA‑CTT‑TCAATCCGCCCCCCGCGGGCTTTTCCTTCTCGGGGGGGGGTTCTCCCACAAACAAACAaaaaaaaaaaaa                                                         >  2:1030548/1‑149 (MQ=9)
                                                                                            ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTCCCCCCCC‑CCCCCCCGGCCGGTGCTTTGCAAC‑CTT‑TAAATCCGCCCCCCGCGGGTTTGTCCTTCTCCGGGGAGGGGTTACCCACCAAC‑aaaaaaaaaaaaaaaaa                                                        >  1:261699/1‑148 (MQ=9)
                                                                                               ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGCCGGTTCTTCGCAAA‑GTT‑CAAATCCGCTCCCCGCGGGTTTTTCCTTCTCCGGAGAGGGTTTACCCCCAAAC‑AAAAGAAAAAAAAAAAAGccc                                                    >  1:592315/1‑150 (MQ=9)
                                                                                               ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCA‑CCCCCCCGGCCGTTGCTTCGCAAC‑GTTCCAAATCC‑CCCCCCGCGGGTTTTTCCTACTCCGGAGAGGGGTTACCGACAAAA‑AACAGATAAAACAAAAAGccc                                                    >  1:498564/1‑150 (MQ=14)
                                                                                                 cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGGCCCCCA‑CCCCCCCGGCCGCTGCCTCGCAAC‑ATT‑CAAAACCGCTCCCCGCCGGGTTTTCCTACTCCGGCGGGGGGTTACCGCCAAAC‑AAAAGATAaaaaaaaaagacc                                                    >  2:1082448/1‑143 (MQ=2)
                                                                                                     tgatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCt                                               <  1:897819/149‑1 (MQ=40)
                                                                                                      gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTcc                                            <  2:619219/151‑1 (MQ=39)
                                                                                                      gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTcc                                            <  1:716756/151‑1 (MQ=39)
                                                                                                      gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTcc                                            <  2:308865/151‑1 (MQ=39)
                                                                                                      gatttgatGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCA‑CCCCCCCGGCCGGTTCTTTGCAAC‑GTT‑CAAAATCGCCCCCGGCGGGTTTGTCCTACTCGGGGGAGTGTTTCCCCACCAAA‑AACAACAAAAAAAAAAAGCCCAGCCTTcc                                            >  2:968094/1‑151 (MQ=255)
                                                                                                       atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCg                                           <  4:219035/151‑1 (MQ=38)
                                                                                                        tttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCCCCCGGCCGGTGCCTCGCAAC‑GTT‑TAAAACCGCTCCCCGCGGGTTTTTCCTACTCGGGAGAAGGTTCACCCACCAAC‑AACAGGTAAAACAAAAAGcc                                                     >  1:554541/1‑140 (MQ=11)
                                                                                                            atGTCTGGCAGTTTATGGCGGGGGGCCCTCCCGCCC‑CCCCCCCGGCCGTTGTTTCTCAAA‑GTT‑TCAAACCCCTCCCGGCGGGTTTTTCCTTCTTCGGGGAGGGGTTACCGCCAAAC‑AACAGATAAAAAAAAAAACCCAGTCTTTCg                                           >  2:212797/1‑146 (MQ=255)
                                                                                                            atGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGa                                      <  2:546743/151‑1 (MQ=38)
                                                                                                            atGCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCC‑CCCCCCCGGCCCTTGCTTCGCCAC‑GTT‑CAAAACCGCCCCCCGCGGGTTTTTTCCACTCAAGAGAGCGGTCACCCCCAAAA‑AAAAAAAAAAAAAAAAAGGCCCGGCTTCCCAc                                         >  1:27714/1‑148 (MQ=0)
                                                                                                               tCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCC‑CCCCCCCGGCCGTTGATTCGCAAA‑GTT‑TCAATCCGTTCCCCGCGGGTTTTTTCTATTTCGGACGGTGTTTAACGACAAAA‑AACAAAAAAAACAAAAAGCCCAGTCTTCCgtgtgacgc                                   >  1:136390/1‑143 (MQ=2)
                                                                                                                  ggCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                 <  2:330381/150‑1 (MQ=35)
                                                                                                                    cAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGCCGTTGCTTTGCAAC‑CTT‑CCAAACCCCCCCCCGCGGGTTTTTCCTACTCGGGGGAGGGGTCACCGACAAAC‑AACAGAAAAAACAAAAAGACCAGTCTTTCCACTGCGCCtttt                               >  2:306704/1‑150 (MQ=2)
                                                                                                                      gTTTATGGCGGGCGTCCTGCCCGCCC‑CCCCCCGGGCCGGTGCTTTGCAAC‑GCT‑TTAATCCCCCCCCCGCGGGGTTTTTCTTCTCAGGGGAGCGTTTACCCCCAAAA‑AAAAAAAAAAAAAAAAAACCCCCGCTTTCTGCTGAGCCCTTTGt                             >  2:120643/1‑150 (MQ=0)
                                                                                                                      gTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                            <  1:253134/151‑1 (MQ=38)
                                                                                                                          aTGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt                        <  3:115381/151‑1 (MQ=38)
                                                                                                                            ggcgggcgTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                      >  2:1077025/1‑151 (MQ=38)
                                                                                                                              cgggcgTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                    <  3:18271/151‑1 (MQ=38)
                                                                                                                               gggcgTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                   <  2:502623/151‑1 (MQ=38)
                                                                                                                                     ccTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCta              <  2:1031657/150‑2 (MQ=255)
                                                                                                                                     ccTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg             >  3:88254/1‑151 (MQ=40)
                                                                                                                                       tGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc            >  4:443834/1‑150 (MQ=40)
                                                                                                                                        gCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc            >  3:232196/1‑149 (MQ=40)
                                                                                                                                        gCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc            >  4:435711/1‑149 (MQ=40)
                                                                                                                                        gCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc            <  2:1282703/149‑1 (MQ=40)
                                                                                                                                        gCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAg          <  4:47763/151‑1 (MQ=40)
                                                                                                                                            gacA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGACCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc      >  3:246117/3‑151 (MQ=255)
                                                                                                                                            gCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc      >  2:72164/1‑151 (MQ=40)
                                                                                                                                               a‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTGTGTTTTATTTGATGTCTGGCAGTTCCCTa   <  2:79767/151‑1 (MQ=255)
                                                                                                                                                 cccTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAc  <  1:678176/151‑1 (MQ=40)
                                                                                                                                                    |                                                                                                                                                      
CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTGCTTCGCAAC‑GTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAAC‑AACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAC  >  NC_000913/3423147‑3423441

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