Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,234,331	T→G	7.8%	K319Q (AAG→CAG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,234,331	0	T	G	7.8%	235.5 / 9.3	77	K319Q (AAG→CAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base G (4/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTTCATCAA  >  NC_002947/1234183‑1234478                                                                                                                                                      |                                                                                                                                                    ccACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCtt                                                                                                                                                     <  2:176462/149‑1 (MQ=255)  cACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCtt                                                                                                                                                     >  2:316533/1‑148 (MQ=255)     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                                                                                                              <  5:67458/149‑1 (MQ=255)                                        tCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCAtt                                                                                                              <  8:269485/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                                             <  7:22508/149‑1 (MQ=255)                                                                                                          cTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCa                                             <  3:235279/149‑1 (MQ=255)                                                                                                              cTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGc                                          <  6:121118/148‑1 (MQ=255)                                                                                                               tGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAg                                                                                            <  7:349840/97‑1 (MQ=255)                                                                                                               tGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAg                                                                                            >  8:349840/1‑97 (MQ=255)                                                                                                               tGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAACCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCg                                        >  7:178243/1‑149 (MQ=255)                                                                                                                 gCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGt                                                  >  1:289518/1‑137 (MQ=255)                                                                                                                 gCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGt                                                  <  2:289518/137‑1 (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)                                                                                                                                                      |                                                                                                                                                    CCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGGCTGGCACCAGGTA‑CACGCCATTGCTGTCCTTGACCTTGCTGGCGAAATACTCGGTATCCAGGGCGTCGTTGACGATCTTCAGTTCGTCACGCAGCCATTGCACGGTAGAGCCACCATTGAATACCGCGCCTTCCAGTGCATAGGCCACCTCACCACGCGGGCCGCAGGCAATGGTCGTGAGCAGGCCGTGGGACGACTTCACCGCCTGGTCGCCCGTGTTCATCAA  >  NC_002947/1234183‑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: ‑