Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2870GC41.1% 76.9 / 34.9 73intergenic (+93/+136)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (15/28);  new base C (30/0);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 1.85e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGC  >  NC_000913/3423141‑3423431
                                                                                                                                                  |                                                                                                                                                    
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)
                     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)
                                                                                                                                                  |                                                                                                                                                    
GCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGC  >  NC_000913/3423141‑3423431

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: