Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	4,276,215	C→A	R362L (CGC→CTC)	acs ←	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,276,215	0	C	A	96.9%	209.6 / ‑2.3	65	R362L (CGC→CTC)	acs	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (16/47);  total (16/49)
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.07e-01

CGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGCGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACG  >  CP009273/4276150‑4276279                                                                  |                                                                 cGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAgg                                                                >  1:3502818/1‑68 (MQ=255)  gagaATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGa                                                               <  1:4080739/68‑1 (MQ=255)       tGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGTGCGCGGATCGcc                                                          <  1:1185155/68‑1 (MQ=255)       tGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGcc                                                          <  1:4894044/68‑1 (MQ=255)          gcAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTg                                                       >  1:2023500/1‑68 (MQ=255)          gcAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTg                                                       >  1:5727492/1‑68 (MQ=255)           cAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTgg                                                      <  1:487856/68‑1 (MQ=255)             gCGACGAACTGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGGGAGGATCGCCgtgtgt                                                    <  1:3772317/68‑1 (MQ=255)             gCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCATCCGCCATCAGCGCGAGGATCGCCGTGGGt                                                    <  1:5431744/68‑1 (MQ=255)             gCGACGAACGGTCGGAGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGt                                                    <  1:3172291/68‑1 (MQ=255)               gacgaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCCGCGCGAGGATCGCCGTGGGTGc                                                  <  1:2227650/68‑1 (MQ=255)               gacgaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGc                                                  >  1:1776348/1‑68 (MQ=255)                 cgaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCgg                                                <  1:5180173/68‑1 (MQ=255)                  gatcggtcggtGCCTTCAATCGGTTTTTCGCCTTCCGCCATCAGCGCGAGGTTCGCCGTGGGTGCGGt                                               <  1:4664240/65‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCCGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:4598746/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:3560111/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:390387/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:497135/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:2933466/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:449556/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:5191950/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:1113569/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:1546675/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:2177897/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:5465122/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:209419/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:1966859/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:5467463/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGt                                               <  1:5921204/68‑1 (MQ=255)                  gaACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGAGGGTGCGGt                                               <  1:507609/68‑1 (MQ=255)                   aaCGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGta                                              <  1:3401382/68‑1 (MQ=255)                    aCGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGCGGATCGCCGTGGGTGCGGtat                                             <  1:4848122/68‑1 (MQ=255)                    aCGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGtat                                             <  1:2464556/68‑1 (MQ=255)                    aCGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGtat                                             <  1:5844594/68‑1 (MQ=255)                    aCGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGtat                                             <  1:2216229/68‑1 (MQ=255)                      ggtcggtGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATag                                           <  1:1208927/68‑1 (MQ=255)                      ggtcggtGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATag                                           <  1:2898180/68‑1 (MQ=255)                          ggtGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGaa                                       <  1:5695532/68‑1 (MQ=255)                          ggtGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGaa                                       <  1:1120617/68‑1 (MQ=255)                          ggtGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGaa                                       <  1:4558645/68‑1 (MQ=255)                           gtGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCTCGAGGATCGCCGTGGGTGCGGTATAGAGAat                                      <  1:4781641/68‑1 (MQ=255)                           gtGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAat                                      >  1:309376/1‑68 (MQ=255)                           gtGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAat                                      <  1:163008/68‑1 (MQ=255)                           gtGCCTTCGATCGCTTTATCGCCTGCCGCCATCCGCGCGATGATCGCCGTGGGTGCGGTATAGAGAat                                      <  1:1715684/68‑1 (MQ=255)                           gtGCCTTCGATCGCTTTAACGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAat                                      <  1:5882962/68‑1 (MQ=255)                             gCCGTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAatat                                    <  1:3247253/68‑1 (MQ=255)                              ccTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATAtt                                   >  1:2131424/1‑68 (MQ=255)                                 tCGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGAc                                <  1:2718996/68‑1 (MQ=255)                                  cGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGAcc                               <  1:4228559/68‑1 (MQ=255)                                  cGATCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGAcc                               <  1:1891515/68‑1 (MQ=255)                                    aTCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGCCCTg                             >  1:5924684/1‑68 (MQ=255)                                     tCGCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGa                            <  1:6075630/68‑1 (MQ=255)                                       gCTTTATCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATg                          >  1:5477815/1‑68 (MQ=255)                                       gCTTTATCGCCTTCCGACATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATg                          >  1:3047240/1‑68 (MQ=255)                                            aTCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGt                     <  1:5064755/68‑1 (MQ=255)                                             tCGCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTc                    >  1:2740484/1‑68 (MQ=255)                                               gCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTcca                  >  1:185474/1‑68 (MQ=255)                                               gCCTTCCGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTcaa                  >  1:4577492/1‑66 (MQ=255)                                                     cGCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCt            <  1:655893/68‑1 (MQ=255)                                                      gCCATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTg           >  1:2852767/1‑68 (MQ=255)                                                       ccATCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTgc          >  1:286002/1‑68 (MQ=255)                                                          tCAGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGcc       >  1:6206020/1‑68 (MQ=255)                                                          tCAGCGCGAGGATCGCCGTGGGAGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGcc       <  1:1881138/68‑1 (MQ=255)                                                           cAGCGCGAGGATCGCCGGGGGTGCGGTATAGAGAATATTGACCTGATGCTTGACCACCACCTGCGCCa      <  1:5065616/68‑1 (MQ=255)                                                            aGCGCGAGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCAt     >  1:3460794/1‑68 (MQ=255)                                                               gcgagGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACg  <  1:3843155/68‑1 (MQ=255)                                                                  |                                                                 CGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGCGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACG  >  CP009273/4276150‑4276279

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