Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,995,991	G→T	10.8%	P46P (CCC→CCA)	PP_4399 ←	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,995,991	0	G	T	10.8%	119.3 / 3.3	37	P46P (CCC→CCA)	PP_4399	glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base G (17/16);  new base T (2/2);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AGCTGGCCCCGTGTGGACAGTGCCGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAGCCTTCG  >  NC_002947/4995884‑4996120                                                                                                            |                                                                                                                                  aGCTGGCCCCGTGTGGACAGTGCCGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCcgt                                                                                                       >  2:216112/1‑136 (MQ=255)              tgGACAGTGCCGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTa                                                                                          >  1:374967/1‑136 (MQ=255)                        cGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTTGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCt                                                                                >  1:301436/1‑136 (MQ=255)                        cGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTTGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCt                                                                                >  1:340383/1‑136 (MQ=255)                        cGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAAttt                                                                                    >  2:339437/1‑132 (MQ=255)                        cGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAAttt                                                                                    <  1:339437/132‑1 (MQ=255)                        cGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCt                                                                                >  2:193855/1‑136 (MQ=255)                        cGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCt                                                                                >  1:197098/1‑136 (MQ=255)                        cGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCAAGGGAAGGATGTGGAAATTACTGGGTTCGGCCACCAGTGCCCGGTCGCCGTCGTCACTGCCGTAGAATTTCGCt                                                                                >  1:244459/1‑136 (MQ=255)                        cGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCACCACGGCACGCGGTGGACCAGACTGGGTTCGGCCACCAGTGCCAGGTAGCCGTCCTCACTGACGTAGAATTTCGCt                                                                                >  1:279426/1‑136 (MQ=255)                              ccggccATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGa                                                                          <  1:88447/136‑1 (MQ=255)                                  ccATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCCTCGCGCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCCCCAGTGCCCGGTCGCCGTCGTCACAGCCGTAGAATTTCGCTGGTGGATaaa                                                                      >  1:384580/1‑135 (MQ=255)                                  ccATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTATAATTTCGCTGGTGGATaaa                                                                      >  1:102462/1‑135 (MQ=255)                                  ccATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAAc                                                                      >  2:320703/1‑136 (MQ=255)                                  ccATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAAc                                                                      >  2:195945/1‑136 (MQ=255)                                     tCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGc                                                                   <  1:320703/136‑1 (MQ=255)                                     tCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGc                                                                   <  1:422397/136‑1 (MQ=255)                                          cTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCatg                                                              <  2:69536/136‑1 (MQ=255)                                          cTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGGGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCatg                                                              <  2:46695/136‑1 (MQ=255)                                             gTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCatgatg                                                           <  2:335559/136‑1 (MQ=255)                                                     aacaGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTgcag                                                   <  2:339109/135‑1 (MQ=255)                                                                           gTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTgg                             >  2:105612/1‑135 (MQ=255)                                                                                         gCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTc               <  1:10731/136‑1 (MQ=255)                                                                                          cAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCg              <  1:79956/136‑1 (MQ=255)                                                                                          cAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCg              <  1:277496/136‑1 (MQ=255)                                                                                            ggCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAg            <  2:239302/136‑1 (MQ=255)                                                                                                gTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAg        >  2:161153/1‑136 (MQ=255)                                                                                                gTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAg        >  2:327616/1‑136 (MQ=255)                                                                                                gTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCAATCTGCAGACGACGGCCTTCGAGGAAg        >  2:103679/1‑136 (MQ=255)                                                                                                gTGGACCTGACTGGGTTCGGACACCAGTCCCAGGTCGCCGTCGTCACTGCCGTAGAAATTCGCTGGTGGATAACCGCCCATGAAGCATTGCAGCAGCACGCCCCGGGCCATCTCCAGACGACGGCCTACGAGGAAg        >  2:36315/1‑136 (MQ=255)                                                                                                    aCCTGACTTGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAGCCtt    <  2:301436/136‑1 (MQ=255)                                                                                                    aCCTGACTTGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAGCCtt    <  2:340383/136‑1 (MQ=255)                                                                                                    aCCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAGCCtt    <  2:197098/136‑1 (MQ=255)                                                                                                    aCCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAGCCtt    <  2:244459/136‑1 (MQ=255)                                                                                                    aCCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAGCCtt    <  1:22974/136‑1 (MQ=255)                                                                                                    aCCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAGCCtt    <  1:207924/136‑1 (MQ=255)                                                                                                      cTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAGCCTTcg  >  1:287219/1‑136 (MQ=255)                                                                                                            |                                                                                                                                  AGCTGGCCCCGTGTGGACAGTGCCGAGGGCCGGCCATCCAGCTCGTTGGCATCACAGATGGCCAAGGCGCGGCCGTCATCGCTCCAGGGCAGGCGGTGGACCTGACTGGGTTCGGCCACCAGTGCCAGGTCGCCGTCGTCACTGCCGTAGAATTTCGCTGGTGGATAACCGCCCATGATGCATTGCAGCAGCACGCCCCGGGCCATCTGCAGACGACGGCCTTCGAGGAAGCCTTCG  >  NC_002947/4995884‑4996120

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