Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,564	0	A	C	44.1%	113.4 / 18.6	68	*559E (TAA→GAA)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (34/4);  new base C (0/30);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 5.19e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAAT  >  NC_000913/501419‑501714                                                                                                                                                   |                                                                                                                                                       gggATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCg                                                                                                                                                   >  2:80506/1‑151 (MQ=255)   gATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTc                                                                                                                                                 <  1:254546/151‑1 (MQ=255)     tCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTcac                                                                                                                                               >  2:1006682/1‑151 (MQ=255)      cGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTcacc                                                                                                                                              >  1:1226418/1‑151 (MQ=255)       gTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTcacc                                                                                                                                              >  1:503870/1‑150 (MQ=255)             ttATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCAc                                                                                                                                       >  2:1061739/1‑151 (MQ=255)             ttATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCAc                                                                                                                                       >  2:536374/1‑151 (MQ=255)               aTCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCAccc                                                                                                                                     <  2:60282/151‑1 (MQ=255)                 cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACccgc                                                                                                                                   <  1:662802/151‑1 (MQ=255)                 cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCATACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACccgc                                                                                                                                   >  1:220223/1‑151 (MQ=255)                     ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGt                                                                                                                               >  2:879785/1‑151 (MQ=255)                        tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGc                                                                                                                            >  2:667266/1‑151 (MQ=255)                          gTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCAGGC‑AACCACCGTTTACCCCGTCACCACCTAACCCGCCGGTGGCGt                                                                                                                          >  1:717893/1‑151 (MQ=255)                             ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                                       >  1:922847/1‑151 (MQ=255)                             ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                                       >  1:922857/1‑151 (MQ=255)                             gccaaGAAAAGGCTTTCCCGCCGCTTCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATTTTTTACATCCGGC‑AACCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGtt                                                                                                                         <  1:80500/145‑1 (MQ=25)                                 aGATAGGGCGTTCACCCCGCATCTGGCTTTTGGCTCTCGATGCATGATGCGACGCGGGCGCGTTTTATCTTCCCTACATATTTTTCATATTTTACTCCCGGC‑AACCCCCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCcag                                                                                                                    <  2:503915/150‑1 (MQ=255)                                 aGATAAGGCTTTCACGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGC‑ACCCACCTTTTCCCCCGTCCCCACCCCACCCGCCGGTGGCGTTTCcagc                                                                                                                   <  2:616790/151‑1 (MQ=38)                                   aaagggCGTCCCCGCCCCTTCTGGTTTTTGGCTCTCGTTGCCTGATGCGACGCGGGCGCTTCTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGC‑ACCCCCCGTTCCCCCCTCCACCCCCCCACCCGCCGGTGGCGTTTCcagca                                                                                                                  <  1:1214980/146‑1 (MQ=255)                                     agggcttcccccccgcATCGGGCTTTGGGCTCTCGTTGCCTGAGGGGCCGCGGGCGCTTTTTATCTTCCCTAAATTTTTTTAATTTTTTCCACCCGGC‑ACCCCCCTTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCcag                                                                                                                    <  1:1063560/134‑1 (MQ=255)                                          gTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGc                                                                                                          >  1:539138/1‑151 (MQ=255)                                             cccgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGC‑ACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                                                                        <  2:751185/148‑1 (MQ=255)                                              aCGCCGCTTCTGGTTTTGGCCTTTCTATGCCGGATGTGCCGCTGGCGCGTTTTATCAGGCCTACTTTTTTTTTATTTTTTACACCCGCC‑ACCCACCTTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                         <  1:595812/148‑1 (MQ=255)                                               cgccgcTTCTGGCTTTTGCCTCTCGATGCTTGAGGCGACGCGGGCGCGTCTTATCATCCCTCAATATTTTTCATATTTTACATCCGGC‑ACCCCCCTTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                     <  2:387264/151‑1 (MQ=255)                                               cgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                                      >  1:188072/1‑150 (MQ=255)                                               cgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                     >  2:67451/1‑151 (MQ=255)                                                  cgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACg                                                                                                  >  2:108548/1‑151 (MQ=255)                                                  cgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACg                                                                                                  >  1:235649/1‑151 (MQ=255)                                                        tGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                            >  2:609675/1‑151 (MQ=255)                                                             tttGTCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCTTCCCTACTTTTTTTTCATATTTTCCTCCCGCC‑ACCCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                                          <  2:916856/148‑1 (MQ=255)                                                             tttGGCTCTCGATGCCTGATGCGACGCTGGCGCTTCTTATCATCCCTACATATTTTTCATATTTTCCACCCGGC‑ACCCCCCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                                          <  2:223173/148‑1 (MQ=255)                                                             tttGCCTCTCGAGGCCTGATGCGCCGCTGGCGCGTTTTATCACGCCTACATATTTTTCATTTTTTACA‑CCGGCAAACCCCCGTTTCCCCCGGCACCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                         <  1:49903/149‑1 (MQ=255)                                                                     tcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                                                                               >  2:447716/1‑151 (MQ=255)                                                                      cGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                                                                               >  1:966787/1‑150 (MQ=255)                                                                           ccTGATGCGCCGCCGGCGCTTCTTATCCTGCCTATATTTTTTTCTTATTTTCCCCCCGAC‑ACCCACCTTTTCCCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                 <  1:214470/143‑1 (MQ=255)                                                                                 gCGACGCGGGCGCTTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGCC‑ACCCCCCTTTCCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                       <  1:705650/147‑1 (MQ=255)                                                                                  cgccgcgGGCGCGTTTTATCATGCCTACATATTTTTCATATTTTACACCCGGC‑ACCCCCCTTTTCCCCCCTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTGCCt                                                                     <  2:922952/145‑1 (MQ=255)                                                                                   gACGCTGGCGCGTTTTATCATGCCTACATTTTTTTCATTTTTTACATCCGGC‑ACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                                    <  1:692765/148‑1 (MQ=255)                                                                                    aCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                >  2:291695/1‑151 (MQ=255)                                                                                     cGCTGGCGCTTTTTTTCTTGCCTCCATTTTTTTCATATTTTCCACCCGCC‑ACCCCCCTTTTCCCCCCCCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCg                                                               <  1:891735/151‑1 (MQ=255)                                                                                     cGCTGGCGCATCTTATCGTCCCTACATATTTTTCATTTTTCACATCCGGC‑ACCCCCCGTTTACCCCTCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                <  1:67446/150‑1 (MQ=255)                                                                                      gCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTCCACCCGGC‑ACCCCCCGTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCt                                                                     <  2:922942/144‑1 (MQ=255)                                                                                         ggCGCTTCTTATCATCCCTACATTTTTTTCATTTTTTCCATCCGGC‑ACCCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTa                                                                        <  2:678881/138‑1 (MQ=255)                                                                                             cgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                                       >  2:312205/1‑151 (MQ=255)                                                                                               ttttATCATGCCTACATTTTTTTCATATTTTACCCCCGCC‑ACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                      <  1:247609/148‑1 (MQ=255)                                                                                                  tATCATGCCTACATTTTTTTCATTTTTTACACCCGGC‑ACCCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                                   <  2:993067/150‑1 (MQ=255)                                                                                                     cATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                               >  2:1154718/1‑151 (MQ=255)                                                                                                     cATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                               >  1:280614/1‑151 (MQ=255)                                                                                                       tGCCTACATATTTTTAATATTTTCCACCCGGC‑ACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                              <  1:834593/150‑1 (MQ=255)                                                                                                        gCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                             >  2:693835/1‑150 (MQ=255)                                                                                                              ttttttttCATTTTTAACACCCGGC‑AACCCCTTTTCCCCCCCTCCCCACCGCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                               <  2:1184226/139‑1 (MQ=255)                                                                                                              atttttttCATTTTTTCCACCCGGC‑ACCCACCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCa                                       <  1:1006582/147‑1 (MQ=255)                                                                                                              atatTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                      >  2:60052/1‑151 (MQ=255)                                                                                                               tatTTTTCATATTTTACACCCGGC‑ACCCCCCGTTTCCCCCGCCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                      <  2:863648/150‑1 (MQ=255)                                                                                                                 tttttCATTTTTTCCACCCGCC‑ACCCCCCTTTTCCCCCTCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                                    <  1:1061632/150‑1 (MQ=255)                                                                                                                  ttttCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                                   >  2:377711/1‑150 (MQ=255)                                                                                                                          ttttCCACCCGGC‑ACCCCCCTTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCa                            <  2:590195/149‑1 (MQ=255)                                                                                                                           tttACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                          >  1:510162/1‑150 (MQ=255)                                                                                                                            ttacccccGGC‑ACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATaa                         <  1:667206/144‑1 (MQ=255)                                                                                                                              aCATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                      >  1:391660/1‑151 (MQ=255)                                                                                                                               cacccGGC‑ACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                      <  1:693774/147‑1 (MQ=255)                                                                                                                               cATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATggg                     >  1:359426/1‑151 (MQ=255)                                                                                                                                ccccGCC‑ACCCCCCTTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                    <  1:377673/149‑1 (MQ=255)                                                                                                                                aTCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                    >  1:174648/1‑151 (MQ=255)                                                                                                                                 tCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg                   >  1:540153/1‑151 (MQ=255)                                                                                                                                    ggC‑ACCCACCCTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg                <  2:717957/151‑1 (MQ=255)                                                                                                                                        aaCCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc              >  1:1190271/1‑150 (MQ=255)                                                                                                                                            accGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata          >  2:577888/1‑150 (MQ=255)                                                                                                                                                   aCCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAat  >  1:670572/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAAT  >  NC_000913/501419‑501714

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