Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,678,043	0	C	T	41.8%	80.0 / 27.7	67	E107K (GAA→AAA)	yphE	putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (7/32);  new base T (28/0);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.15e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677903‑2678187                                                                                                                                             |                                                                                                                                                  cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTc                                                                                                                                         <  2:242136/151‑1 (MQ=255)  gcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGGCCCAGGCAGAGGGGTTTCCCCACTGt                                                                                                                                          >  1:1149027/1‑149 (MQ=255)   ccGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAg                                                                                                                                       <  2:402192/151‑1 (MQ=255)    cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑c                                                                                                                                     <  2:339636/151‑1 (MQ=255)      ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCCGGCCGAGGGTTTTCGCCCCCGTCCG‑ccc                                                                                                                                   >  1:1579736/1‑151 (MQ=255)           aGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTcc                                                                                                                               <  2:57328/150‑1 (MQ=255)             cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAG‑CCCTTccacc                                                                                                                            >  2:178646/1‑151 (MQ=255)             cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTTACCGAGGCAGAGGGTTTTCGCCACTGTCCG‑CCCCTccccc                                                                                                                            >  2:221020/1‑151 (MQ=255)               ttGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCCCCCCCGTCAG‑CCCTTccccc                                                                                                                            >  2:818019/1‑149 (MQ=255)                tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAg                                                                                                                          <  2:251603/150‑1 (MQ=255)                  aaaCAAGTTGTTCAGGCCTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAg                                                                                                                          <  1:907814/148‑1 (MQ=255)                    aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCAGAGGGTTTTCGCCCCTGTCAG‑CCCCTCCACCCGACtt                                                                                                                      >  2:1033902/1‑150 (MQ=255)                    aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTGCCGCCACTGTCAG‑CCCTTCCACCAGACTTa                                                                                                                     <  1:300584/151‑1 (MQ=255)                     cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTaa                                                                                                                    <  2:1594258/151‑1 (MQ=255)                      aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTaa                                                                                                                    <  2:1231217/150‑1 (MQ=255)                       aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACt                                                                                                                  <  1:1127095/151‑1 (MQ=255)                       aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAG‑CCCTTCCCCCAGACTTaa                                                                                                                    >  2:1070877/1‑149 (MQ=255)                       aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACCGTCAG‑CCCTTCCCCCAGACttta                                                                                                                    >  2:129421/1‑147 (MQ=255)                               aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAATGACCGAGGGAGAGGGTTTTCGCCCCTGTCAG‑CCCTTCCACCAGACCTAACTCCTGa                                                                                                             >  2:834861/1‑148 (MQ=255)                               aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGCAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCAGAGGGTTTCCGCCCCTGTCAG‑CCCCTCCCCCCCACTTAACTCCTGAt                                                                                                            >  1:1205727/1‑149 (MQ=255)                                ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATaaa                                                                                                         <  1:1045572/151‑1 (MQ=255)                                    tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCGCCACTGTCAG‑CCCCTCCCCCACACTTAACTCCCGATaaa                                                                                                         >  2:490458/1‑147 (MQ=255)                                     aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACcgcg                                                                                                    <  2:1312540/151‑1 (MQ=255)                                     aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCGGAGGGTTTTCGCCCCTGCTAG‑CCCCTCCCCCCGCCTTAACACCTGATaaa                                                                                                         >  1:313073/1‑146 (MQ=255)                                        gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCAGGCAGAGGGTTTTCGCCCCCGTCAG‑CCCCTCCACCCGACTTAAATCCCGATAAACcgcgc                                                                                                   >  1:1672602/1‑149 (MQ=255)                                           aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGAACGCGGCAGAGGGTTTTCGCCCCCGTCAA‑CCCTTCCACCAACCTTAAATCCTCATAAAACCCCCGCAcc                                                                                              >  1:1103327/1‑151 (MQ=255)                                                ccAGCGCCTGTAAGCAGCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCAGGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAg                                                                                         <  1:1676357/151‑1 (MQ=255)                                                ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGAGCCAGGCAGAGGGTTTCCCGCCCTGTCAA‑CCCCTCCACCCGAATTAAATCCCGATAAAACGCGCGCAccccc                                                                                           >  2:914769/1‑149 (MQ=255)                                                  aGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCCCCGTCAG‑CCCTTCCACCAGACTTAAATCCTGATAAACCGCGCGCAccccccggt                                                                                       >  2:614091/1‑147 (MQ=255)                                                       cTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCCGAGGGTTTTCGCCACTGTCAG‑CCCTTCCCCCAGACTTAACTCCTGATAAACCCCGCGCCCCCCCAGTTCAGCg                                                                                  >  2:1490226/1‑151 (MQ=255)                                                        tGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGGGGCAGGGGGTTTCCGCCCCCGTCAG‑CCCTTCCCCCCCACCTTACTCCTGGTAAAACGCGCGCccccccccttccac                                                                                   >  1:444914/1‑141 (MQ=255)                                                           aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCAc                                                                              <  1:26246/151‑1 (MQ=255)                                                           aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCAc                                                                              <  1:1590550/151‑1 (MQ=255)                                                           aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCTACCCCGGGCGAGGGTTTTCGCCCCTGTCAG‑CCCCTCCACCCCACCTAAATCCTGCTAAAACCCGCccccccccccgt                                                                                       >  1:1672542/1‑138 (MQ=255)                                                            aGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGGCCCAGGCAGAGGGTTTTCGCCACTGTCAG‑CCCCTCCCCCACACTTAACTCCTGATAAACCGCGCGCCCCCCCCGCTCAGCGGccc                                                                              >  1:556490/1‑148 (MQ=255)                                                              cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACgg                                                                            <  1:602762/150‑1 (MQ=255)                                                              cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAAGCAGAGGGGTTTCCCCCCGCACCG‑CCCTTCCCCCCGACTTAAATCCTGATAAAACGCGCGCCCCCCCCGGTCacccggccg                                                                             >  2:1647416/1‑141 (MQ=255)                                                               aaCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAAACGCGCGCACCCCCCGTTCAGCgg                                                                                 >  2:984719/1‑144 (MQ=255)                                                                    tGGGCATCCTGCGCCATTTGCAGGAAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTc                                                                      <  2:1277380/150‑1 (MQ=255)                                                                        cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGCGGCAGAGGGTTTTCCCCACCGTCCG‑CCCCTCCCCCCCACTTAAATCCTTATAAACCCCCCCCACCCCCCGTTCAACGGCACCGCGAGTCacc                                                                   >  2:822968/1‑147 (MQ=255)                                                                          tCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGCCCCAGGCCGAGGTTTTTCCCCACCGTCAG‑CCCCTCCACCCGACTTAACTCCTGGTAAACCCCGCGCACCCCCCCTTCAGCGGCCCCGCGAGTCAGCg                                                                  >  1:302499/1‑148 (MQ=255)                                                                             tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCGCCCCTGTCAG‑CCCTTCCACCCGACTTAAATCCTGATAAAACGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCagcttagctt                                                            >  2:131281/1‑151 (MQ=255)                                                                             tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCCCCGTCAG‑CCCTTCCACCCCACTTAACTCCCGATAAAACGCCCGCACCCCCCCTTCAGCGGCACCGCGGGTCAGCGCAGc                                                              >  1:1096144/1‑149 (MQ=255)                                                                                     ttGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCCGAGGGTTTTCGCCCCTGTCCGTCCCTTCC‑CCACACTTAACACCTGTTAAACCGCCCCCACCCCCCCGTCAGCGGCACGGCGGGGCACGGGAGCTTCGTCAcc                                                     >  1:176941/1‑150 (MQ=255)                                                                                       gCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCtt                                                   <  1:284158/150‑1 (MQ=255)                                                                                         aGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCTAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCa                                                <  1:659934/151‑1 (MQ=255)                                                                                             aatcaatcATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcg                                            <  1:1007672/151‑1 (MQ=255)                                                                                             aatcaatcATGCCGTTGCGGCGGGGCCACTGACCCAGGCAGAGGGTTTTCCCCCCCCTCAG‑ACCTTCCCCCCGACTTAAATCCTGATAAAACCCGCGCACCCCCCCTTCAACGGGACGGCGAGGCAACGTAGCCTCGTCACCCTcccgcc                                             >  1:1233782/1‑146 (MQ=255)                                                                                                  atcaTGCCGTTGCGGCGGGGCCACTGACCCAGGGAGAGGGTTTTCGCCCCTGTCAG‑CCCCTCCCCCCCACTTAACTCCTGATAAAACCCGCGCACCCCCCGTTCAGCGGCACCGGGGGGCAGCGGAGCTTCGTCACCCTcccctc                                             >  1:1148811/1‑141 (MQ=255)                                                                                                   tcaTGCCGTTGCGGCGGGGCCCCAGACCGAGGGAGAGGGTTTTCCCCCCTGTCAG‑CCCCTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCCCTTCAGCGGCACGGCGAGGCAGCGGAGCTTCGTCACCCTCCCGTCGCGgctc                                      >  1:887460/1‑151 (MQ=255)                                                                                                    caTGCCGTTGCGGCGGGGCCCCCGACCCGGGCAGAGGGTTTTCGCCCCCGTCCG‑CCCCTCCCCCCGACTTAACTCCTGATAAAACGCGCGCACCCCCCGGTCAGCGGGACGGCGAGTCAGCGGAGCTTCGTCACCCTCCCGt                                              >  1:421638/1‑142 (MQ=255)                                                                                                    caTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGGCTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCg                                     <  1:948509/151‑1 (MQ=255)                                                                                                    caTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtctc                                      <  1:1329503/150‑1 (MQ=255)                                                                                                     aTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtctc                                      <  1:59480/149‑1 (MQ=255)                                                                                                      tGCCGTTGCGGCGGGGCCACTCACCGAGGGAGAGGGTTTCCGCCACTGTCCG‑CCCCTCCCCCAGACTTAAATCCTGATTAACCACCCGCACCCCCAGCTCTGCGGCACGGCGGGGCAGCCGAGCTTCCTCACCTTCCAGTCCCGTCCCGcc                                   >  1:23328/1‑151 (MQ=255)                                                                                                              cggcggGGCCCCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAAACGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCCCCTTCCAGTCGCGTCTCGCCAATcc                              >  2:530247/1‑148 (MQ=255)                                                                                                                 cggGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                        >  1:1418507/1‑151 (MQ=255)                                                                                                                         cTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                 <  1:656533/150‑1 (MQ=255)                                                                                                                           gACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc              <  1:796730/151‑1 (MQ=255)                                                                                                                             ccGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg            <  1:921527/151‑1 (MQ=255)                                                                                                                              cGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCgg           <  1:723496/151‑1 (MQ=255)                                                                                                                                  gCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGt       >  2:453465/1‑151 (MQ=255)                                                                                                                                     gggattttCCGCCACTGTCAG‑CCCCGCCGCCAGACTTAACTCCGGATAAGCCGCGCACACCCCCAGTTCAGCGGCACGGGGAGTTAGCGTAGCTTCGTCACCTTCCAGTGGCGTCTCGCCAATCCAGACATCACCGCTATCCGGGCGtt      <  1:1185999/145‑1 (MQ=255)                                                                                                                                     gagTTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc     <  2:542621/150‑1 (MQ=255)                                                                                                                                       ggTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  <  2:4346/151‑1 (MQ=255)                                                                                                                                       ggTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  <  1:746049/151‑1 (MQ=255)                                                                                                                                       ggTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  <  2:928715/151‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                  CGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677903‑2678187

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