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

AGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCT  >  NC_000913/3423155‑3423444
                                                                                                                                               |                                                                                                                                                      
agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCTCCCGGcc                                                                                                                                                    >  2:663232/1‑147 (MQ=255)
    cGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                           <  1:116897/151‑1 (MQ=255)
         gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTCGTCCTCGCCGCCACCCTCCCGCCCGGTGATTccc                                                                                                                                         >  1:791964/1‑146 (MQ=255)
         gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCTCaa                                                                                                                                      >  2:97787/1‑150 (MQ=255)
            ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCGTTGCTTCGcca                                                                                                                                      >  1:854244/1‑145 (MQ=255)
            ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGAGGCCTGCCCGCCCCCACCCCGGCCGCTGCtttttaa                                                                                                                                       >  1:11671/1‑142 (MQ=255)
             gCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                <  2:653167/151‑1 (MQ=255)
               aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                              <  2:183062/151‑1 (MQ=255)
                 ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCTCCCGGCCGGTTCTTTGCAACGTTCaaa                                                                                                                             >  1:384300/1‑150 (MQ=25)
                          aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTCAACCTCCCCCCCACCCCCCCGGCCGATGCGTTCCAACGCTTCATTTCGCTcccc                                                                                                                   >  2:1108529/1‑150 (MQ=255)
                            gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                 <  1:53023/151‑1 (MQ=255)
                              aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                               <  1:526732/151‑1 (MQ=255)
                              aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCACCCCCCCGGCCGTTGCTTCGCCAAGTTCAAATCCCCCCCccgcgg                                                                                                               >  1:685721/1‑151 (MQ=21)
                              aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTTCTTCGCAACCTTCCAAACCCCCCCccgcgg                                                                                                               >  1:673979/1‑151 (MQ=255)
                                   aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGGTGCTTCGCCACGGTTAAATCCCCTCCCCGCGGGttt                                                                                                           >  1:4586/1‑150 (MQ=18)
                                   aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCCTTCTCCCCACCGACCGCCCGGGCGGTGCCTTACCAACGTTCAATACCCTCCCCGCCGAtt                                                                                                            >  1:540839/1‑149 (MQ=255)
                                          gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                    <  2:1021940/150‑1 (MQ=255)
                                          gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                    <  1:617827/150‑1 (MQ=255)
                                          gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                    <  1:675318/150‑1 (MQ=255)
                                          gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCGGCCGCTGCTTTGCAACGTTCAAATCCGCCCCCCGCGGGTTTTTcc                                                                                                       >  1:865124/1‑147 (MQ=18)
                                            aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg                                                                                                 <  2:498449/151‑1 (MQ=255)
                                            aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg                                                                                                 <  2:1016485/151‑1 (MQ=255)
                                             tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGACGTCCTCCCCGCCACCCTCCCGGCCGTTTTTTCGCCACGTTCAAAAACTCTACCCGCGGCGTTGTCCTACt                                                                                                   >  1:166028/1‑148 (MQ=255)
                                               atttgtagAAACGCAAAAAGGCCATCCGTCAGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCCGCCCGTGCTTTGCAAAGTTCAAAATCCCCCCCCGCGGGGTTGTTCTTCTCcgg                                                                                               >  2:820282/8‑150 (MQ=2)
                                                tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCACCCCCCCGGCCGTTGCTTTGCAACCTTCAAAACCCCCCCccgcggggtttttct                                                                                                      >  1:811841/1‑133 (MQ=255)
                                                  gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCCCCCCCCCCCCGGCCGTTCTTTCCCAACCTTCACATTCTCTCCCGGCGGCATTGTCCTCCCCgggcgg                                                                                            >  2:472164/1‑147 (MQ=17)
                                                    tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgt                                                                                          <  2:729421/150‑1 (MQ=255)
                                                     cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                         <  2:471260/150‑1 (MQ=255)
                                                     cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTGCTTTGCAAAGTTCAAATCCCCTCCCCGCGGGTTTTTCCTCTTCGGGGGAgtg                                                                                         >  2:1078295/1‑150 (MQ=17)
                                                       aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGGTGCTTCGCCACGTTCAAAACCCCCCCCCGCGGGTTTTTTCTACTCCGGGGAgtggt                                                                                       >  1:430445/1‑148 (MQ=17)
                                                        aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGAGTTCTGACCGCCAACCACCCCGCCGTTTTTCCTCAACATTCCAATCCGCTCCCGGGGGCtttcttctcatcagct                                                                                              >  1:54654/1‑129 (MQ=255)
                                                         aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCCCCCGCCCCCCCCCCGGCCGGTGCTTTGCAACAGTCAAAACCACCCCCCGCGGGGTTGTCATACTCTGGAGGgtg                                                                                         >  2:787655/1‑146 (MQ=255)
                                                             aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGCTTGTTTGCCAAGGTCAAAAACGCCCCCCGCGGGTTTGTTCTACTCGGGGGGGTGTTCTCCCac                                                                                >  2:602958/1‑151 (MQ=12)
                                                              aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGGCTGCCCCCCACCCCCCCGGCCGGTGCTTTGCAACCTTCAAAACCGCTCCCCGCGGGTTTGTTCTACTCCGGggaggggttcccgac                                                                                >  2:842059/1‑139 (MQ=255)
                                                              aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCTCCCGGGCGTTGCTTTGCCACGTTCCAATCCCCTCCCCGCGGGGTTGTCCTCCTCCGGAGAGTGTTTACCGac                                                                                >  1:90901/1‑150 (MQ=17)
                                                              aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCCACGTTCAAAACCGCCCCCGGCGGGTTTGTTCTTCTCCGGAGAGTGGTTAcc                                                                                   >  1:1009488/1‑147 (MQ=17)
                                                               aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                              <  1:1028694/151‑1 (MQ=255)
                                                               aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                              <  2:997344/151‑1 (MQ=255)
                                                               aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCACCCCCCCGGCCGGTGCTTTTCCAACGTCAAAACCCCCCCCCGCGGGTTTTTCCTCCTCCGGggaggggttccc                                                                                   >  2:556872/1‑138 (MQ=255)
                                                                aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                             <  1:952064/151‑1 (MQ=255)
                                                                aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                             <  1:753343/151‑1 (MQ=255)
                                                                 aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCCCCCGGCCCTTTCTTCGCCAACGTCAAATTCCCCCCCGCCGGGGTTTTCCTCCTCTCGGGGGGGGTTACCGacaaa                                                                             >  2:997956/1‑150 (MQ=255)
                                                                  ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                           <  2:18678/151‑1 (MQ=255)
                                                                       tCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCCCCCCCCCCCCGGCCGTTGCTTCGCAAAGTTCAAATCCCCCCCCCGCGGATTTGTCCTACTCGGGAGAGTGTTTACCCacaaacaa                                                                          >  1:518253/1‑147 (MQ=17)
                                                                          tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCTTTCCAACCGTCCAAATCCCCCCCCGCGGGGTTTTTCTACTCCGGGGGGGGGTTTCCCACAAACAACAAATa                                                                    >  2:631979/1‑150 (MQ=255)
                                                                             cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCGCCCCCCCCCGGGCCGGTGGTTCGCAAACTTTAAAACCGCCCCCGGCGGGCTTGTTCTTCTTCCGGGGGTGTTTGCCCCCAAACAAAAGaaaaaa                                                                 >  2:524414/1‑150 (MQ=9)
                                                                             cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCACCCCCCCGGCCGGTGCTTTGCAACGCTCAAATTCCCCCCCCGCGGGGTTTTCCTACTCCGGGGAGGGGTCACCCCCAAACAACAAATaaaa                                                                 >  1:882921/1‑150 (MQ=255)
                                                                                  ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAg                                                           <  2:143274/151‑1 (MQ=255)
                                                                                  ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCCGTGCTTCGCAAACTTCAAATCCCCACCCCCCGGGATTGTCCTCCTCGGGGGAGGGTGTACCCCCAAACAACAACaaaaaaaaaa                                                            >  2:884071/1‑150 (MQ=255)
                                                                                      tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCGGGCCGCTGCTTCGCAACCTTCAAATCCGCTCCCCGCGGGTTTTTCCTACTCGGGGGAGGGTTTTCCGACCAACAACAGATAaaaaaaaaaaccc                                                       >  2:767331/1‑146 (MQ=255)
                                                                                         cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                     <  1:833455/150‑1 (MQ=255)
                                                                                               atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                              <  1:1079175/151‑1 (MQ=38)
                                                                                                 ttgatGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCCCCCTCCGGGCCCTTTCTTCGCAACGTTCAAATCCTCCCCCGGCGGATTTTTCCTTCTCCGGGGAGGGTTTCCCCTCAAAAAACAAATAAAACAAAAAACCCa                                                      >  2:468826/1‑141 (MQ=2)
                                                                                                     tGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTTCTTTGCAACGTTCAAATCCCCTCCCCGCGGGTTTTTTCTTCTCGGGGGAGGGTTTACCCACAAAAAACAAATAAAAAAAAAAGGCCCGTCTTCCgc                                             >  1:235954/1‑145 (MQ=2)
                                                                                                      gTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                       <  2:960032/151‑1 (MQ=37)
                                                                                                          ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                    <  2:210035/150‑1 (MQ=35)
                                                                                                            cAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                                 <  2:49652/151‑1 (MQ=37)
                                                                                                                    ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                         <  2:1045986/151‑1 (MQ=38)
                                                                                                                       gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                      <  2:793518/151‑1 (MQ=38)
                                                                                                                           gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt                  >  1:220397/1‑151 (MQ=40)
                                                                                                                             ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                 <  1:755428/150‑1 (MQ=40)
                                                                                                                               tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc               >  1:224085/1‑150 (MQ=40)
                                                                                                                                gCCCGCCGCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc               <  2:700438/149‑1 (MQ=255)
                                                                                                                                gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc               >  2:312671/1‑149 (MQ=40)
                                                                                                                                    gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc         >  2:219901/1‑151 (MQ=40)
                                                                                                                                         cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAc     <  2:1069665/151‑1 (MQ=40)
                                                                                                                                            tCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct  >  2:292738/1‑151 (MQ=40)
                                                                                                                                             ccGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct  <  1:897599/150‑1 (MQ=40)
                                                                                                                                               |                                                                                                                                                      
AGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCT  >  NC_000913/3423155‑3423444

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