Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,678,043	0	C	T	45.3%	67.3 / 35.4	64	E107K (GAA→AAA)	yphE	putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (6/29);  new base T (29/0);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.34e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGT  >  NC_000913/2677903‑2678193                                                                                                                                             |                                                                                                                                                        cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCCGAGGGTTTCCCCCCCCGTc                                                                                                                                               >  1:131461/1‑151 (MQ=255)     gggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCTTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGcc                                                                                                                                           <  2:680277/151‑1 (MQ=255)      ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGAGGGTTTTCGCCACCGTCAGccc                                                                                                                                          >  2:153819/1‑151 (MQ=255)       gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGAAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCTCCCCGGGCAGAGGGTTTTCGCCCCTTTTAGccc                                                                                                                                          >  1:378474/1‑150 (MQ=255)         ttAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCCAGGCAGAGGGTTTTCGCCACTGTCCGcccct                                                                                                                                        >  1:374989/1‑148 (MQ=255)          tAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACCGGGGGAGAGGGTTTTCGCCCCTGTCAGACCTTcc                                                                                                                                      >  2:454931/1‑151 (MQ=255)             cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCTACCCAGGCCGGGGGTTTTCCCCCCCCTCAGCCCCTccccc                                                                                                                                   >  1:880833/1‑151 (MQ=255)                  aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑t                                                                                                                             <  2:464186/151‑1 (MQ=255)                      aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTTTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTaa                                                                                                                          <  1:994426/150‑1 (MQ=255)                      aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCAGGCAGAGGGTTTTCGCCCCCGTCAGCCCCTCCACCCGAC‑CTaa                                                                                                                          >  1:777949/1‑150 (MQ=255)                         ttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACCGGGGCGGAGGTTTTTCGCCACTCTCAGCCCCCccaccatccctaaatcc                                                                                                                       >  1:674661/1‑140 (MQ=255)                         ttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCAGGGAGAGGTTTTCCCCCCCCGTCAGCCCCTTCACCCGAC‑TTAACTcc                                                                                                                      >  2:754183/1‑151 (MQ=255)                             tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGAGGGTTTTCGCCCCCGTCAGCCCCTCCACCCAAC‑TTAACTCCTGAt                                                                                                                  >  2:1055535/1‑151 (MQ=255)                                     aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACcgcg                                                                                                          <  1:516097/151‑1 (MQ=255)                                     aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACcgc                                                                                                           <  2:570371/150‑1 (MQ=255)                                     aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGAC‑TTAACTCCTGATAAAAcgc                                                                                                           >  2:1073979/1‑150 (MQ=255)                                      aCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCCCCCCGGCCGCGGGTTTTCTCCCCCGTCCGCCCTTCCCCCAGAC‑TTAAATCCCGTTaaaac                                                                                                             >  1:296736/1‑145 (MQ=255)                                       cGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAGTCAATCATGCCGTTGCGGCGGGGTCACCGACCGAGGCGAAGGGTTTTCGCCACTGTCAGCCCCTCCCCCAGAA‑TTAACCCCTGATAAAAcgc                                                                                                           >  1:928931/1‑148 (MQ=255)                                       cGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACcgc                                                                                                           >  2:433474/1‑148 (MQ=255)                                           aaCGCCCAACGCCTGGAAGCAACGTTGGGAATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGGGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCa                                                                                                      <  1:613762/149‑1 (MQ=255)                                            aCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGGAGAGGGTTTTCGCCACTGTCCGCCCCTCCCCCCGAC‑TTAACTCCTGGTAAACCCCGCGCcccc                                                                                                   >  1:247593/1‑151 (MQ=255)                                               cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCa                                                                                                <  2:5395/151‑1 (MQ=255)                                                ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTTACCCAGGCAGAGGTTTTTCGCCCCTGTCCGCCCTTCCCCCCGAC‑CTAACTCCTGATAAAACGCGCGCcccccc                                                                                                 >  2:541945/1‑149 (MQ=255)                                                ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAg                                                                                               <  1:838517/151‑1 (MQ=255)                                                 cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGCCCCAGGCCGAGGGTTTTCCCCCCCCTTAGCCCCTCCCCCCGTCTTTAACCCCT‑ATAAACCCCGCGCAcccccc                                                                                                >  1:919709/1‑148 (MQ=255)                                                 cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACGGAGGGAGGGGGTTTTCGCCACCTTCAGCCCCTCCCCCAGAC‑TTAACTCCCGATAAAACGCCCGCAccccc                                                                                                 >  1:885405/1‑148 (MQ=255)                                                 cAGCGCCTGTAAGCAACGTTGGGCATCCCGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCAGGCAGAGGTTTTTCCCCACCGTCAGCCCCTCCCCCAGAC‑CTAACTCCTGATAAACCCCGCCCAccccc                                                                                                 >  2:486448/1‑148 (MQ=255)                                                      ccTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAg                                                                                          <  1:1025593/150‑1 (MQ=255)                                                         gTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCg                                                                                        <  2:731902/149‑1 (MQ=255)                                                             gCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTCACCGAGGCAGAGGGTTTTCCCCCCTGTCAGCCCCTCCACCCGAC‑TTAACTCCTGATAAACCCCGCGCACCCCCCGCTCAGccg                                                                                       >  1:516255/1‑144 (MQ=255)                                                                aCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCCACCGAGGCAGAGGGTTTTCTCCCCTGTCCGCCCCTCCCCCCCAC‑TTAACTCCTGATAAAACCCGCGCACCCCCCGTTCAGCCGCCCGGc                                                                                 >  1:921765/1‑149 (MQ=255)                                                                 cGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCCCCCGAC‑TTAACTCCTGATAAACCCCGCGCCCCCCCCGTTCAGCGGCACGGc                                                                                 >  2:259248/1‑148 (MQ=255)                                                                  gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCCCCAGAC‑TTAACTCCTGATAAAACGCGCGCACCCCCCGCTCAGCGGCCCGGCGAg                                                                              >  2:661654/1‑150 (MQ=255)                                                                        cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCAGGCAGAGGGTTTTCGCCCCCGTCAGCCCTTCCACCAGAC‑TTAAATCCTGATAAAACGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCg                                                                        >  2:891191/1‑150 (MQ=255)                                                                        cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCCGGCAGAGGTTTTTCCCCCCTGTCCGCCCCTCCCCCAGAC‑TTAACTCCTGATAAAACGCCCCCACCCCCCCTTCCGCGGCCCGGCGAGTCAGCg                                                                        >  1:941546/1‑150 (MQ=255)                                                                         aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTa                                                                      <  1:249722/151‑1 (MQ=255)                                                                          tCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTCCCCCCCTGTCAGCCCTTCCACCCGAC‑TTAACTCCTGATAAAACCCGCGCCCCCCCAGTTCAGCGGCACCGCGAGGCAGc                                                                         >  2:386238/1‑147 (MQ=255)                                                                               cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCACGGCGCAGGGTTTTCCCCACCCCTTCCCCTTTCACCACAC‑TCCAAACCTGATAAACCGCGCGCACCCCCCGTTCACCTGCCCGCCGGGTCCGCGTa                                                                      >  1:334104/1‑145 (MQ=255)                                                                                 ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGt                                                               <  1:160424/150‑1 (MQ=255)                                                                                 ccATTTGCAGGTAATCAACCCTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc                                                              <  1:24578/151‑1 (MQ=255)                                                                                     ttGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGCGGCAGAGGGTTTTCCCCCCTGTCAGCCCCTCCCCCCGAC‑TTAACTCCCGATAAAACCCGCGCCCCCCCCGGTCAGCGGCACGGCGGGGCAGCGGAGCTTCGTCAcc                                                           >  2:1031889/1‑150 (MQ=255)                                                                                        cAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTcc                                                       <  2:880878/151‑1 (MQ=255)                                                                                          ggTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGAC‑TTAACTCCTGATAAAACCCGCGCACCCCCCCTTCAGCGGCACCGCGGGGCAGCGTAGCCTCGTCACCCTcc                                                       >  2:1053085/1‑149 (MQ=255)                                                                                             aatcaatcATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTc                                                   <  1:1031823/150‑1 (MQ=255)                                                                                                caatcaTGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTTCGCCACCGTCAGCCCTTCCCCCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCCTCCCGTc                                                   >  2:182962/1‑147 (MQ=255)                                                                                                caatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGt                                               <  1:844807/151‑1 (MQ=255)                                                                                                  atcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtct                                             <  2:828022/151‑1 (MQ=255)                                                                                                    caTGCCGTTGCGGCGGGGCCCCCGACCCAGGCAGAGGGTTTTCGCCCCCCTCAGCCCCTCCCCCAGAC‑TTAACTCCTGATAAAACCCGCGCACCCCCCCCTCAGCGGGACCGCGGGGCCGCGGAGCTTCGTCACCTTCCCCTCGCGcct                                             >  1:758007/1‑149 (MQ=255)                                                                                                     aTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCACCCGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGGGGCAGCGGAGCTTCGTCACCCTCCCGTCGCGTCTCGc                                          >  2:950075/1‑151 (MQ=255)                                                                                                         cGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAAt                                      <  2:904492/151‑1 (MQ=255)                                                                                                          gTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATc                                     <  2:716301/151‑1 (MQ=255)                                                                                                             gcggcgGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCa                                   <  2:663117/150‑1 (MQ=255)                                                                                                               ggcggGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGat                                <  1:290276/151‑1 (MQ=255)                                                                                                                         cTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                       <  1:957963/150‑1 (MQ=255)                                                                                                                           gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTa                      <  1:1079489/149‑1 (MQ=255)                                                                                                                           gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc                    <  2:2036/151‑1 (MQ=255)                                                                                                                              cGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCgg                 >  1:246388/1‑151 (MQ=255)                                                                                                                                   cAGAGTTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGt             <  1:636743/150‑1 (MQ=255)                                                                                                                                   cAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGtt            <  1:693336/151‑1 (MQ=255)                                                                                                                                    agagGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc           <  2:375009/151‑1 (MQ=255)                                                                                                                                     gagGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCg          <  2:988800/151‑1 (MQ=255)                                                                                                                                       ggTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt        <  1:757400/151‑1 (MQ=255)                                                                                                                                         ttttCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGc      >  1:760631/1‑151 (MQ=255)                                                                                                                                          tttCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGcc     >  2:957093/1‑151 (MQ=255)                                                                                                                                             ccGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGt  <  1:1343/151‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                        CGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGT  >  NC_000913/2677903‑2678193

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