Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2940TC48.6% 69.8 / 28.4 72intergenic (+100/+129)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base T (9/27);  new base C (33/1);  total (44/28)
Fisher's exact test for biased strand distribution p-value = 1.23e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCC  >  NC_000913/3423147‑3423438
                                                                                                                                                   |                                                                                                                                                    
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)
               aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGc                                                                                                                                    <  2:583623/150‑1 (MQ=255)
               aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCCCCCCCCCCCCCCCGCGCCGGTGCtt                                                                                                                                        >  1:288109/1‑147 (MQ=37)
                cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCGtt                                                                                                                                            >  2:26421/1‑142 (MQ=255)
                 gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCCCCCCCCCCGCGGCTGCTTCTCaa                                                                                                                                 >  2:991633/1‑150 (MQ=255)
                  aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCGGGCCGtttcctttcaa                                                                                                                                   >  2:396369/1‑140 (MQ=255)
                       aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCGCCACCCCCCCGGGCGTTGCTTCCCCACGtt                                                                                                                             >  1:546102/1‑148 (MQ=255)
                       aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCGGCCCCCCCCCAACCGGGGCggggga                                                                                                                                         >  1:1000527/1‑133 (MQ=25)
                        tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTTCTTCCCCCCCCCCCCCCGGCCGTTTCTTCTCAACCTTcta                                                                                                                          >  1:1035614/1‑148 (MQ=25)
                         ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCCGGCCGTTGCTTCGCAACCCTCaaa                                                                                                                         >  1:26654/1‑150 (MQ=25)
                          tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                        <  2:236969/150‑1 (MQ=255)
                          tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGGTTCTTCGCCACGGTCaaa                                                                                                                         >  1:486046/1‑149 (MQ=25)
                           ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTTCCCCCCCCCCCCCCGGCCGGTGCTTTGCCACGtttcaa                                                                                                                         >  2:374790/1‑144 (MQ=255)
                               cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTCCCCCCCCCCCCCCCGGGCGCTGCTTTTCAACGATCTAATCCGc                                                                                                                    >  1:22116/1‑149 (MQ=255)
                                 cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTTGCCGCCCCCCCCCGGGGCCTTGTTTCGCAACGGTTAAAttccccccc                                                                                                                >  2:100915/1‑143 (MQ=255)
                                  aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGTGTCCTGCCCGCCCCCCCCCGCGCCGTCGCTTCCCAAAGATCAAATcccccccc                                                                                                                >  1:114738/1‑150 (MQ=255)
                                     tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGTCCCCCCCCCCCCCGGCCCTTGCTTCGCCACGTTCAAAACCCCCCCccgc                                                                                                             >  1:858905/1‑150 (MQ=21)
                                     tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTTCTTCGCAACGTTCAAAACCTCCCCccgc                                                                                                             >  1:28308/1‑150 (MQ=21)
                                         ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                       <  2:785342/151‑1 (MQ=255)
                                          cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTTCTTTGCCACATTCCAATCCCCCCCCCGCCGGttt                                                                                                      >  2:534834/1‑151 (MQ=18)
                                               aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                 <  1:402058/151‑1 (MQ=255)
                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCTCCCCGGCGTTTCTTCTGCAAGATCCAAACCGCCCCccgcggggtttttct                                                                                                 >  1:352165/1‑140 (MQ=18)
                                                  gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGTTTCTTCTCAACCTTTAAATCCGCCCCCCGCGGGGTTTTCCTa                                                                                                >  1:829708/1‑149 (MQ=18)
                                                        tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTTGCAACGTTCAAAACCGCCCCCCGCGGGTTTTTTCTACTCcgg                                                                                          >  1:1108460/1‑149 (MQ=17)
                                                              gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGtt                                                                                 <  1:667920/151‑1 (MQ=255)
                                                               aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCCGGGCGGTTCTTCGCAACGCTCAAAACCGCCCCCCGCGGGTTTTTCCTCCTCCGGAGAGTGtt                                                                                 >  1:531552/1‑150 (MQ=17)
                                                                  aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCGGGCCGGTGCTTTGCAACGTTTAAAACCCCCCCCCGCGGGGTTTTCCTACCCCGGGGAGGGTTTAcc                                                                             >  1:157171/1‑151 (MQ=14)
                                                                       aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTGCTCCCCCCCCCCCCCCCGGCCGCTGCCTCCCAAACCTCAAAACCCCTCCCCGCGGGTTTTTTCTTCCCCGGAGGGGGTTCACCCac                                                                          >  1:1015384/1‑149 (MQ=255)
                                                                              aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAg                                                                 <  1:862650/151‑1 (MQ=255)
                                                                              aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGCTGCTTTCCAAAGGTCAAATCCCCCCCCCGCGGGATTTTCCTACTCCGGGGAGGGTTCCCCCacaaacaa                                                                    >  1:53739/1‑148 (MQ=9)
                                                                                ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGa                                                                <  2:845207/150‑1 (MQ=255)
                                                                                  tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCACCCCCCGGGCCGCTTCTTCGCAACCTTCAAATCCGCCCCCCGCGGGGTTTTCCTCCTCGCGGGAGGGGTCACCGacaa                                                                        >  1:779334/1‑140 (MQ=14)
                                                                                  gtcaggATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCGTTTCCTCCCAAACTTTTAAACCCCTCCCCGGGGGATTTTTCTACCCCCGAGGGGGGTCACCCacaaa                                                                       >  1:200950/5‑141 (MQ=0)
                                                                                   ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTTGCAACATTCAAAACCGCTCCCCGCGGGGTTTTTCTTCTCGGGGGGGGGTTTACCGACAAaaaacaaaaaa                                                             >  2:442156/1‑145 (MQ=255)
                                                                                      ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAc                                                          <  1:555303/150‑1 (MQ=255)
                                                                                      ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACa                                                         <  1:667981/151‑1 (MQ=255)
                                                                                         tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCGCCCGGCCCGTGTTTCCCCAATATTAATTCCCCCCCCCCCCGATTTGTTCTTCTCGCGGCGGGGTTTTCCCACCAACAACCGATaaaaaaaa                                                       >  2:166305/1‑150 (MQ=9)
                                                                                            ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCACCCCCCGGGCCGTTGCTTTGCCAAGTTCAAAACCCCCCCCCGCGGGGTTTTTCCTCTCGGGGGAGTGGTTCCCGACCAACAACAGATAaaaaaaaa                                                      >  1:216895/1‑148 (MQ=9)
                                                                                               ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGCCGTTGTTTCTCAACCTTCAAATCCGCTCCCCGCGGGTTTTTTCTCCTCGGGGGAGTGGTTACCGACAAACAACAGATAAAAAAAAAAGcc                                                  >  1:486736/1‑149 (MQ=14)
                                                                                               ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGcc                                                  <  1:148442/149‑1 (MQ=255)
                                                                                               ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCTCCCGGCCGGTGCTTCTCAACGTTCAAATCCGCTCCCCGCGGGGTTGTTCTTCTCCGGAGAGGGTTCACCGACAAACAACAGATAAAAAAAAAAGccc                                                 >  1:48902/1‑150 (MQ=14)
                                                                                                 cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                               <  1:314452/150‑1 (MQ=255)
                                                                                                 cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                               <  2:804130/150‑1 (MQ=255)
                                                                                                  gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCTCCCCGCCGCTTCTTTGCAACCTTCAAATCCCCCCCCCGCGGGTTTTTCCTACTCGCGGGAGTGTTCTCCGCCAAAAAAAAGAAAAAAAAAAAAGccccaac                                             >  1:634688/1‑147 (MQ=255)
                                                                                                    ttaatttgatGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCCGCCGGTGCTTTCCAACGTTCCAATCCCCTCCCCGCGGGGTTTTTCTTCTCCAGAGGGGGGTTCCCCCCCAAACAAAAAAAAAAAAAAAAGGccccacct                                            >  2:875560/4‑146 (MQ=0)
                                                                                                     tgatttgatGTCTGGCAGTTTATGGCGGGCGTCCCGTCCCCCCCCCCCCGGGCCGTTTGTTTGCCAAGTTTCAATTCGCGCCCCGCCGCTTTTTTCCTCTCCGGGGAGGGTTTACCCACAAAAAACAAAGAAAAAAAAAAGACCAgg                                              >  2:819625/1‑146 (MQ=2)
                                                                                                      gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGGTGCTTCGCAACGTTTAAATCCCCCCCCCGCGGGTTTTTCCTCCTCCGGGGAGTGTTCACCGCCAAACAACAGATAAAAAAAAAAGCCCAGTCTTcc                                         >  1:550781/1‑151 (MQ=11)
                                                                                                       atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCCGGGCGTTGCTTCGCAACGTTCAAAACCGCCCCCCGCGGGGTTGTCCTACTCCGGAGAGTGGTTACCGACAAACAACAGATAAAAAAAAAAGACCAGTCTTcc                                         >  1:549387/1‑150 (MQ=11)
                                                                                                        tttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                        <  2:367559/150‑1 (MQ=38)
                                                                                                          tgatTTCTGGGAGTTTATGGCGGGCCTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                      <  2:962306/150‑1 (MQ=255)
                                                                                                          tgatGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCGGGCCCGTGCTTTGCAAAGTTTCAATACCCCCCCCGCGGGATTTTTCTACTCCGGAGGGGGTTCACCCACAAACAAAAAATAAAAAAAAAAGccccagctt                                           >  1:1155593/1‑139 (MQ=2)
                                                                                                                 tGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGGCGTTGCCTCGCAACGTTCAAAATCGCCCCCGGCGGGGTTTTTCCCCTCGGGAGAGGGGTTTCCCACCAAAAACAGATAAAACAAAAAGCCCAGTCTTTCCGCCGAGCCtt                              >  1:202411/1‑151 (MQ=2)
                                                                                                                  ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                <  1:529954/148‑1 (MQ=35)
                                                                                                                    cAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                           <  2:367996/151‑1 (MQ=37)
                                                                                                                      gTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGt                          <  1:1160488/150‑1 (MQ=38)
                                                                                                                          aTGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                       <  2:831493/149‑1 (MQ=38)
                                                                                                                            ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                    <  2:2541/150‑1 (MQ=38)
                                                                                                                             gcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                  <  1:333076/151‑1 (MQ=38)
                                                                                                                              cgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCATTCTTCCGACTGAGCCTTTTGTTTTATTTg                  <  2:713749/150‑1 (MQ=38)
                                                                                                                               gggcgTCCTGCCCGCCCCCCCCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                <  1:484659/151‑1 (MQ=255)
                                                                                                                               gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                <  1:551811/151‑1 (MQ=38)
                                                                                                                                 gcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt              >  2:1066078/1‑151 (MQ=39)
                                                                                                                                    tCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg           >  1:28788/1‑151 (MQ=40)
                                                                                                                                     ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg           <  2:183991/150‑1 (MQ=40)
                                                                                                                                       tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc         >  1:608085/1‑150 (MQ=40)
                                                                                                                                         cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGACTTTTGTTTTATTTGATGCCTGGc         <  2:31342/148‑1 (MQ=25)
                                                                                                                                             ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc  <  1:633273/151‑1 (MQ=40)
                                                                                                                                                   |                                                                                                                                                    
CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCC  >  NC_000913/3423147‑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: