Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,141,493	Δ1 bp	10.0%	coding (410/972 nt)	PP_2757 →	sugar‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,141,493	0	T	.	10.0%	139.5 / 11.0	50	coding (410/972 nt)	PP_2757	sugar‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (25/18);  new base . (2/3);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

GCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGAT  >  NC_002947/3141347‑3141629                                                                                                                                                   |                                                                                                                                          gCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTcg                                                                                                                                         >  2:350269/1‑149 (MQ=255)  cAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTcgc                                                                                                                                        <  1:55528/149‑1 (MQ=255)   aaaaGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGTGCGTTATCGACCCCTGGCTCAAGCGCATCCAGGCAGGAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTcgcg                                                                                                                                       <  2:365866/149‑1 (MQ=255)    aaaGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGcgg                                                                                                                                      <  4:135407/149‑1 (MQ=255)                ggTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGt                                                                                                                         >  5:300045/1‑149 (MQ=255)                ggTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGt                                                                                                                         >  8:156541/1‑149 (MQ=255)                 gTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGtt                                                                                                                        >  5:247816/1‑149 (MQ=255)                  tAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTg                                                                                                                       <  7:156541/149‑1 (MQ=255)                    aTTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGCCAACG‑CGCGGCCGGCCGCGACGTTgtc                                                                                                                     <  6:163284/149‑3 (MQ=255)                    aTTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGc                                                                                                                     <  6:300045/149‑1 (MQ=255)                        aGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGCCAACG‑CGCGGCCGGCCGCGACGTTGTcc                                                                                                                    >  4:17677/1‑146 (MQ=255)                        aGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGCCAACG‑CGCGGCCGGCCGCGACGTTGTcc                                                                                                                    <  3:17677/146‑1 (MQ=255)                            ccTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGc                                                                                                             <  6:306216/149‑1 (MQ=255)                            ccTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGc                                                                                                             >  8:166188/1‑149 (MQ=255)                              tCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTg                                                                                                           >  5:214129/1‑149 (MQ=255)                                    cccTGAGCGGGATCGCCCCCTGGCTCAAGCGCAGCCAGGCACCAGGTATCCGTCTTTTCACCATTGATGTGCCTTCGAAGTACAGCCTCAACACCACCACCTCGGACAACGCCGCGGCCGGCCGCG‑GGTTGGCCGAACAGCTGGTgcgc                                                                                                     <  6:247816/149‑1 (MQ=255)                                                 cGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCggcggc                                                                                         >  8:194525/1‑148 (MQ=255)                                                    ccccTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGCCAACG‑CGCGGCCGGCCGCGACGTTgtc                                                                                                                     >  3:90169/1‑115 (MQ=38)                                                    ccccTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGCCAACG‑CGCGGCCGGCCGCGACGTTgtc                                                                                                                     <  4:90169/117‑3 (MQ=38)                                                            tCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATc                                                                             >  3:14534/1‑149 (MQ=255)                                                                 cgcATCCAGGCGGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGgt                                                                        <  4:50070/149‑1 (MQ=255)                                                                       cAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGTCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCaa                                                                  >  4:209653/1‑149 (MQ=255)                                                                       cAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCaa                                                                  <  3:38474/149‑1 (MQ=255)                                                                        aGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACcggcgg                                                                                          >  6:222259/1‑124 (MQ=255)                                                                        aGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACcggcgg                                                                                          <  5:222259/124‑1 (MQ=255)                                                                                ggTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCaa                                                                                       >  5:144684/1‑119 (MQ=255)                                                                                ggTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCaa                                                                                       <  6:144684/119‑1 (MQ=255)                                                                                     cccGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACgggg                                                    >  1:320319/1‑149 (MQ=255)                                                                                                 cATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGAGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCtgtg                                        >  5:148451/1‑149 (MQ=255)                                                                                                  aTTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGgccggccg                                                                                                                               >  2:49850/1‑62 (MQ=255)                                                                                                  aTTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGgccggccg                                                                                                                               <  1:49850/62‑1 (MQ=255)                                                                                                  aTTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGc                                       >  5:273039/1‑149 (MQ=255)                                                                                                      aggtgCGATCGCAGTACAGCCTCAACAACACCACCTCGGACTACTTCGCGGCCGGCCGCG‑CGTTGTCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCTCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATc                                   <  6:107432/147‑1 (MQ=255)                                                                                                         tgCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGc                                >  2:297466/1‑149 (MQ=255)                                                                                                           ccATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTa                              <  6:273039/149‑1 (MQ=255)                                                                                                                gCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTacgac                          <  7:166188/148‑1 (MQ=255)                                                                                                                   gTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGt                      <  7:194525/149‑1 (MQ=255)                                                                                                                      cAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTga                   <  1:350269/149‑1 (MQ=255)                                                                                                                         ccTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTgaag                 >  5:159477/1‑148 (MQ=255)                                                                                                                           tCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGGCTGTGCGATCCGCTACGACGAGTTGaagaag              >  7:258943/1‑149 (MQ=255)                                                                                                                            caacaaCACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAgg             >  8:284296/1‑149 (MQ=255)                                                                                                                                 acacCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTg        >  8:314394/1‑149 (MQ=255)                                                                                                                                    ccaccTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGcc     >  7:252713/1‑149 (MQ=255)                                                                                                                                     caccTCGGACAACGTCGCGGCCGGTCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCg    <  3:209653/149‑1 (MQ=255)                                                                                                                                     caccTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGcc     >  8:268243/1‑148 (MQ=255)                                                                                                                                        cTCGGACAACGTCTCGGCCTGCCGCG‑CGTTGGCAGAACAGCTGTTGCGCTACACCGGCGGCAAGGGCAGCATCCTGTTGTTCAACGGCTTCTACTGGGTGACGGTCTGTTCGATCAGCTACGACGAGTTGAAGAAGGTACTGGCCGat  >  4:223799/1‑148 (MQ=255)                                                                                                                                             aCAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCtgt                                         >  1:301900/1‑104 (MQ=255)                                                                                                                                             aCAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCtgt                                         <  2:301900/104‑1 (MQ=255)                                                                                                                                               aaCGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTAc          >  4:362848/1‑133 (MQ=255)                                                                                                                                               aaCGGCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTAc          <  3:362848/133‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                          GCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCG‑CGTTGGCCGAACAGCTGGTGCGCGACACCGGCGGCAAGGGCCGCATCCTGGTGTTCAACGGCTTCTACGGGGTGCCGGTCTGTGCGATCCGCTACGACGAGTTGAAGAAGGTACTGGCCGAT  >  NC_002947/3141347‑3141629

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