Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2970GC41.4% 83.0 / 36.4 70intergenic (+103/+126)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (14/27);  new base C (29/0);  total (43/27)
Fisher's exact test for biased strand distribution p-value = 2.01e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTC  >  NC_000913/3423147‑3423445
                                                                                                                                                      |                                                                                                                                                        
ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                                                                                          <  2:49117/151‑1 (MQ=255)
   aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCGCCCCCCCCCcgg                                                                                                                                                        >  2:551791/1‑150 (MQ=255)
   aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCggg                                                                                                                                                        <  1:359449/150‑1 (MQ=255)
       cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGGCTGCCCCCCCCCCTCCCGGCCGt                                                                                                                                                    >  2:1093235/1‑150 (MQ=255)
        agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTGGTCCTGCCCCCCCCCCCCCCGGcc                                                                                                                                                      >  1:662483/1‑147 (MQ=255)
               aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTTCCCCCCCCCCCCCCGGCCGGTGGTTcccc                                                                                                                                           >  2:862598/1‑148 (MQ=255)
               aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑Gc                                                                                                                                           <  2:583623/150‑1 (MQ=255)
               aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCCCCCCCCCCCCCCCGCGCCGGTGCtt                                                                                                                                               >  1:288109/1‑147 (MQ=37)
                cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCGtt                                                                                                                                                   >  2:26421/1‑142 (MQ=255)
                 gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCCCCCCCCCCGCGGCTGCTTC‑TC‑aa                                                                                                                                        >  2:991633/1‑150 (MQ=255)
                  aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCGGGCCGtttcctttcaa                                                                                                                                          >  2:396369/1‑140 (MQ=255)
                       aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCGCCACCCCCCCGGGCGTTGCTTC‑CC‑CACGtt                                                                                                                                    >  1:546102/1‑148 (MQ=255)
                       aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCGGCCCCCCCCCAACCGGGGCggggga                                                                                                                                                >  1:1000527/1‑133 (MQ=25)
                        tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTTCTTCCCCCCCCCCCCCCGGCCGTTTCTTC‑TC‑AACCTTcta                                                                                                                                 >  1:1035614/1‑148 (MQ=25)
                         ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCCGGCCGTTGCTTC‑GC‑AACCCTCaaa                                                                                                                                >  1:26654/1‑150 (MQ=25)
                          tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAAt                                                                                                                               <  2:236969/150‑1 (MQ=255)
                          tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGGTTCTTC‑GC‑CACGGTCaaa                                                                                                                                >  1:486046/1‑149 (MQ=25)
                           ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTTCCCCCCCCCCCCCCGGCCGGTGCTTT‑GC‑CACGtttcaa                                                                                                                                >  2:374790/1‑144 (MQ=255)
                               cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTCCCCCCCCCCCCCCCGGGCGCTGCTTT‑TC‑AACGATCTAATCCGc                                                                                                                           >  1:22116/1‑149 (MQ=255)
                                 cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTTGCCGCCCCCCCCCGGGGCCTTGTTTC‑GC‑AACGGTTAAAttccccccc                                                                                                                       >  2:100915/1‑143 (MQ=255)
                                  aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGTGTCCTGCCCGCCCCCCCCCGCGCCGTCGCTTC‑CC‑AAAGATCAAATcccccccc                                                                                                                       >  1:114738/1‑150 (MQ=255)
                                     tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGTCCCCCCCCCCCCCGGCCCTTGCTTC‑GC‑CACGTTCAAAACCCCCCCccgc                                                                                                                    >  1:858905/1‑150 (MQ=21)
                                     tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTTCTTC‑GC‑AACGTTCAAAACCTCCCCccgc                                                                                                                    >  1:28308/1‑150 (MQ=21)
                                         ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGAtt                                                                                                              <  2:785342/151‑1 (MQ=255)
                                          cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTTCTTT‑GC‑CACATTCCAATCCCCCCCCCGC‑CGGttt                                                                                                             >  2:534834/1‑151 (MQ=18)
                                               aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCt                                                                                                        <  1:402058/151‑1 (MQ=255)
                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCTCCCCGGCGTTTCTTCTGC‑AA‑GATCCAAACCGCCCCccgc‑ggggtttttct                                                                                                        >  1:352165/1‑140 (MQ=18)
                                                  gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGTTTCTTC‑TC‑AACCTTTAAATCCGCCCCCCGC‑GGGGTTTTCCTa                                                                                                       >  1:829708/1‑149 (MQ=18)
                                                        tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTT‑GC‑AACGTTCAAAACCGCCCCCCGC‑GGGTTTTTTCTACTCcgg                                                                                                 >  1:1108460/1‑149 (MQ=17)
                                                              gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGtt                                                                                        <  1:667920/151‑1 (MQ=255)
                                                               aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCCGGGCGGTTCTTC‑GC‑AACGCTCAAAACCGCCCCCCGC‑GGGTTTTTCCTCCTCCGGA‑GAGTGtt                                                                                        >  1:531552/1‑150 (MQ=17)
                                                                  aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCGGGCCGGTGCTTT‑GC‑AACGTTTAAAACCCCCCCCCGC‑GGGGTTTTCCTACCCCGGG‑GAGGGTTTAcc                                                                                    >  1:157171/1‑151 (MQ=14)
                                                                       aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTGCTCCCCCCCCCCCCCCCGGCCGCTGCCTC‑CC‑AAACCTCAAAACCCCTCCCCGC‑GGGTTTTTTCTTCCCCGGA‑GGGGGTTCACCCac                                                                                 >  1:1015384/1‑149 (MQ=255)
                                                                              aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAg                                                                        <  1:862650/151‑1 (MQ=255)
                                                                              aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGCTGCTTT‑CC‑AAAGGTCAAATCCCCCCCCCGC‑GGGATTTTCCTACTCCGGG‑GAGGGTTCCCCCacaaacaa                                                                           >  1:53739/1‑148 (MQ=9)
                                                                                ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGa                                                                       <  2:845207/150‑1 (MQ=255)
                                                                                  tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCACCCCCCGGGCCGCTTCTTC‑GC‑AACCTTCAAATCCGCCCCCCGC‑GGGGTTTTCCTCCTCGCGG‑GAGGGGTCACCGacaa                                                                               >  1:779334/1‑140 (MQ=14)
                                                                                  gtcaggATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCGTTTCCTC‑CCAAACTTTTAAACCCCTCCCCGGG‑GGATTT‑TTCTACCCCCGA‑GGGGGGTCACCCacaaa                                                                              >  1:200950/5‑141 (MQ=0)
                                                                                   ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTT‑GC‑AACATTCAAAACCGCTCCCCGC‑GGGGTTTTTCTTCTCGGGG‑GGGGGTTTACCGACAAaaaacaaaaaa                                                                    >  2:442156/1‑145 (MQ=255)
                                                                                      ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAAc                                                                 <  1:555303/150‑1 (MQ=255)
                                                                                      ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACa                                                                <  1:667981/151‑1 (MQ=255)
                                                                                         tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCGCCCGGCCCGTGTTTC‑CC‑CAATATTAATTCCCCCCCCCCC‑CGATTTGTTCTTCTCGCGG‑CGGGGTTTTCCCACCAACAACCGATaaaaaaaa                                                              >  2:166305/1‑150 (MQ=9)
                                                                                            ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCACCCCCCGGGCCGTTGCTTT‑GC‑CAAGTTCAAAACCCCCCCCCGCGGGGTTTTTCCT‑CTCGGGG‑GAGTGGTTCCCGACCAACAACAGATAaaaaaaaa                                                             >  1:216895/1‑148 (MQ=9)
                                                                                               ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGCCGTTGTTTC‑TC‑AACCTTCAAATCCGCTCCCCGC‑GGGTTTTTTCTCCTCGGGG‑GAGTGGTTACCGACAAACAACAGATAAAAAAAAAAGcc                                                         >  1:486736/1‑149 (MQ=14)
                                                                                               ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGcc                                                         <  1:148442/149‑1 (MQ=255)
                                                                                               ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCTCCCGGCCGGTGCTTC‑TC‑AACGTTCAAATCCGCTCCCCGC‑GGGGTTGTTCTTCTCCGGA‑GAGGGTTCACCGACAAACAACAGATAAAAAAAAAAGccc                                                        >  1:48902/1‑150 (MQ=14)
                                                                                                 cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                      <  2:804130/150‑1 (MQ=255)
                                                                                                 cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                      <  1:314452/150‑1 (MQ=255)
                                                                                                  gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCTCCCCGCCGCTTCTTT‑GC‑AACCTTCAAATCCCCCCCCCGC‑GGGTTTTTCCTACTCGCGG‑GAGTGTTCTCCGCCAAAAAAAAGAAAAAAAAAAAAGccccaac                                                    >  1:634688/1‑147 (MQ=255)
                                                                                                    ttaatttgatGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCCGCCGGTGCTTT‑CC‑AACGTTCCAATCCCCTCCCCGC‑GGGGTTTTTCTTCTCCAGAGGGGGGTTCCCCCCCAAACAAAAAAAAAAAAAAAAGGccccacct                                                   >  2:875560/4‑146 (MQ=0)
                                                                                                     tgatttgatGTCTGGCAGTTTATGGCGGGCGTCCCGTCCCCCCCCCCCCGGGCCGTTTGTTT‑GC‑CAAGTTTCAATTCGCGCCCCGC‑CGCTTTTTTCCTCTCCGGG‑GAGGGTTTACCCACAAAAAACAAAGAAAAAAAAAAGACCAgg                                                     >  2:819625/1‑146 (MQ=2)
                                                                                                      gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGGTGCTTC‑GC‑AACGTTTAAATCCCCCCCCCGC‑GGGTTTTTCCTCCTCCGGG‑GAGTGTTCACCGCCAAACAACAGATAAAAAAAAAAGCCCAGTCTTcc                                                >  1:550781/1‑151 (MQ=11)
                                                                                                       atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCCGGGCGTTGCTTC‑GC‑AACGTTCAAAACCGCCCCCCGC‑GGGGTTGTCCTACTCCGGA‑GAGTGGTTACCGACAAACAACAGATAAAAAAAAAAGACCAGTCTTcc                                                >  1:549387/1‑150 (MQ=11)
                                                                                                        tttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                               <  2:367559/150‑1 (MQ=38)
                                                                                                          tgatTTCTGGGAGTTTATGGCGGGCCTCCGGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                             <  2:962306/150‑1 (MQ=255)
                                                                                                          tgatGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCGGGCCCGTGCTTT‑GC‑AAAGTTTCAATACCCCCCCCGC‑GGGATTTTTCTACTCCGGA‑GGGGGTTCACCCACAAACAAAAAATAAAAAAAAAAGccccagctt                                                  >  1:1155593/1‑139 (MQ=2)
                                                                                                                 tGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGGCGTTGCCTC‑GC‑AACGTTCAAAATCGCCCCCGGC‑GGGGTTTTTCCCCTCGGGA‑GAGGGGTTTCCCACCAAAAACAGATAAAACAAAAAGCCCAGTCTTTCCGCCGAGCCtt                                     >  1:202411/1‑151 (MQ=2)
                                                                                                                  ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                       <  1:529954/148‑1 (MQ=35)
                                                                                                                    cAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                                  <  2:367996/151‑1 (MQ=37)
                                                                                                                      gTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGt                                 <  1:1160488/150‑1 (MQ=38)
                                                                                                                          aTGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                              <  2:831493/149‑1 (MQ=38)
                                                                                                                            ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                           <  2:2541/150‑1 (MQ=38)
                                                                                                                             gcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                         <  1:333076/151‑1 (MQ=38)
                                                                                                                              cgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCATTCTTCCGACTGAGCCTTTTGTTTTATTTg                         <  2:713749/150‑1 (MQ=38)
                                                                                                                               gggcgTCCTGCCCGCCCCCCCCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                       <  1:484659/151‑1 (MQ=255)
                                                                                                                               gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                       <  1:551811/151‑1 (MQ=38)
                                                                                                                                 gcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt                     >  2:1066078/1‑151 (MQ=39)
                                                                                                                                    tCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                  >  1:28788/1‑151 (MQ=40)
                                                                                                                                     ccTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                  <  2:183991/150‑1 (MQ=40)
                                                                                                                                       tGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                >  1:608085/1‑150 (MQ=40)
                                                                                                                                         cccGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGACTTTTGTTTTATTTGATGCCTGGc                <  2:31342/148‑1 (MQ=25)
                                                                                                                                             ccACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc         <  1:633273/151‑1 (MQ=40)
                                                                                                                                                    ccGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACtctc  <  2:924182/151‑1 (MQ=40)
                                                                                                                                                      |                                                                                                                                                        
CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTC  >  NC_000913/3423147‑3423445

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