Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,443	0	A	C	42.6%	112.9 / 12.2	61	Y40S (TAC→TCC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (33/2);  new base C (0/26);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 3.10e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421301‑2421592                                                                                                                                                    |                                                                                                                                                      aaGTAAAACAACGGAGACCCTGCAATGGCAAAACCTGCCGCTATTATTCGTATAAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCccgccg                                                                                                                                               <  1:670806/151‑1 (MQ=255)                 aaCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                                 >  1:249812/1‑149 (MQ=255)                   ccgaaaaggccacaCCCTCCCGTTATTTTTCGTAAAAAAACCC‑TT‑CGTTTGCTTACTTTTATCGAAATAAA‑GAAAAAAAAAATAAC‑CACCCGTCAGA‑ATTTTGTTATAAAGGTGACGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             <  1:76907/141‑1 (MQ=255)                       cAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAT‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                            >  2:909788/1‑148 (MQ=255)                        actgccccacccggCCGTTTTTTTTCGTATAAAGACCC‑TT‑CTTTTGCGTCCTTTTATCGAAATTAA‑GAAAAAAGAAATAAA‑CAACCGTAAGT‑TTTTTTTTAAAAATGAGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                        <  2:515029/138‑1 (MQ=255)                          tGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                                      >  2:2378/1‑151 (MQ=255)                           gccacacccgcccGTTATTTTTCG‑AAAAAAAACC‑TTCCTTTTTCGTCTGTTTATCGAAATAAG‑GAAAAAAAAAATAAA‑CACCCGTAAGG‑ATTTTTTTATAAAGGGGACGTCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                          <  1:598113/135‑1 (MQ=255)                                acccTGCCTTTATTATTCTTAAAAAAACCC‑TT‑CTTTTGCGTTCTTTTACCGAAATAAA‑GAAAAAAGAAATAAA‑CCCCCG‑CAGGAATTTTTTTATAATGGTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATc                                                                                                                     <  2:5259/144‑1 (MQ=255)                                acccTGCCGTTTTTTTTCGTAAAAAAATCC‑TT‑CTTTTGCTTACTTTTATGGAAATAAA‑GAAAGAAAAAATAAA‑CACCCGTCAGA‑ATTTTTTTATAATTGTGACCATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGa                                                                                                                <  2:361505/149‑1 (MQ=255)                                 cccTGCCTTTATTTTTC‑TAAAAAAAACCTTT‑CTTTGGCTTACGTTTA‑CGGAAATAA‑AGAAAAAAAAAATAACCACCCGTCAGA‑ATATTTTAATAAATATGACAACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                                  <  1:220882/147‑1 (MQ=255)                                 cccTGCCGCTATTTTTGTTAAAAAGACCT‑TT‑CGTTTGGGTACTTTTATCGAAATAAA‑GAAAAAAAAAATAAA‑CACCCGTAAGA‑ATATGTTAATAAAGGTGACATCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGa                                                                                                                <  1:654432/149‑1 (MQ=255)                                    tGCCGCTATTATTCGTATAAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                            >  1:894559/1‑151 (MQ=255)                                     ccatttattatTTGTAAAAAGACCT‑TT‑CTTTTGCTGAGTTTTATCGAAAAAAA‑GAAAGAAAAAAAAAA‑CACCCGGAAGA‑ATATTTTTATAAAGGCGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcgc                                                                                                           <  1:771536/146‑1 (MQ=255)                                      cctttattttTCTTAAAAAGACCT‑TT‑TGTTTGCGTAGTTTTTGCGAAATAAA‑GAAAAAAAAAATAAA‑CACCCTTAAGT‑ATTTTTTAATAATTGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                          <  2:609894/147‑1 (MQ=255)                                               ttCGTATAAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                                 >  2:559871/1‑151 (MQ=255)                                               ttCGTAAAAAGAACC‑TT‑CTTTTGCGTCCTTTTATCGAAATTAA‑GGAAAAAAAAATTAA‑CACCCGGAAGA‑ATTTTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                                  <  2:451762/150‑1 (MQ=255)                                                 cGTATAAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                                >  2:185966/1‑150 (MQ=255)                                                    ataAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            >  1:141866/1‑151 (MQ=255)                                                    ataAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            >  2:33072/1‑151 (MQ=255)                                                     taAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           >  2:276328/1‑151 (MQ=255)                                                        aGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                        >  2:659955/1‑151 (MQ=255)                                                                tttttttGCTTCCTTTTATAGAAATTAA‑GGAAAAAAAAATTAC‑CAACCGTCAGG‑ATTTTGTTAACAATGAGAAGATCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                                   <  2:174920/146‑1 (MQ=255)                                                                          ttgttttATCGAAATAAA‑GAAAGAAAAAATAAC‑CACCCGGAAGG‑ATTTTTTTATAAATGTGCCGACCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt                                                                           <  1:393353/144‑1 (MQ=255)                                                                            cGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                       >  2:148135/1‑150 (MQ=255)                                                                             gTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                     >  1:187319/1‑151 (MQ=255)                                                                              tttacggaaaTAAA‑GGAAAAAAAAATAAA‑CACCCGTAAGA‑ATATTTTAATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCa                                                                      <  2:584827/141‑1 (MQ=255)                                                                              tttATGGAAATAAA‑GAAAAAAAAAATAAA‑CACCCTTAAGT‑TTTTTTTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCa                                                                      <  2:591997/149‑1 (MQ=255)                                                                                tATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                   >  2:235734/1‑150 (MQ=255)                                                                                tATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                   >  2:830693/1‑150 (MQ=255)                                                                                 aTCGGAATAAA‑GAAAAAAAAAATAAA‑CACCCGTCAGG‑ATATTGTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                   <  2:862691/149‑1 (MQ=255)                                                                                  tCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                   >  2:519411/1‑148 (MQ=255)                                                                                   cGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 >  1:850438/1‑149 (MQ=255)                                                                                     gAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                              >  1:652403/1‑150 (MQ=255)                                                                                      aaTTAAGGGAAGAAGAAATAAA‑CACCCCCCAGG‑ATTTTGTTATAAAGGACGCCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                             <  2:505662/151‑1 (MQ=255)                                                                                         tAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                         >  1:853200/1‑151 (MQ=255)                                                                                         tAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                         >  2:576604/1‑151 (MQ=255)                                                                                          aa‑gaaagaAAAAATAAA‑CACCCGAAAGA‑ATTTTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATaa                                                          <  1:700612/149‑1 (MQ=255)                                                                                              gaagaagaaATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGtt                                                     >  1:260232/1‑151 (MQ=255)                                                                                              gaagaagaaATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGtt                                                     >  1:88639/1‑151 (MQ=255)                                                                                               aaaaaaaaaTAAA‑CACCCGTAAGA‑ATATTTTTATAAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                    <  1:360955/145‑1 (MQ=255)                                                                                                  aagaaATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtctc                                                 >  1:312728/1‑151 (MQ=255)                                                                                                  aaaaaaTAAC‑CACCCGTCAGA‑ATTTTGTTATCAATGTAACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                                   <  2:723920/146‑1 (MQ=255)                                                                                                    aaaaTAAA‑CACCCGTAAGT‑ATTTTGTTATCAATGTGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               <  2:438804/150‑1 (MQ=255)                                                                                                     aaaTAAA‑CACCCGCCAGG‑ATTTTGTTATCAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               <  2:821297/150‑1 (MQ=255)                                                                                                             caacCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                         >  2:734014/1‑149 (MQ=255)                                                                                                               acCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGAAGAGAATAACCGTTTCTCTTTGCTGGaaa                                      >  2:529736/1‑150 (MQ=255)                                                                                                                ccGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                                     >  1:741163/1‑150 (MQ=255)                                                                                                                ccGTAAGG‑ATTTTGTTATAAATGGGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaaa                                    <  2:702365/151‑1 (MQ=255)                                                                                                                 cGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAt                                   >  1:898845/1‑151 (MQ=255)                                                                                                                 cGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAt                                   >  2:592172/1‑151 (MQ=255)                                                                                                                          tatTTTTATAAATGTGACCATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCa                            <  1:925867/150‑1 (MQ=255)                                                                                                                          tatTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             >  2:823554/1‑149 (MQ=255)                                                                                                                          tatTGTTATAAATGTGCCGTCCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACt                              <  2:569894/148‑1 (MQ=255)                                                                                                                            tttttATAAATGGAACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAg                           <  1:372697/146‑1 (MQ=255)                                                                                                                                 ataaaTGTGACGCCCCACACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgt                      <  1:576349/146‑1 (MQ=255)                                                                                                                                  tCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATCATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCt                   >  2:536868/1‑151 (MQ=255)                                                                                                                                    aaTGGGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCt                   <  2:948789/149‑1 (MQ=255)                                                                                                                                        tgACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat             >  2:469398/1‑151 (MQ=255)                                                                                                                                               ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt      >  2:70506/1‑151 (MQ=255)                                                                                                                                                 aCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  1:251477/1‑150 (MQ=255)                                                                                                                                                   tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  2:117263/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                      AAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TT‑CGTTTGCGTACGTTTATCGGAATTAA‑GGAAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421301‑2421592

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑