Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2900AC43.8% 50.8 / 29.3 73intergenic (+96/+133)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base A (15/26);  new base C (30/2);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 4.47e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCC  >  NC_000913/3423141‑3423438
                                                                                                                                                     |                                                                                                                                                        
gcgcATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCa                                                                                                                                                          <  1:288928/150‑1 (MQ=255)
 cgcATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTGCCTCCCCCCAcc                                                                                                                                                        >  2:724345/1‑151 (MQ=255)
      ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                                                                                   <  2:49117/151‑1 (MQ=255)
      ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCCaccccccacccccc                                                                                                                                                    >  1:922915/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                        
GCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTC‑GC‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGA‑GAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCC  >  NC_000913/3423141‑3423438

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: ‑