Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,579,154	A→C	6.3%	A163A (GCA→GCC)	talA →	transaldolase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,154	0	A	C	6.3%	259.9 / 8.2	95	A163A (GCA→GCC)	talA	transaldolase A
Reads supporting (aligned to +/- strand):  ref base A (41/48);  new base C (2/4);  total (43/52)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

CTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGACTACTA  >  NC_000913/2579005‑2579300                                                                                                                                                      |                                                                                                                                                   cTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCa                                                                                                                                                    >  2:240790/1‑150 (MQ=255)   ggtggACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGCAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACg                                                                                                                                                  >  2:942353/1‑150 (MQ=255)     tggACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACggg                                                                                                                                                >  2:1675276/1‑150 (MQ=255)     tggACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACggg                                                                                                                                                >  2:368045/1‑150 (MQ=255)       gACTTGTATCAGCAACAAGTCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGTAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGcc                                                                                                                                              <  1:182178/150‑1 (MQ=255)       gACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGcc                                                                                                                                              <  2:390228/150‑1 (MQ=255)           tGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCtgtg                                                                                                                                          <  1:1256035/150‑1 (MQ=255)           tGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCtgtg                                                                                                                                          >  1:470220/1‑150 (MQ=255)            ttATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGc                                                                                                                                         <  1:542046/149‑1 (MQ=255)            gTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGc                                                                                                                                         >  2:961471/1‑150 (MQ=255)             tATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCg                                                                                                                                        >  1:1109564/1‑150 (MQ=255)               tCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaa                                                                                                                                      <  2:921393/150‑4 (MQ=255)               tCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaa                                                                                                                                      <  1:1064695/150‑4 (MQ=255)               tCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGa                                                                                                                                      >  2:1044063/1‑150 (MQ=255)                 aGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAg                                                                                                                                    >  1:200702/1‑150 (MQ=255)                 aGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGATGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGGTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAg                                                                                                                                    >  1:35168/1‑150 (MQ=255)                  gCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGc                                                                                                                                   <  2:324423/150‑1 (MQ=255)                     acaaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAgg                                                                                                                                <  2:252802/150‑1 (MQ=255)                      caaGGCGTTGAGAAATCACGCATTCTGATCACGCTGGCTTCGACCTGGGAAGGAATGCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGc                                                                                                                               <  1:191445/150‑1 (MQ=255)                      caaGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGc                                                                                                                               >  1:981710/1‑150 (MQ=255)                         ggCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGACGGTAGTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGtt                                                                                                                            <  1:1087967/150‑1 (MQ=255)                         ggCGGTGAGAAAGCACGCCTTCTGCTCAAGCGGGCTTCGACCGGGGCAGGACTTCGCGCGGCAGAAGAGCGGGAAAACGAAGGTATTACCTGCAACCTGACGCTGCGGGGTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGtt                                                                                                                            <  2:1278544/150‑1 (MQ=255)                            gTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGttttt                                                                                                                         <  2:200702/150‑1 (MQ=255)                              tGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCt                                                                                                                       >  1:349448/1‑150 (MQ=255)                               gagaAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTg                                                                                                                      <  1:1010605/150‑1 (MQ=255)                               gagaAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTg                                                                                                                      <  2:349448/150‑1 (MQ=255)                                agaAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGa                                                                                                                     >  1:1196385/1‑150 (MQ=255)                                      cACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGc                                                                                                               >  1:1488226/1‑150 (MQ=255)                                      cACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGc                                                                                                               >  2:809358/1‑150 (MQ=255)                                       gcGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGcc                                                                                                              >  1:1417200/2‑150 (MQ=255)                                          cATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGtt                                                                                                           >  1:1577115/1‑150 (MQ=255)                                          cATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCGGTTTTCTTTTGCACAGGTACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGtt                                                                                                           <  1:732204/150‑1 (MQ=255)                                          cATTCTGATAAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCTACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGtt                                                                                                           >  1:1207040/1‑150 (MQ=255)                                            ttCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTg                                                                                                         >  2:1480817/1‑150 (MQ=255)                                                gATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCgg                                                                                                     >  2:311630/1‑150 (MQ=255)                                                 aTCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCggg                                                                                                    <  2:1089137/150‑1 (MQ=255)                                                 aTCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCggg                                                                                                    <  1:691439/150‑1 (MQ=255)                                                  tCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGc                                                                                                   >  2:1409434/1‑150 (MQ=255)                                                   cAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTTGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCATGCGTTTTTCTGATTTCGCCGTTTGTCGGGCg                                                                                                  <  2:903057/150‑1 (MQ=255)                                                   cAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGTAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCg                                                                                                  <  2:909194/150‑1 (MQ=255)                                                       cTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTAtt                                                                                              >  1:235087/1‑150 (MQ=255)                                                         ggCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTa                                                                                            <  1:1409434/150‑1 (MQ=255)                                                             tCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGAc                                                                                        >  2:1233409/1‑150 (MQ=255)                                                                  cTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTa                                                                                   <  2:470220/150‑1 (MQ=255)                                                                     ggaaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCa                                                                                <  1:85132/150‑1 (MQ=255)                                                                       aaggaaTTCGCGCTGCAGAAGAGCTGGAAAAAGAAGTTATTAACTGCCACCTGACGCTGCTGTTTTCTTTTGCACAGTCACGGGCCTGTGCGGAAGCAGGCTTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAgg                                                                              <  2:1093442/150‑1 (MQ=255)                                                                        aggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGATTGTCGGGCGTATTTATGACTGGTATCAGGc                                                                             >  2:872174/1‑150 (MQ=255)                                                                         ggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCa                                                                            <  2:1441381/150‑1 (MQ=255)                                                                           aaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCCCAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACg                                                                          <  1:698814/150‑1 (MQ=255)                                                                           aaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACg                                                                          >  1:559134/1‑150 (MQ=255)                                                                               cgcgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGTGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAg                                                                      <  2:429335/150‑1 (MQ=255)                                                                               cgcgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAg                                                                      <  1:1675276/150‑1 (MQ=255)                                                                               cgcgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAg                                                                      <  1:1044063/150‑1 (MQ=255)                                                                               cgcgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAg                                                                      >  2:144127/1‑150 (MQ=255)                                                                               cgcgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAg                                                                      <  1:368045/150‑1 (MQ=255)                                                                               cgcgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGAGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAg                                                                      >  2:833323/1‑150 (MQ=255)                                                                                 cgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTa                                                                                            <  2:754227/126‑1 (MQ=255)                                                                                 cgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTa                                                                                            >  1:754227/1‑126 (MQ=255)                                                                                         aagaTCTGGAAAACGTATGTATTATCTGCTACCAGAGGCTGCGGTTTTCTTTTGCACAGGCACGGGACTGTGCGGATGGAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTTTTTATTAATGGTTTCAGGGACGCAAGCCGATggtgc                                                            <  2:733041/150‑4 (MQ=255)                                                                                          agagCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAccc                                                           >  1:1293032/1‑150 (MQ=255)                                                                                           gcgcaggAAAAAGAACGTATTAACTGCAACCTGACGCTTCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCATGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCg                                                          <  2:216494/145‑1 (MQ=255)                                                                                           gagCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATgg                                                               >  1:1580253/1‑145 (MQ=255)                                                                                           gagCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATgg                                                               <  2:1580253/145‑1 (MQ=255)                                                                                               tGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGAGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGAT‑GACCCGTAtgt                                                     <  1:833323/150‑1 (MQ=255)                                                                                                  aaaaaGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATgtgg                                                   <  1:1569568/150‑1 (MQ=255)                                                                                                     aagaagGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGtggtgg                                                <  1:961471/150‑1 (MQ=255)                                                                                                      agaagGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGga                                               >  1:325215/1‑150 (MQ=255)                                                                                                      agaagGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGga                                               <  1:240790/150‑1 (MQ=255)                                                                                                         agGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaaga                                            <  1:1233409/150‑1 (MQ=255)                                                                                                          ggTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGgaagaa                                           <  1:311630/150‑1 (MQ=255)                                                                                                            tataAACTTCAACCTGACGCTGCTGTTTTCTTTTGCACAGTCACGGGCCTGTTCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAagaaga                                         <  1:1079194/150‑1 (MQ=255)                                                                                                               tAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATcc                                      <  1:740913/150‑1 (MQ=255)                                                                                                                 aCTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCgg                                    >  1:1051069/1‑150 (MQ=255)                                                                                                                  cTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCggg                                   >  2:1074085/1‑150 (MQ=255)                                                                                                                   tGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGc                                  <  2:1577115/150‑1 (MQ=255)                                                                                                                   tGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGc                                  <  1:1480817/150‑1 (MQ=255)                                                                                                                     cAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGt                                <  2:207991/150‑1 (MQ=255)                                                                                                                      aaCCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGtt                               <  2:235087/150‑1 (MQ=255)                                                                                                                       aCCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTa                              <  2:1196385/150‑1 (MQ=255)                                                                                                                        ccTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTaa                             <  1:1517824/150‑1 (MQ=255)                                                                                                                             cGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCgg                        <  2:559134/150‑1 (MQ=255)                                                                                                                                 gctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTgcg                    >  2:909029/1‑150 (MQ=255)                                                                                                                                    gTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCaa                 >  2:480650/1‑150 (MQ=255)                                                                                                                                     tgttCTTTTGCACAGGCACGGGCCTGTGTGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCTTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAat                <  1:134539/148‑1 (MQ=255)                                                                                                                                         cTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCt            >  1:1254431/1‑150 (MQ=255)                                                                                                                                         cTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCt            >  2:1442189/1‑150 (MQ=255)                                                                                                                                           tttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGTAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTAc          >  1:1046023/1‑150 (MQ=255)                                                                                                                                           tttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCTGGCGTTAAATCGGTGCGCAATATCTAc          <  1:144127/150‑1 (MQ=255)                                                                                                                                           tttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTAc          >  1:1582929/1‑150 (MQ=255)                                                                                                                                            ttGGACAGGCACGGGCCGGTGCGGATGCAGGCGTTTTTCTGAATTCGCCGGTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACg         <  1:20055/150‑1 (MQ=255)                                                                                                                                            ttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACg         <  1:441772/150‑1 (MQ=255)                                                                                                                                            ccgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACg         >  2:1560861/3‑150 (MQ=255)                                                                                                                                             tGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCCTATTTATGACTGGTATCAGGCACGCAAGCCGAAGGACCCGTATGTGGTGGAAGAACATCCGGGCGTTAAATCGGTGCGCAATATCTACGa        >  1:805076/1‑150 (MQ=255)                                                                                                                                             cgCACAGGCCCGTGCTTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGa        >  2:1691801/2‑150 (MQ=255)                                                                                                                                                 caGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCGGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCTATATCTACGactac    >  2:131563/1‑150 (MQ=255)                                                                                                                                                   ggCACGGGCCTGTGCGGAAGCAGGAGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGActacta  >  2:1214675/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGAAGAAGATCCGGGCGTTAAATCGGTGCGCAATATCTACGACTACTA  >  NC_000913/2579005‑2579300

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