Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,851,260	C→A	A115S (GCG→TCG)	acsA1 ←	acetate‑‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,851,260	0	C	A	96.5%	90.1 / ‑2.2	29	coding (343/1707 nt)	acsA1	acetate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (12/16);  minor base . (1/0);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CTACAGGTACTCGAGGTAACAATGCCTTATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGCAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGTATACGATTCTTTTCTGTGC  >  CP000730/1851116‑1851399                                                                                                                                                 |                                                                                                                                            cTACAGGTACTCGAGGTAACAATGCCTTATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGaaaa                                                                                                                                          <  2:139094/148‑1 (MQ=255) cTACAGGTAATCGAGGTAACAATGCCTTATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGaaaa                                                                                                                                          <  2:82892/148‑1 (MQ=255)     aGGTACTCGAGGTAACAATGCCTTATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAatat                                                                                                                                      >  1:23854/1‑148 (MQ=255)                            tattaGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACataaa                                                                                                               <  1:52965/148‑1 (MQ=255)                            tattaGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACataaa                                                                                                               <  1:231332/148‑1 (MQ=255)                             attaGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACataaat                                                                                                              <  1:387005/148‑1 (MQ=255)                             attaGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACataaat                                                                                                              <  1:184700/148‑1 (MQ=255)                                aGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACataaatata                                                                                                           >  1:25265/1‑148 (MQ=255)                                                            ttataaagatatgCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTg                                                                               >  2:196869/13‑148 (MQ=255)                                                              tAATCTATCCGCAACTGCCTTTTCCATAAATGTTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCa                                                                             >  2:398463/1‑148 (MQ=255)                                                              tAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCa                                                                             >  2:414071/1‑148 (MQ=255)                                                              tAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCa                                                                             >  2:336839/1‑148 (MQ=255)                                                                 tCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGt                                                                          >  1:402702/1‑148 (MQ=255)                                                                  cTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACG‑AAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTc                                                                        >  1:277982/1‑148 (MQ=255)                                                                                cTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTc                                                           >  1:353852/1‑148 (MQ=255)                                                                                     ccATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGTCATAATTATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCgctgc                                                      <  1:95628/148‑1 (MQ=255)                                                                                     ccATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCgctgc                                                      <  1:58300/148‑1 (MQ=255)                                                                                         aaaTGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTtta                                                  >  1:281589/1‑148 (MQ=255)                                                                                         aaaTGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTtta                                                  >  1:340304/1‑148 (MQ=255)                                                                                            tGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTAAGGTGTACTCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATTTTCAGACAAAACATTCGCTGCTttattt                                               >  2:290412/1‑147 (MQ=255)                                                                                                     tAACGTCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGt                                      <  2:162127/148‑1 (MQ=255)                                                                                                     tAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCTACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGt                                      >  2:314261/1‑148 (MQ=255)                                                                                                      aaCGGCCCAACAATTTCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGtt                                     <  2:165937/148‑1 (MQ=255)                                                                                                                    tGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATaa                       <  1:196869/148‑1 (MQ=255)                                                                                                                    tGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATaa                       <  1:414021/148‑1 (MQ=255)                                                                                                                       aCCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGta                    <  1:48167/148‑1 (MQ=255)                                                                                                                       aCCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGta                    <  1:43054/148‑1 (MQ=255)                                                                                                                                   aaCACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGTATACGATTCtttt        <  2:297631/148‑1 (MQ=255)                                                                                                                                   aaCACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGTATACGATTCtttt        <  1:30078/148‑1 (MQ=255)                                                                                                                                     cacCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGTATACGATTCTTTTCt      <  1:130551/148‑1 (MQ=255)                                                                                                                                         tAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGTATACGATTCTTTTCTGTGc  <  2:402702/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                            CTACAGGTACTCGAGGTAACAATGCCTTATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGCAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGTATACGATTCTTTTCTGTGC  >  CP000730/1851116‑1851399

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