Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	8,588	A→G	75.7%	K117K (AAA→AAG)	talB →	transaldolase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	8,588	0	A	G	75.7%	68.3 / 19.5	37	K117K (AAA→AAG)	talB	transaldolase B
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base G (13/15);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAG  >  CP009273/8442‑8726                                                                                                                                                   |                                                                                                                                           cagATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGctg                                                                                                                                         <  2:47658/150‑2 (MQ=255)  agATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAActct                                                                                                                                        <  2:106107/150‑1 (MQ=255)          gatgcttccgatAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGa                                                                                                                                 <  1:54146/142‑1 (MQ=255)             gcttccgatAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGcg                                                                                                                             >  2:8436/5‑149 (MQ=255)                     taagcTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAg                                                                                                                     <  1:175660/146‑1 (MQ=255)                                    catcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCt                                                                                                      >  1:85418/5‑150 (MQ=255)                                      tcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGa                                                                                                    <  1:126184/148‑1 (MQ=255)                                                   ccTGAAACTGGTTCCGGGCCGGATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTc                                                                                       >  1:87062/1‑150 (MQ=255)                                                        aaCTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCt                                                                                  >  1:166432/1‑150 (MQ=255)                                                         aCTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTg                                                                                 <  2:34287/150‑1 (MQ=255)                                                         aCTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTg                                                                                 <  2:156867/150‑1 (MQ=255)                                                         aCTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTg                                                                                 <  2:136966/150‑1 (MQ=255)                                                         aCTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTg                                                                                 >  2:163992/1‑150 (MQ=255)                                                                      cGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGt                                                                    <  2:56898/150‑1 (MQ=255)                                                                             ccaccgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGcgccgct                                                              >  1:26209/6‑148 (MQ=255)                                                                                     gTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTg                                                     >  1:159872/1‑150 (MQ=255)                                                                                         aTGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGaa                                                  >  2:158248/1‑149 (MQ=255)                                                                                            cacgTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGaaa                                                 >  2:169713/3‑147 (MQ=255)                                                                                              cgTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAaagaa                                            <  1:163992/150‑1 (MQ=255)                                                                                               gTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGATGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCCCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAaagaag                                           <  1:86026/150‑1 (MQ=255)                                                                                                   ctCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAt                                       >  2:42370/2‑150 (MQ=255)                                                                                                      cTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTaa                                    <  2:47724/150‑1 (MQ=255)                                                                                                         tGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAAc                                   <  1:144326/148‑1 (MQ=255)                                                                                                           acacCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGt                                <  1:58100/149‑1 (MQ=255)                                                                                                           acacCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTa                               <  2:85418/150‑1 (MQ=255)                                                                                                             cccGAAGCGTCTATCGCGAAAGCCACGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAAc                             <  1:61782/149‑1 (MQ=255)                                                                                                                gAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCt                           >  2:97631/1‑149 (MQ=255)                                                                                                                gAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCt                           >  1:15725/1‑149 (MQ=255)                                                                                                                  aGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCCGAc                        <  1:167412/150‑1 (MQ=255)                                                                                                                   gCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGAcg                       >  1:59341/1‑149 (MQ=255)                                                                                                                     gTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGct                     >  1:140583/1‑150 (MQ=255)                                                                                                                      tctatcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctg                    >  1:19050/7‑150 (MQ=255)                                                                                                                            gCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGtt               <  2:71850/149‑1 (MQ=255)                                                                                                                            gCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTtc              >  2:38708/1‑150 (MQ=255)                                                                                                                             cGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtct             >  2:134773/1‑150 (MQ=255)                                                                                                                             cGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTtc              <  2:159872/149‑1 (MQ=255)                                                                                                                                        cgcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAg  >  1:38725/7‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                           CAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAG  >  CP009273/8442‑8726

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