Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,234,332	C→A	8.5%	L318L (CTG→CTT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,234,332	0	C	A	8.5%	243.5 / 7.1	78	L318L (CTG→CTT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base C (35/36);  new base A (5/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTTCATCAA  >  NC_002947/1234187‑1234478                                                                                                                                                   |                                                                                                                                                   gtgtCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGt                                                                                                                                                 <  4:292888/149‑1 (MQ=255)  tgtCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGtt                                                                                                                                                >  6:305877/1‑149 (MQ=255)      aGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGt                                                                                                                                             <  3:258397/148‑1 (MQ=255)        gCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCAc                                                                                                                                          <  5:273942/149‑1 (MQ=255)        gCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTCGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCAc                                                                                                                                          >  4:80893/1‑149 (MQ=255)          cGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGc                                                                                                                                        <  5:365090/149‑1 (MQ=255)              caacGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGGCGGTAAAGGCTGGCTTGAGGTACCACGCTTTTGCTGGCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGc                                                                                                                                     <  1:354734/149‑1 (MQ=255)               aacGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCa                                                                                                                                   >  6:435024/1‑149 (MQ=255)                acGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCAt                                                                                                                                  >  6:67458/1‑149 (MQ=255)                acGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCAt                                                                                                                                  <  5:197263/149‑1 (MQ=255)                    accacGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCa                                                                                                                              >  2:187306/1‑149 (MQ=255)                      cacGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACg                                                                                                                            >  5:393470/1‑149 (MQ=255)                      cacGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACg                                                                                                                            >  5:121118/1‑149 (MQ=255)                      cacGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGAACTGAAGATCGTCaac                                                                                                                                         >  6:134859/1‑134 (MQ=255)                      cacGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGAACTGAAGATCGTCaac                                                                                                                                         <  5:134859/136‑3 (MQ=255)                       acGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCg                                                                                                                                              <  4:133116/130‑1 (MQ=255)                       acGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACgg                                                                                                                           >  6:174551/1‑149 (MQ=255)                           gCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTa                                                                                                                         >  8:22508/1‑147 (MQ=255)                             aTACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGc                                                                                                                     >  4:235279/1‑149 (MQ=255)                               aCGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGATGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACCAAGCCATTTCCCGTAAGAGcca                                                                                                                   >  8:114371/1‑149 (MQ=255)                                cGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGccac                                                                                                                  <  1:187306/149‑1 (MQ=255)                                  ggTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGccacca                                                                                                                <  2:29200/149‑1 (MQ=255)                                    tCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCAtt                                                                                                              <  8:269485/149‑1 (MQ=255)                                    tCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCAtt                                                                                                              <  5:67458/149‑1 (MQ=255)                                      ccAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTg                                                                                                             <  4:32293/148‑1 (MQ=255)                                       cAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGaa                                                                                                           >  4:271771/1‑149 (MQ=255)                                          tACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGCCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATAc                                                                                                        >  3:386455/1‑149 (MQ=255)                                              gggCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACcgcg                                                                                                    <  7:99376/149‑1 (MQ=255)                                               ggCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGcc                                                                                                                                    <  8:380635/116‑1 (MQ=255)                                               ggCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGcc                                                                                                                                    >  7:380635/1‑116 (MQ=255)                                                   cccAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTc                                                                                               <  5:435024/149‑1 (MQ=255)                                                    ccAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTc                                                                                               >  3:156363/1‑148 (MQ=255)                                                     cAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCa                                                                                             <  2:301821/149‑1 (MQ=255)                                                          cGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCa                                                                                        <  7:114371/149‑1 (MQ=255)                                                             tAAAGGCTGGCAGCAGGTA‑CACGCCATTGCTGTACTTGACCTTTCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAg                                                                                     <  8:145602/149‑1 (MQ=255)                                                             tAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAg                                                                                     >  3:71646/1‑149 (MQ=255)                                                               aaGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGc                                                                                   <  2:80750/149‑1 (MQ=255)                                                                  gCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCAc                                                                                <  6:393470/149‑1 (MQ=255)                                                                       cACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTcac                                                                           <  4:156363/149‑1 (MQ=255)                                                                       cACCAGGTA‑CACGACATTGCTGTCCTTGCCCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTTACGATCTTCAGTTCGTCCCGCAGGAATTGCACGGTAGAGACACAATTGAATACCGCGCCTTCCAGTGCATAGGCGACATcac                                                                           >  8:213468/1‑149 (MQ=255)                                                                        aCCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGATGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTAAGTTCGTCAAGCAGCCATTGCACTGTAGAGACACCATTGAATAGCGCGCCTTCCAGTGCATAGGCAACCTcac                                                                           >  2:339234/1‑148 (MQ=255)                                                                           aGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCAcg                                                                       >  4:394603/1‑149 (MQ=255)                                                                           aGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCAcg                                                                       >  1:100978/1‑149 (MQ=255)                                                                                  acGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGcc                                                                 <  3:80893/149‑1 (MQ=255)                                                                                      cATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAg                                                             <  1:316533/149‑1 (MQ=255)                                                                                      cATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAg                                                             >  4:140569/1‑149 (MQ=255)                                                                                                 ttGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGt                                                  >  6:255485/1‑149 (MQ=255)                                                                                                    aCCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGa                                                <  5:174551/148‑1 (MQ=255)                                                                                                    aCCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAg                                               <  6:159649/149‑1 (MQ=255)                                                                                                      cTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCa                                             <  3:235279/149‑1 (MQ=255)                                                                                                      cTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCa                                             <  7:22508/149‑1 (MQ=255)                                                                                                          cTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGc                                          <  6:121118/148‑1 (MQ=255)                                                                                                           tGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAg                                                                                            >  8:349840/1‑97 (MQ=255)                                                                                                           tGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAg                                                                                            <  7:349840/97‑1 (MQ=255)                                                                                                           tGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAACCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCg                                        >  7:178243/1‑149 (MQ=255)                                                                                                             gCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGt                                                  <  2:289518/137‑1 (MQ=255)                                                                                                             gCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGt                                                  >  1:289518/1‑137 (MQ=255)                                                                                                              cGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGt                                                  <  4:57580/136‑1 (MQ=255)                                                                                                              cGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTgg                                     <  5:305877/149‑1 (MQ=255)                                                                                                              cGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGACGCAGGCAATGGTCGt                                                  >  3:57580/1‑136 (MQ=255)                                                                                                               gAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCTGGCCGTgg                                     >  3:270645/1‑148 (MQ=255)                                                                                                               gAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTgg                                     >  8:51460/1‑148 (MQ=255)                                                                                                                    aCTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGTGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGgacgac                               >  8:241136/1‑149 (MQ=255)                                                                                                                     cTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGgacgac                               <  5:255485/148‑1 (MQ=255)                                                                                                                      tCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACt                              >  4:258274/1‑148 (MQ=255)                                                                                                                          tATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCTGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCAcc                         <  1:265069/149‑1 (MQ=255)                                                                                                                          tATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCAcc                         <  3:258274/149‑1 (MQ=255)                                                                                                                          tATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCAcc                         <  3:271771/149‑1 (MQ=255)                                                                                                                             ccAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGcc                      >  5:279060/1‑149 (MQ=255)                                                                                                                               aGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTg                    >  2:276580/1‑149 (MQ=255)                                                                                                                                     tgcgtGACGAACTGAAGATCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGcc                      >  6:351132/5‑141 (MQ=255)                                                                                                                                     tgcgtGACGAACTGAAGATCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGAAGGCCGTGGGACGACTTCACCGcc                      <  5:351132/137‑1 (MQ=255)                                                                                                                                        ttGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCgt           >  5:43256/1‑149 (MQ=255)                                                                                                                                        gtGACGAACTGAAGATCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCgt           >  4:243425/2‑149 (MQ=255)                                                                                                                                        gtGACGAACTGAAGATCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCgt           >  4:237849/2‑149 (MQ=255)                                                                                                                                           aCGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGtt        >  3:16256/1‑149 (MQ=255)                                                                                                                                             gATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTtca      >  3:420420/1‑149 (MQ=255)                                                                                                                                                 ttCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTTCATCaa  >  3:340516/1‑149 (MQ=255)                                                                                                                                                   cAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTtcatc    <  3:140569/145‑1 (MQ=255)                                                                                                                                                   cAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCACAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTtcat     <  3:394603/144‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTTCATCAA  >  NC_002947/1234187‑1234478

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