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

AGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCT  >  NC_000913/3423155‑3423444
                                                                                                                                               |                                                                                                                                                      
agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCC‑CCCCTCCCGGcc                                                                                                                                                    >  2:663232/1‑147 (MQ=255)
    cGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTc                                                                                                                                           <  1:116897/151‑1 (MQ=255)
         gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTCGTCCTCGCCGCC‑ACCCTCCCGCCCGGTGA‑TTccc                                                                                                                                         >  1:791964/1‑146 (MQ=255)
         gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCC‑CCCCCCCCGGCCGGTGC‑TTCT‑Caa                                                                                                                                      >  2:97787/1‑150 (MQ=255)
            ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCC‑CCCCCCCCGGCCGTTGC‑TTCG‑cca                                                                                                                                      >  1:854244/1‑145 (MQ=255)
            ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGAGGCCTGCCCGCC‑CCCACCCCGGCCGCTGC‑tttttaa                                                                                                                                       >  1:11671/1‑142 (MQ=255)
             gCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTc                                                                                                                                <  2:653167/151‑1 (MQ=255)
               aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCaa                                                                                                                              <  2:183062/151‑1 (MQ=255)
                 ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCC‑CCCCTCCCGGCCGGTTC‑TTTG‑CAA‑CGTTCaaa                                                                                                                             >  1:384300/1‑150 (MQ=25)
                          aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGTCAACCTCCCCCCC‑ACCCCCCCGGCCGATGC‑GTTC‑CAA‑CGCTTCATTTCGCTcccc                                                                                                                   >  2:1108529/1‑150 (MQ=255)
                            gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCcggc                                                                                                                 <  1:53023/151‑1 (MQ=255)
                              aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCcggcgg                                                                                                               <  1:526732/151‑1 (MQ=255)
                              aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCC‑ACCCCCCCGGCCGTTGC‑TTCG‑CCA‑AGTTCAAATCCCCCCCccgcgg                                                                                                               >  1:685721/1‑151 (MQ=21)
                              aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCC‑CCCCCCCCGGCCGTTTC‑TTCG‑CAA‑CCTTCCAAACCCCCCCccgcgg                                                                                                               >  1:673979/1‑151 (MQ=255)
                                   aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCC‑CCCCCCCCGGCCGGTGC‑TTCG‑CCA‑CGGTTAAATCCCCTCCCCGCGGGttt                                                                                                           >  1:4586/1‑150 (MQ=18)
                                   aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCCTTCTCCCCACCGACCGCCCGGGCGGTGCC‑TTAC‑CAA‑CGTTC‑AATACCCTCCCCGCCGAtt                                                                                                            >  1:540839/1‑149 (MQ=255)
                                          gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                    <  2:1021940/150‑1 (MQ=255)
                                          gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                    <  1:617827/150‑1 (MQ=255)
                                          gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                    <  1:675318/150‑1 (MQ=255)
                                          gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCCCCCGGCCGCTGC‑TTTG‑CAA‑CGTTCAAATCCGCCCCCCGCGGGTTTTTcc                                                                                                       >  1:865124/1‑147 (MQ=18)
                                            aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg                                                                                                 <  2:498449/151‑1 (MQ=255)
                                            aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg                                                                                                 <  2:1016485/151‑1 (MQ=255)
                                             tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGACGTCCTCCCCGCC‑ACCCTCCCGGCCGTTTT‑TTCG‑CCA‑CGTTCAAAAACTCTACCCGCGGCGTTGTCCTACt                                                                                                   >  1:166028/1‑148 (MQ=255)
                                               atttgtagAAACGCAAAAAGGCCATCCGTCAGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCC‑CCCCCCCCCGCCCGTGC‑TTTG‑CAA‑AGTTCAAAATCCCCCCCCGCGGGGTTGTTCTTCTCcgg                                                                                               >  2:820282/8‑150 (MQ=2)
                                                tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCC‑ACCCCCCCGGCCGTTGC‑TTTG‑CAA‑CCTTCAAAACCCCCCCccgcggggtttttct                                                                                                      >  1:811841/1‑133 (MQ=255)
                                                  gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCCCC‑CCCCCCCCGGCCGTTCT‑TTCC‑CAA‑CCTTCACATTCTCTCCCGGCGGCATTGTCCTCCCCgggcgg                                                                                            >  2:472164/1‑147 (MQ=17)
                                                    tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgt                                                                                          <  2:729421/150‑1 (MQ=255)
                                                     cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                         <  2:471260/150‑1 (MQ=255)
                                                     cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCC‑ACCCCCCCGGCCGTTGC‑TTTG‑CAA‑AGTTCAAATCCCCTCCCCGCGGGTTTTTCCTCTTCGGGGGAgtg                                                                                         >  2:1078295/1‑150 (MQ=17)
                                                       aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCC‑ACCCCCCCGGCCGGTGC‑TTCG‑CCA‑CGTTCAAAACCCCCCCCCGCGGGTTTTTTCTACTCCGGGGAgtggt                                                                                       >  1:430445/1‑148 (MQ=17)
                                                        aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGAGTTCTGACCGCC‑AACCACCCCGCCGTTTT‑TCCT‑CAA‑CATTCCAATCCGCTCCCGGGGGCtttcttctcatcagct                                                                                              >  1:54654/1‑129 (MQ=255)
                                                         aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCCCCCGCC‑CCCCCCCCGGCCGGTGC‑TTTG‑CAA‑CAGTCAAAACCACCCCCCGCGGGGTTGTCATACTCTGGAGGgtg                                                                                         >  2:787655/1‑146 (MQ=255)
                                                             aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCC‑CCCCCCCCGGCCGCTTG‑TTTG‑CCA‑AGGTCAAAAACGCCCCCCGCGGGTTTGTTCTACTCGGGGGGGTGTTCTCCCac                                                                                >  2:602958/1‑151 (MQ=12)
                                                              aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGGCTGCCCCCC‑ACCCCCCCGGCCGGTGC‑TTTG‑CAA‑CCTTCAAAACCGCTCCCCGCGGGTTTGTTCTACTCCGGggaggggttcccgac                                                                                >  2:842059/1‑139 (MQ=255)
                                                              aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCC‑CCCCTCCCGGGCGTTGC‑TTTG‑CCA‑CGTTCCAATCCCCTCCCCGCGGGGTTGTCCTCCTCCGGAGAGTGTTTACCGac                                                                                >  1:90901/1‑150 (MQ=17)
                                                              aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCC‑CCCCCCCCGGCCGTTGC‑TTCG‑CCA‑CGTTCAAAACCGCCCCCGGCGGGTTTGTTCTTCTCCGGAGAGTGGTTAcc                                                                                   >  1:1009488/1‑147 (MQ=17)
                                                               aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                              <  1:1028694/151‑1 (MQ=255)
                                                               aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                              <  2:997344/151‑1 (MQ=255)
                                                               aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCC‑ACCCCCCCGGCCGGTGCTTTTC‑CAA‑CG‑TCAAAACCCCCCCCCGCGGGTTTTTCCTCCTCCGGggaggggttccc                                                                                   >  2:556872/1‑138 (MQ=255)
                                                                aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                             <  1:952064/151‑1 (MQ=255)
                                                                aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                             <  1:753343/151‑1 (MQ=255)
                                                                 aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCC‑ACCCCCCCGGCCCTTTC‑TTCGCCAA‑CG‑TCAAATTCCCCCCCGCCGGGGTTTTCCTCCTCTCGGGGGGGGTTACCGacaaa                                                                             >  2:997956/1‑150 (MQ=255)
                                                                  ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                           <  2:18678/151‑1 (MQ=255)
                                                                       tCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCCCC‑CCCCCCCCGGCCGTTGC‑TTCG‑CAA‑AGTTCAAATCCCCCCCCCGCGGATTTGTCCTACTCGGGAGAGTGTTTACCCacaaacaa                                                                          >  1:518253/1‑147 (MQ=17)
                                                                          tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCC‑CCCCCCCCGGCCGTTGC‑TTTC‑CAACCGTCCAAATCC‑CCCCCCGCGGGGTTTTTCTACTCCGGGGGGGGGTTTCCCACAAACAACAAATa                                                                    >  2:631979/1‑150 (MQ=255)
                                                                             cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCGCC‑CCCCCCCGGGCCGGTGG‑TTCG‑CAA‑ACTTTAAAACCGCCCCCGGCGGGCTTGTTCTTCTTCCGGGGGTGTTTGCCCCCAAACAAAAGaaaaaa                                                                 >  2:524414/1‑150 (MQ=9)
                                                                             cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCC‑ACCCCCCCGGCCGGTGC‑TTTG‑CAA‑CGCTCAAATTCCCCCCCCGCGGGGTTTTCCTACTCCGGGGAGGGGTCACCCCCAAACAACAAATaaaa                                                                 >  1:882921/1‑150 (MQ=255)
                                                                                  ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAg                                                           <  2:143274/151‑1 (MQ=255)
                                                                                  ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCC‑CCCCCCCCGGCCCGTGC‑TTCG‑CAA‑ACTTCAAATCCCCACCCCCCGGGATTGTCCTCCTCGGGGGAGGGTGTACCCCCAAACAACAACaaaaaaaaaa                                                            >  2:884071/1‑150 (MQ=255)
                                                                                      tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCC‑CCCCCCCGGGCCGCTGC‑TTCG‑CAA‑CCTTCAAATCCGCTCCCCGCGGGTTTTTCCTACTCGGGGGAGGGTTTTCCGACCAACAACAGATAaaaaaaaaaaccc                                                       >  2:767331/1‑146 (MQ=255)
                                                                                         cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                     <  1:833455/150‑1 (MQ=255)
                                                                                               atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                              <  1:1079175/151‑1 (MQ=38)
                                                                                                 ttgatGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCC‑CCCCTCCGGGCCCTTTC‑TTCG‑CAA‑CGTTCAAATCCTCCCCCGGCGGATTTTTCCTTCTCCGGGGAGGGTTTCCCCTCAAAAAACAAATAAAACAAAAAACCCa                                                      >  2:468826/1‑141 (MQ=2)
                                                                                                     tGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCC‑CCCCCCCCGGCCGTTTC‑TTTG‑CAA‑CGTTCAAATCCCCTCCCCGCGGGTTTTTTCTTCTCGGGGGAGGGTTTACCCACAAAAAACAAATAAAAAAAAAAGGCCCGTCTTCCgc                                             >  1:235954/1‑145 (MQ=2)
                                                                                                      gTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                       <  2:960032/151‑1 (MQ=37)
                                                                                                          ggCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                    <  2:210035/150‑1 (MQ=35)
                                                                                                            cAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                                 <  2:49652/151‑1 (MQ=37)
                                                                                                                    ggcgggcgTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                         <  2:1045986/151‑1 (MQ=38)
                                                                                                                       gggcgTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                      <  2:793518/151‑1 (MQ=38)
                                                                                                                           gTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt                  >  1:220397/1‑151 (MQ=40)
                                                                                                                             ccTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                 <  1:755428/150‑1 (MQ=40)
                                                                                                                               tGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc               >  1:224085/1‑150 (MQ=40)
                                                                                                                                gCCCGCC‑GCCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc               <  2:700438/149‑1 (MQ=255)
                                                                                                                                gCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc               >  2:312671/1‑149 (MQ=40)
                                                                                                                                    gCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc         >  2:219901/1‑151 (MQ=40)
                                                                                                                                         cccTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAc     <  2:1069665/151‑1 (MQ=40)
                                                                                                                                            tCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct  >  2:292738/1‑151 (MQ=40)
                                                                                                                                             ccGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct  <  1:897599/150‑1 (MQ=40)
                                                                                                                                               |                                                                                                                                                      
AGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCC‑ACCCTCCGGGCCGTTGC‑TTCG‑CAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCT  >  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: ‑