Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,679,656	C→T	11.1%	G960S (GGC→AGC)	bcsC ←	cellulose synthase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,679,656	0	C	T	11.1%	97.4 / 2.5	36	G960S (GGC→AGC)	bcsC	cellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base C (16/16);  new base T (2/2);  total (18/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.73e-01

GCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGGCCGATTGTGACGCGGCGATTGTTCTGGTTAATGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCC  >  CP009273/3679408‑3679895                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                aagcgacagCCCCATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCTACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGGCCGATTGTGACGCGGCGATTGTTCTGGTTAATGACCTTCTAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTc                                                                                                                                                                                                                        <  1:249856/266‑1 (MQ=255)        cGCCACATCTCCTGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGNCCGATTGTGCTACGGCGATTGTTCTGGTTAATGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACt                                                                                                                                                                                                                   <  1:603688/272‑1 (MQ=255)                     gTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGGGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGGCCGAATGTGGCGCGGCGATTGTTCTGGTTAATGACCGTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTGGTCTTCGACGTTTTTACAGGTTAACGGGTCGCCACTAAGCGATGCCCAGACGCCGTTGGCGTCTCCTTTCGCGTACCTCGGACTTAGcccncgcctctag                                                                                                                                                                                                   >  2:249852/1‑265 (MQ=255)                                                               ccTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGGCCGATTGTGACGCGGCGATTGTTCTGGTTAATGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCACGTTATGGCACGGGGAg                                                                                                                                                         >  1:398471/1‑275 (MQ=255)                                                               ccTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCAAGATCATGTTGTTCAGGCCGATTGTGACGCGGCGATTGTTCTGGTTAATGACCTTATAGTAATAGCCCGTCATCCAGCGCCAGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCAACTAAGCGATGCCCAGACCTCGTTTGCTTCACCTTTATCGTCGCTCAGACTTATCACCCAACCGCCGGCACGTAAGCCCCCCCATTTTTTCCCGGTAtgtctcccgtc                                                                                                                                                          >  2:168606/1‑264 (MQ=255)                                                                   aCCGAGTTAGTAGACACTCAGATCTTGGTCGTAATGCCAGATCATGTTGTTCAGGCCGATTGTGACGCGGCGATTGTTCTGGTTAAGGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTcc                                                                                                                                                                          <  1:168609/254‑1 (MQ=255)                                                                           agtAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGGCCGATTGTGACGCGGCGATTGTTCTGGTTAATGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGt                                                                                                                                                                                                                                               >  1:72666/1‑177 (MQ=255)                                                                               gCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGGCCGATTGTGACGCGGCGATTGTTCTGGTTAATGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGc                                                                                                                                                                                                                                           <  2:72666/177‑1 (MQ=255)                                                                                        ctCTTTGTCGAAACNCCAGATCATTTTGTTCAGGCCGATTGTGACGCGGCGATTGTTCTGGTTAATGACCTTATAGTAATAGCCCGTCAGCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCa                                                                                                                                <  2:493022/274‑1 (MQ=255)                                                                                                                                  gACGCGGCGATTGTTCTGGTTAATGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCaaa                                                                                                                              <  1:200311/235‑1 (MQ=255)                                                                                                                                                 cTGGTTAATGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCtt                                                                                                                                                                                                                                   >  2:534482/1‑119 (MQ=255)                                                                                                                                                     ttAATGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATc                                                                                                                                                                                                                               <  1:534493/119‑1 (MQ=255)                                                                                                                                                         tGACCTTCTAGTAATAGCCCGTCATCCAGCGCACGCGCCACTTATCTTCGACCTTTTTACCGGTTAACTGGTCGCCACTACGCGATGCCCAACCGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTCATCCCCACCCCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCACCAAATAACTGGAGATGGGCCGACGGTGGCCGTTAACGGCGTAACCACGCGGCCCGCTNTCATCTCTg                                                               >  2:249734/1‑275 (MQ=255)                                                                                                                                                         tGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCCACGGCCAGCAAAGCACTTGACATGGGCCGACGGTGGGCGGTCACGGTGGCACCCAGCGGccccacagcatgccac                                                               >  2:41115/1‑262 (MQ=255)                                                                                                                                                                ataGTAGTAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGAGGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGa                                                        <  1:249738/275‑1 (MQ=255)                                                                                                                                                                      aTAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCtttt                                                                                                                                                   >  2:543719/1‑178 (MQ=255)                                                                                                                                                                      aTAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTCTTTGCCCACCAAAGGCCAGCAAAGACCTGGAAATGGGCCGCCGGTCGGCGTTACCGGTGTACCCCAACGCCCCGCTATCACCCCGGTAACTGATtccc                                                   >  2:624177/1‑273 (MQ=255)                                                                                                                                                                       tAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCg                                                                                                                                                                                                                              >  2:584165/1‑102 (MQ=255)                                                                                                                                                                          cccGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGt                                                                                                                                                                                                                             >  1:6204/1‑100 (MQ=255)                                                                                                                                                                          cccGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGcc                                                                                                                                               <  1:543730/178‑1 (MQ=255)                                                                                                                                                                           ccGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGc                                                                                                                                                                                                                          <  1:584176/102‑1 (MQ=255)                                                                                                                                                                            cGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGt                                                                                                                                                                                                                             >  1:541440/1‑98 (MQ=255)                                                                                                                                                                              tCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCg                                                                                                                                                                                                                              >  2:547322/1‑95 (MQ=255)                                                                                                                                                                              tCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCt                                                                                                                                                                                                                         <  2:6204/100‑1 (MQ=255)                                                                                                                                                                                aTCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCt                                                                                                                                                                                                                         <  2:541429/98‑1 (MQ=255)                                                                                                                                                                                  ccAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGc                                                                                                                                                                                                                          <  1:547333/95‑1 (MQ=255)                                                                                                                                                                                     gcgcACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCTGTTTGCTTCACCTTTATCGTAGCTCAGAc                                                                                                                                                                                                                    >  2:585648/1‑98 (MQ=255)                                                                                                                                                                                         aCGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCTGTTTGCTTCACCTTTATCGTAGCTCAGACTTAg                                                                                                                                                                                                                <  1:585659/98‑1 (MQ=255)                                                                                                                                                                                                                      tgTTACCTGGTCGCCACTAAGCGATTCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATccc  <  1:41115/274‑1 (MQ=255)                                                                                                                                                                                                                        ttAACTGGTCGCCACTAAGCGATGCCCAGACGCTGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCTGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGata                                                                        >  2:165871/1‑203 (MQ=255)                                                                                                                                                                                                                         tAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCtt                                                                                                                                                                                                                                   >  2:36648/1‑47 (MQ=255)                                                                                                                                                                                                                            cTGGTCGCCACTAAGCGATGCCCAGACGCTGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCTGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATatcat                                                                    <  1:165874/203‑1 (MQ=255)                                                                                                                                                                                                                             tGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATc                                                                                                                                                                                                                               <  1:36648/47‑1 (MQ=255)                                                                                                                                                                                                                                 cGCCACTAAGCGATGCCCAGACGACGTATGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCCCACCGTCGGCACGTACGACACCCCATTTTTTCCCGGGATTGCTCGGGGCGTCTTTTTGCCCACCAAAGGCCAGCAAAGACCTGGAGATGGGCCGACGGAGGGCGTTAACGGTGTAACACAGCGGCCCGCTATCATCGCTGTAACTGATgccg                                                   >  1:403176/1‑215 (MQ=255)                                                                                                                                                                                                                                  gCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCt                                                                                                                                                                                                                         >  1:260677/1‑48 (MQ=255)                                                                                                                                                                                                                                      cTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGa                                                                                                                                                                                                                     <  2:260673/48‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                GCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGGCCGATTGTGACGCGGCGATTGTTCTGGTTAATGACCTTATAGTAATAGCCCGTCATCCAGCGCACGCGCCAGTTATCTTCGACATTTTTACCGGTTAACTGGTCGCCACTAAGCGATGCCCAGACGCCGTTTGCTTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCC  >  CP009273/3679408‑3679895

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