Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,688,515	C→A	14.4%	G28C (GGT→TGT)	bcsA ←	cellulose synthase, catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,688,515	0	C	A	14.4%	102.4 / 2.9	41	G28C (GGT→TGT)	bcsA	cellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (18/17);  new base A (3/3);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

ACCAGCATGTTTGAATGAGATAACGGACCGGGTCCAGCGGACGCGGACGCGAGGCGTTGATATGCGGATACAGGTTTTTATGTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCATTCAGCCATCGCTTCATCACGATGAA  >  CP009273/3688332‑3688750                                                                                                                                                                                        |                                                                                                                                                                                                                                            aCCAGCATGTTTGAATGAGATAACGGACCGGGTCCAGCGGACGCGGACGCGAGGCGTTGATATGCGGATACAGGTTTTTATGTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCTACGCAcc                                                                                                                                                                                                                                             <  2:204272/184‑1 (MQ=255)                    aTAACGGACCGGGTCCAGCGGACGCGGACGCGAGGCGTTGATATGCGGATACAGGTTTTTATGTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAacca                                                                                                                                                                              >  1:84545/1‑228 (MQ=255)                        cgtaCCGGGTCCAGCGGACGCGGACGCGAGGCGTTGATATGCGGATACAGGTTTTTATTTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAacc                                                                                                                                                                               >  2:250393/4‑223 (MQ=255)                        cGGACCGGGTCCAGCGGACGCGGACGCGAGGCGTTGATATGCGGATACAGGTTTTTATGTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCgg                                                                                                                                                                          <  2:84545/228‑1 (MQ=255)                            ccGGGTCCAGCGGACGCGGACGCGAGGCGTTGATATGCGGATACAGGTTTTTATTTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCg                                                                                                                                                                           <  1:250397/223‑1 (MQ=255)                                                                 cGGATACAGGTTTTTATGTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGcc                                                                                                                                                                                                                   >  2:75147/1‑146 (MQ=255)                                                                     tACAGGTTTTTATGTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGAt                                                                                                                                                                                                               <  1:75147/146‑1 (MQ=255)                                                                     tACAGGTTTTTANGTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGATCGTCGCGCTGAAAGCCGACGTACCGTGATGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGCCATACGCAACCACTGGGTCAGGATACTCa                                                                                                                                                             <  1:479574/196‑1 (MQ=255)                                                                        aGGTTTTTATGTTCTGCGCGAATACGCGGCAAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATAGCGATAACGCCCGATAAGc                                                                                                                                                                                                           <  1:532114/147‑1 (MQ=255)                                                                           tttttATGTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTc                                                                                                                                                                       >  2:572053/1‑180 (MQ=255)                                                                               tATGTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGa                                                                                                                                                                   <  1:572064/180‑1 (MQ=255)                                                                                  gTTCTGCGTGAATACGCTGCCAGAGCGGGTGCTCCAGAGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGt                                                                                                                                                                        <  1:217563/172‑1 (MQ=255)                                                                                                  cTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGTCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCaa                                                                                                                                                                                 >  2:601219/1‑147 (MQ=255)                                                                                                  cTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCg                                                                                                                                                                                                                                            >  1:155879/1‑88 (MQ=255)                                                                                                      cAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGTCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAaccac                                                                                                                                                                             <  1:601230/147‑1 (MQ=255)                                                                                                      cAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTgac                                                                                                                                                                                                                                        <  2:155876/88‑1 (MQ=255)                                                                                                             ggTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACCCCCGACTGGCGTATTCAGCCCGGGGTAGGTCAACAGGCACCAGTTCGCCAGCGTCAGGATCTCTTCAGCCGCCAGCGCATCACTGCGATATTAAACTACTGGTTGcagt                                       >  1:367345/1‑271 (MQ=255)                                                                                                                            tAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACAGTGACGACGATAATCGCGATa                                                                                                                                                                                                                        >  2:433086/1‑82 (MQ=255)                                                                                                                                aaTCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACAGTGACGACGATAATCGCGATAACGc                                                                                                                                                                                                                    <  1:433092/82‑1 (MQ=255)                                                                                                                                       gccaggATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACTCACCGTGACGACGATAATCGCGa                                                                                                                                                                                                                          >  1:302763/1‑69 (MQ=255)                                                                                                                                       gccaggAGCATCCAGAAACAGCCGAGCGTAGCGCTGAAAGCCGACGAACCGTGACGACGATAATCGCGAGAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGCCTGGCGTTTTCAGCTCTGAGTAGTTCAACAGTCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCAt                                                                    <  2:34360/219‑1 (MQ=255)                                                                                                                                                 tcCAGAACCATCCGAGCTTCGCGCTGAAAGCCCACGCACCGTGACGACCATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGTTAAGCAACCACCGGGTCAGGATACTCATTGTCCACTCCCGACTGGCTTGTTCCGCCCGGCACAGTTACACAGCCACCAGTTCTCAACCGCCCGTATCTCTCCAGCCCCCAGCGCATCCCAGCGCTACTCACCCCCTCCTCGCTGCTCCTCCAGGAATTCCGCCATCTCTacacccaactcaa  >  2:203292/1‑265 (MQ=255)                                                                                                                                                 tcCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGATGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACAACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCCACAGGCCCCTGTTCGc                                                                                                    >  1:450689/1‑177 (MQ=255)                                                                                                                                                 tcCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACcggcg                                                                                                                                                                                           >  1:435270/1‑90 (MQ=255)                                                                                                                                                 tcCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGAt                                                                                                                                                                                       >  2:570368/1‑94 (MQ=255)                                                                                                                                                 tcCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCaa                                                                                                                                                                                 >  2:241495/1‑100 (MQ=255)                                                                                                                                                 tcCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTc                                                                                                                                                                       >  2:622919/1‑110 (MQ=255)                                                                                                                                                 tcCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGCGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCATTCACCCATCGCTtccacaccacga   >  1:590655/1‑264 (MQ=255)                                                                                                                                                 tcCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACAGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTg                                                                                                                                                                                                  >  1:142247/1‑83 (MQ=255)                                                                                                                                                 tcCAGAAACAGCCGAGCGTCGCGCTGAAAGCAGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCaa                                                                                                                                                                                 >  2:652856/1‑100 (MQ=255)                                                                                                                                                     gAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGAt                                                                                                                                                                                       <  2:435264/90‑1 (MQ=255)                                                                                                                                                     gAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGc                                                                                                                                                                                   <  1:570379/94‑1 (MQ=255)                                                                                                                                                     gAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAaccac                                                                                                                                                                             <  1:241499/100‑1 (MQ=255)                                                                                                                                                     gAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGa                                                                                                                                                                   <  1:622930/110‑1 (MQ=255)                                                                                                                                                     gAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACAGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACcg                                                                                                                                                                                              <  2:142246/83‑1 (MQ=255)                                                                                                                                                     gAAACAGCCGAGCGTCGCGCTGAAAGCAGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAaccac                                                                                                                                                                             <  1:652867/100‑1 (MQ=255)                                                                                                                                                            ccGAGCGTCGCGCTGAAAGCCGACGCACAGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCa                                                                                                                                                                                  >  1:608269/1‑88 (MQ=255)                                                                                                                                                                gCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCAGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAgccgcc                                 <  1:357566/229‑1 (MQ=255)                                                                                                                                                                gCGTCGCGCTGAAAGCCGACGCACAGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAacca                                                                                                                                                                              <  2:608258/88‑1 (MQ=255)                                                                                                                                                                      cgcTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAaccac                                                                                                                                                                             >  1:95248/1‑83 (MQ=255)                                                                                                                                                                          gAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCggg                                                                                                                                                                         <  2:95248/83‑1 (MQ=255)                                                                                                                                                                                   cGCACCGTGACGACGATAATCTCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGatat                                                         >  2:608739/1‑186 (MQ=255)                                                                                                                                                                                       ccGTGACGACGATAATCTCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCAc                                                     <  1:608750/186‑1 (MQ=255)                                                                                                                                                                                        |                                                                                                                                                                                                                                            ACCAGCATGTTTGAATGAGATAACGGACCGGGTCCAGCGGACGCGGACGCGAGGCGTTGATATGCGGATACAGGTTTTTATGTTCTGCGCGAATACGCTGCCAGCGCGGGTGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCATTCAGCCATCGCTTCATCACGATGAA  >  CP009273/3688332‑3688750

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