Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,678,039	0	T	G	45.3%	104.6 / 16.3	64	N108T (AAC→ACC)	yphE	putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (6/29);  new base G (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

GCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTG  >  NC_000913/2677893‑2678188                                                                                                                                                   |                                                                                                                                                      gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGgttt                                                                                                                                                    >  2:916765/1‑150 (MQ=255) gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCTACCCAGGCCGAGGgttt                                                                                                                                                    >  1:1202052/1‑150 (MQ=255) gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGAGGgttt                                                                                                                                                    >  2:283159/1‑150 (MQ=255)   tgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGCAGGGGgttttc                                                                                                                                                  >  1:758678/1‑148 (MQ=255)    gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGCCCGCGGCAGAGGgtttt                                                                                                                                                   >  1:445190/1‑147 (MQ=255)         tGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGCGGCAGAGGGTTTCCGCCAc                                                                                                                                             >  1:348105/1‑149 (MQ=255)           cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGGAGAGGGTTTTCGcccccgt                                                                                                                                          >  2:741867/1‑145 (MQ=255)               gggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTTACCGAGGCAGAGGGTTTCCGCCACCGTCAGcc                                                                                                                                     >  2:951950/1‑151 (MQ=255)               gggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGcc                                                                                                                                     <  1:679453/151‑1 (MQ=255)                 gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGCAGAGGGTTTTCGCCCCCGTCCGccc                                                                                                                                    >  1:428368/1‑150 (MQ=255)                    tAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCAGGCAGAGGGTTTCCGCCACTGTTAAcccct                                                                                                                                  >  2:953233/1‑147 (MQ=255)                     aGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTcca                                                                                                                               <  1:238036/151‑1 (MQ=255)                        gTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTccacca                                                                                                                            <  1:1093759/151‑1 (MQ=255)                           gAAACAAGTTGTTCAGGACTAACGTCAACGCCCGGCGCCTGTAAGCAACGTTGGGCTCCCGGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGa                                                                                                                          <  1:819007/150‑1 (MQ=255)                           gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGa                                                                                                                          <  2:58701/150‑1 (MQ=255)                           gAAACAAGTTGCTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCAGGGAGAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCAGa                                                                                                                          >  1:186693/1‑150 (MQ=255)                            aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCGACTGACCGAGGCAAAGGGTTTTCGGCACCGGCAACCCTTctccc                                                                                                                             >  2:657461/1‑142 (MQ=255)                                  gttgtgCAGGTCTATCGTCAACGCCCAGCGCCTGTAAGCACCGTTGGGCATCCTGCGCCCTTTGCGGGTAATCAACCATGCCGTTGCGGCGGGGCCACTGACCGTGGCCGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTc                                                                                                                  <  1:261011/151‑1 (MQ=255)                                   ttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTc                                                                                                                  <  1:890730/150‑1 (MQ=255)                                     gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTg                                                                                                               <  2:647019/151‑1 (MQ=255)                                            aCTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGGCCGCGGCAGAGGGTTTCCGCCCCTGTCCGCCCTTTCACCACACTTAACTCCTGTTaaa                                                                                                          >  1:264751/1‑149 (MQ=255)                                                 cGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCAGGCAGAGGGTTTCCCCCACTTTGAGCCCTTCCACCCGACTTAACTCCCCATAAAAcgcgc                                                                                                    >  2:174252/1‑150 (MQ=255)                                                   tCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGACTTAACTCCTGATAAACcgcgc                                                                                                    >  2:93894/1‑148 (MQ=255)                                                   tCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCgcgcgc                                                                                                  <  1:688631/150‑1 (MQ=255)                                                    cAACCCCCAGCGCCTATAAGCACCGTTGGCCATCCTGCGCCATTTGCAGGTAATCAAGCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGGGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCa                                                                                                 <  1:968530/150‑1 (MQ=255)                                                        gCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGGGGGAGAGGGTTTTCGCCCCCGTCAGCCCCTCCACCCGACCTAACTCCTGATAAAACCCGCGCAccccc                                                                                            >  1:1073467/1‑151 (MQ=255)                                                          ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCTTCCCCCCGACTTAACTCCTGATAAAACGCGCGCcccccc                                                                                            >  2:171612/1‑149 (MQ=255)                                                                ccTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTCCGCCACTGTCAGCCCTTCCCCCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGc                                                                                    >  2:1072010/1‑151 (MQ=255)                                                                  tGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCAGGGCGAGGGTTTTTGCCCCTGTCAGCCCCTCCCCCAGACTTAAATCCCGGTAAAACCCGCGCCCCCCCCGCTCAGCg                                                                                   >  1:254961/1‑150 (MQ=255)                                                                  tGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCCGGCAGAGGGTTTTCCCCACTGTCAGCCCCTCCCCCCGACTTTACTCCCTATAAAACGCGCGCCCCCCCCCGTCAGCgg                                                                                  >  2:756717/1‑151 (MQ=255)                                                                    tAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCg                                                                                   <  1:933515/148‑1 (MQ=255)                                                                     aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGGAGAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCCGACCTTACTCCTGATAAAACGCCCGCCCCCCCCGCTCAGCgg                                                                                  >  1:227733/1‑148 (MQ=255)                                                                     aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCCCCAGACTTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCAc                                                                               >  2:467362/1‑151 (MQ=255)                                                                         aaCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCATGAACCAAGCAAAGGGTTTTCCCCCCTGTTAGGGCCTCCACTAGGCTTAAAGTCTGATAAACTGCGCGCACCCCCTGGTTAGGGGCAggggg                                                                           >  1:426829/1‑149 (MQ=255)                                                                            gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGGGGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCACCCCACTTTACTCCTGATAAACCCCCCCCACCCCCCCTTCCGCGGCACGGCgcg                                                                         >  2:243501/1‑148 (MQ=255)                                                                            gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTGGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGt                                                                        <  1:297380/151‑1 (MQ=255)                                                                             ttGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGt                                                                        <  1:839739/150‑1 (MQ=255)                                                                             ttGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGCGCCACTTACCGAGGCAAAGGGTTTTCCCCACCGGCAGCCCTTCCACCCGACTTTACTCCTGATAAAACGCGCGCCCCCCCCATTCCGCGGcc                                                                                >  2:346267/1‑141 (MQ=255)                                                                              tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTc                                                                       <  1:147606/150‑1 (MQ=255)                                                                               gggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCa                                                                      <  1:804662/150‑1 (MQ=255)                                                                                ggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGc                                                                    <  1:346260/151‑1 (MQ=255)                                                                                  cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGt                                                                  <  1:516973/151‑1 (MQ=255)                                                                                  cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGt                                                                  <  2:1117045/151‑1 (MQ=255)                                                                                        gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCt                                                              <  2:426298/149‑1 (MQ=255)                                                                                        gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCACAGCAGAGGGGTTCCCCCACTCTAACCCATTCCCCCAGACCTTACTCCTTACAAAACTCGCGCACCCCCGATTCATCGGCACGGCGCGTCAACggtg                                                                >  1:981091/1‑144 (MQ=255)                                                                                         cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCAGGCAGAGGGTTTTCGCCACTGTCCGCCCCTCCACCAGACTTAAATCCTGATAAACCGCGCGCCCCCCCCGCTCAGCGGccccg                                                                             >  1:93933/1‑129 (MQ=255)                                                                                           ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc                                                         <  1:284790/151‑1 (MQ=255)                                                                                            cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCAGAGGGTTTTCCCCACAGTCAGCCCCTCCCCCCGACTTAAATCCTGATAAACCGCGCGCACCCCCCCTTCCGCGGCGCGGCGGGTCAGCGGGGCttcttc                                                         >  1:573002/1‑150 (MQ=255)                                                                                              tttGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGGAGAGGGTTTCTGCCCCCGTCAGCCCCTCCACCCGACCTAACTCCTGATAAAACGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCAGCGTAGCCTCTTcccc                                                      >  1:924657/1‑151 (MQ=255)                                                                                              tttGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGAACGAGGGAGAGGGTTTTCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCCCGCGCACCCCCCGTTCAGCGGCACGGCGAGGCACCGGCGCCTCGTCAcc                                                      >  1:700458/1‑151 (MQ=255)                                                                                                  cAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTcc                                                  <  2:671213/151‑1 (MQ=255)                                                                                                     gTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGt                                               <  1:156776/151‑1 (MQ=255)                                                                                                         tcaatcaTGCCGTTGCGGCGGGGCCACTGAACCCGGGAGAGGGTTTTCGCCCCCGTCCGCCCCTCCACCAGACTTAAATCCCGATAAAACCCGCGCCCCCCCCACTCAGCGGCCCGGCGAGGCCGCGTAGCTTCGTCACCTTCCCGTcgc                                            >  1:350263/1‑150 (MQ=255)                                                                                                          caatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCCGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcgcg                                           <  1:632641/150‑1 (MQ=255)                                                                                                           aatcaTGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCAAACCTTACTCCTGGTAAACCGCGCGCCCCCCCCGTTCAGCGGCACGGCGGGGCAGCGGAGCTTCGTCACCCTCCAGTcgc                                            >  2:1020390/1‑148 (MQ=255)                                                                                                            atcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCCCACTTTACTCCCGATAAAACACGCGCACCCCCCGTTCAGCGGCACCGCGAGGGAGCGGAGCCTCGTTCCCTTCCAGGcgcg                                           >  1:1190503/1‑148 (MQ=255)                                                                                                                tGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCTTCCACCAGACTTAAATCCTGGTAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGAGGCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGcc                                    >  2:199561/1‑151 (MQ=255)                                                                                                                       gcggcgGGGCCACTGCCCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCa                              <  1:257181/150‑1 (MQ=255)                                                                                                                             gggCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATc                        >  2:814250/1‑150 (MQ=255)                                                                                                                             gggCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCa                       <  2:1117263/151‑1 (MQ=255)                                                                                                                              ggCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCAc                      >  2:928250/1‑151 (MQ=255)                                                                                                                                  aCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                  >  1:596550/1‑151 (MQ=255)                                                                                                                                     gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTa                 <  2:1078674/149‑1 (MQ=255)                                                                                                                                     gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTAt                <  1:467351/150‑1 (MQ=255)                                                                                                                                     gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc               >  1:681421/1‑151 (MQ=255)                                                                                                                                              agagGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc      >  1:338018/1‑151 (MQ=255)                                                                                                                                              agagGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc      <  2:215200/151‑1 (MQ=255)                                                                                                                                               gagGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCg     >  1:799210/1‑151 (MQ=255)                                                                                                                                                  gTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATGACCGCTATCCGGGCGTTCGCTg  <  1:839156/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTG  >  NC_000913/2677893‑2678188

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