Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,575	0	A	C	44.3%	102.1 / 18.6	61	V555G (GTG→GGG)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (30/4);  new base C (0/27);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 2.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATT  >  NC_000913/501428‑501723                                                                                                                                                     |                                                                                                                                                     aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCt                                                                                                                                                   >  2:599921/1‑151 (MQ=255)   gTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTc                                                                                                                                                  >  1:350533/1‑150 (MQ=255)   gTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTCGCTCTCGATGCCTGATGCGACCCTTGCGCGTCTTCTCATGCCTCCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTc                                                                                                                                                  <  2:2876/150‑1 (MQ=255)      aTCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCAccc                                                                                                                                              <  1:109206/151‑1 (MQ=255)          tcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACccgcc                                                                                                                                           <  2:528543/150‑1 (MQ=255)             gTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGt                                                                                                                                        >  2:656861/1‑150 (MQ=255)                cGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCg                                                                                                                                    >  1:428454/1‑151 (MQ=255)                      ccAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTcc                                                                                                                               >  1:278469/1‑150 (MQ=255)                       cAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTc                                                                                                                                >  2:581305/1‑148 (MQ=255)                              ggggTTCACGCCTCATCTGTCATTTGGCTCTCGGTGCCTGATGCGCCGCGGGCGCTTCTTATCATCCCTACTTTTTTTTCATATTTTCCCTCCGGCAACCACCTTTTCCCCCGT‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                                      <  2:811185/148‑1 (MQ=255)                              ggCTTTCACGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCACCCACCTTTTCCCCCGC‑CACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                                      <  2:444829/151‑1 (MQ=37)                              ggCTTTCACGCCGCATCTGGCTTTTGGCTCCCGATGCCTGATGCGACGCGGGCGCGTCTTATCATGCCTAAATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGC‑CCCCCCCTCCCCCGCCGGTGGCGTTTCcagcag                                                                                                                          <  1:837803/147‑1 (MQ=255)                                    cACGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATCCCTACATTTTTTTAATTTTTTCCACCCGGCACCCACCGTTTCCCCCGT‑CCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATg                                                                                                                <  2:608277/151‑1 (MQ=255)                                    cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATg                                                                                                                >  2:612332/1‑151 (MQ=255)                                      cccctcATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTCTTACCATTCCTACATATTTTTCATATTTTACATCCGGCAACCGCCTTTTACCCCGT‑CCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                              <  1:122013/146‑1 (MQ=255)                                        ccgcATTTGGCTTTTGGCTCTGGTGGCTTGATGCGACGCGGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCACCCACCGTTTACCCCGC‑CACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                              <  1:257861/149‑1 (MQ=255)                                        ccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAc                                                                                                            >  2:164953/1‑151 (MQ=255)                                        ccgcATCTGGCATTTGGCTCTCAATGCCTGATGCGCCGCGGGCGCTTTTTATCATGCCTACATATTTTTCATATTTTACATCCGCCAACCACCTTTTCCCCCGT‑CACCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                                               <  1:503184/148‑1 (MQ=255)                                        cccatttgggCTTTTGCCTCCCTATCCCGAATGCGCCCCTGGCGCGTCTTCCCTTCCCTACCTTTTTTACATTTTTTCCACCCGGCAACCCCCTTTCCCCCCGC‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                                  <  2:390734/137‑1 (MQ=255)                                              cTGACATTGGGCTCTCGATCCTTGTGGCGACGCGGGGGCGTTTTATCAGCCCTACATTTTTTTCATATTTCACATCCGGCAACCACCTTTTCCCCCTT‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATg                                                                                                                <  1:802447/141‑1 (MQ=255)                                                   aTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                                  >  1:804569/1‑150 (MQ=255)                                                    tttGGCCCTCGATGCCTGATGCGCCGCTGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCACCTTTTCCCCCGC‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                                  <  2:759348/149‑1 (MQ=255)                                                       ggCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                              >  2:567266/1‑150 (MQ=255)                                                       ggCTCTCGATGCCTGATGCGACGCGGGCGCGTCTTATCATCCCTACATTTTTTTCATATTTTACATCCGGCACCCACCGTTTCCCCCGT‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                              <  1:749074/150‑1 (MQ=255)                                                               aTGCCTGATGCGCCGCTGGCGCTTCTTATCATCCCTACATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGC‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                                                      <  1:913639/150‑1 (MQ=255)                                                                tGCCTGATGCGCCGCTGGCGCTTCTTTTCATCCTTAAATTTTTTTCATTTTTACCACCCGGCACCCCCCGTTTCCCCCGT‑CCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                                                                                        <  2:228201/147‑1 (MQ=255)                                                                tGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                                                     >  1:286693/1‑150 (MQ=255)                                                                     gATGCGACGCTGGCGCTTTTTATCAGGCCTACATATTTTTCATTTTTTACACCCGCAAACCCCCTTTTCCCCCTT‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                                                     <  2:852111/145‑1 (MQ=255)                                                                        gCGACGCTGGCGCGTCTTATCATGCTTACATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGT‑CCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                               <  1:453943/148‑1 (MQ=255)                                                                           aCGCTGGCGGGTTTTTTCATGCCTCCATATTTTACATATTTAACACCCGGCACCCCCCTTTTCCCCCGT‑CCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                                <  2:544364/144‑1 (MQ=255)                                                                            cGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCg                                                                        >  2:929664/1‑151 (MQ=255)                                                                             gcggccgctTCTTATCATCCCTACATTTTTTTCATATTTTCCACCCGGCACCCCCCTTTCCCCCCTC‑CCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                         <  2:632877/144‑1 (MQ=255)                                                                             gCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGc                                                                       >  2:256363/1‑151 (MQ=255)                                                                                ggCGCGTTTTATCATGCCTACATTTTTTTCATTTTTTCCACCCGGCAACCCCCTTTTCCCCCGT‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGca                                                                      <  1:757528/149‑1 (MQ=255)                                                                                ggCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                                                    >  2:353823/1‑151 (MQ=255)                                                                                 gcgcTTTTTATCATCCCTACTTTTTTTTCATTTTTTCCCCCCGGC‑ACCCCCGTTTTCCCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCt                                                                              <  2:453492/140‑1 (MQ=255)                                                                                    cgTTTTATCATCCCTACATTTTTTTCTTTTTTTCCATCGGCCACCCCCCGTTTCCCCCCT‑CCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcaccac                                                                  <  2:52063/149‑1 (MQ=255)                                                                                      ttttATCCTGCCTACATTTTTTTCATATTTAACATCCGGCACCCCCCGTTTCCCCCGT‑CCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                               <  2:853312/148‑1 (MQ=255)                                                                                       cTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                             >  1:585388/1‑151 (MQ=255)                                                                                              tGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                      >  2:305636/1‑151 (MQ=255)                                                                                              tGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                      >  2:403891/1‑151 (MQ=255)                                                                                               gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                      >  1:599959/1‑150 (MQ=255)                                                                                                  tACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                                  >  1:194988/1‑151 (MQ=255)                                                                                                    cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                                 >  1:402562/1‑150 (MQ=255)                                                                                                      tttttttCATTTTTTCCACCCGGCACCCCCCGTTTCCCCCGC‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                      <  1:524919/141‑1 (MQ=255)                                                                                                      tatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                              >  2:542332/1‑151 (MQ=255)                                                                                                        tttttCATATTTTACATCCGGCAACCCCCGTTTCCCCCGC‑CCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                              <  1:524539/149‑1 (MQ=255)                                                                                                        tttttCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                                             >  1:815654/1‑150 (MQ=255)                                                                                                         ttttCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                                            >  2:168114/1‑150 (MQ=255)                                                                                                                 ttttACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                                    >  1:677074/1‑150 (MQ=255)                                                                                                                     aCATCCGGCACCCCCCGTTTCCCCCGT‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                                <  1:350554/150‑1 (MQ=255)                                                                                                                     aCATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                               >  2:557161/1‑151 (MQ=255)                                                                                                                       aTCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                             >  1:167901/1‑151 (MQ=255)                                                                                                                         ccGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc                           >  1:289463/1‑151 (MQ=255)                                                                                                                            caaccccccttttCCCCCGC‑CCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg                         <  1:656954/140‑1 (MQ=255)                                                                                                                                cccccGTTTCCCCCGC‑CCCCCCCCCCCCCGCCGGTGGTGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc                       <  2:632408/145‑1 (MQ=255)                                                                                                                                      tttCCCCCGC‑CCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct               <  2:705536/150‑1 (MQ=255)                                                                                                                                          ccccGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAat           >  1:596214/1‑150 (MQ=255)                                                                                                                                           cccGT‑CCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAatat         <  1:600011/151‑1 (MQ=255)                                                                                                                                            ccGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATc        >  2:583571/1‑151 (MQ=255)                                                                                                                                             cGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATcc       >  1:12663/1‑151 (MQ=255)                                                                                                                                                   ccaccTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAtt  >  2:302319/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGT‑CACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATT  >  NC_000913/501428‑501723

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: ‑