Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,010,969	A→C	6.9%	V10V (GTA→GTC)	bcsG →	endoglucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,010,969	0	A	C	6.9%	167.9 / 2.5	55	V10V (GTA→GTC)	bcsG	endoglucanase
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base C (2/2);  total (29/26)
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.86e-01

CTTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGC  >  NC_002947/3010824‑3011117                                                                                                                                                  |                                                                                                                                                     ctTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTg                                                                                                                                                   >  7:329666/1‑149 (MQ=255) ctTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTg                                                                                                                                                   >  6:184985/1‑149 (MQ=255)    gttgCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGcc                                                                                                                                                >  7:368/1‑149 (MQ=255)       gCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGc                                                                                                                                             >  4:202173/1‑149 (MQ=255)              cgcCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTg                                                                                                                                      <  8:329666/149‑1 (MQ=255)               gcCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACAGCTGGGCGTACGCCGCCGCGCTGCACGCCACGAACCCCAGAATGACTAATACAACCATCACCATCCCGGTCCTGGCCCGCTGGGCTg                                                                                                                                      >  3:3150/1‑148 (MQ=255)                 cTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGAc                                                                                                                                   >  5:197889/1‑149 (MQ=255)                 cTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGccagtac                                                                                                                                   <  4:285774/149‑6 (MQ=255)                    tgtTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCAGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                <  2:149957/149‑1 (MQ=255)                         ggCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAggg                                                                                                                           >  2:200777/1‑149 (MQ=255)                         ggCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCAGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTTGGAggg                                                                                                                           >  4:279401/1‑149 (MQ=255)                         ggCATTACGTGGCGTGTTGCGCCGTCGACTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAgg                                                                                                                            <  6:365997/148‑1 (MQ=38)                          gCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTCCCTGAAACCGGTGGGGGTACGCCGCCGCGCTGCACTCCAGGAACCCCAGTATGACTAATACCACCATCACCGTCCCGGTACTGGCCCGCTGGGCTGGACTGGGAGGGt                                                                                                                          <  3:138652/149‑1 (MQ=255)                           cATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTg                                                                                                                         >  7:350855/1‑149 (MQ=255)                             ttACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGa                                                                                                                       <  2:43080/149‑1 (MQ=255)                                  tGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTCCGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGGTGACCctc                                                                                                                  >  4:193477/1‑149 (MQ=255)                                  tGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACctc                                                                                                                  <  5:185377/149‑1 (MQ=255)                                  tGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACctc                                                                                                                  >  6:185377/1‑149 (MQ=255)                                    gCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTa                                                                                                                >  7:326802/1‑149 (MQ=255)                                        gtTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCTCGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAtttt                                                                                                            >  4:256408/1‑149 (MQ=255)                                             gcCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTg                                                                                                        >  5:123682/1‑148 (MQ=255)                                                   gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCTCCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGtt                                                                                                 >  3:130032/1‑149 (MQ=255)                                                     cTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTc                                                                                               <  8:326802/149‑1 (MQ=255)                                                     cTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTc                                                                                               <  6:197889/149‑1 (MQ=255)                                                         ccccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTctgc                                                                                            <  6:123682/148‑1 (MQ=255)                                                           ccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTgg                                                                                         >  3:36/1‑149 (MQ=255)                                                              cGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTcg                                                                                      >  4:242732/1‑149 (MQ=255)                                                                    cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                                <  8:368/149‑1 (MQ=255)                                                                      gAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                               >  6:226107/1‑148 (MQ=255)                                                                       aaaCCGCTGGTCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGTAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                               <  4:3150/147‑1 (MQ=255)                                                                            gCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCTCCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTc                                                                        <  4:130032/149‑1 (MQ=255)                                                                                ggCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACt                                                                    <  5:184985/149‑1 (MQ=255)                                                                                      cgccgccgcGCTGCTCGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAg                                                               <  3:256408/148‑1 (MQ=255)                                                                                           ccgcGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTg                                                         >  5:345557/1‑149 (MQ=255)                                                                                                        cAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCtgttgt                                            >  6:81176/1‑149 (MQ=255)                                                                                                          ggAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTc                                          >  8:140133/1‑149 (MQ=255)                                                                                                             aCCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCt                                       <  3:202173/149‑1 (MQ=255)                                                                                                               cccAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGc                                     <  6:345557/149‑1 (MQ=255)                                                                                                                 cAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGt                                   >  5:275328/1‑149 (MQ=255)                                                                                                                        aCTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTggg                                                                                                                               <  5:25486/50‑1 (MQ=255)                                                                                                                        aCTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTggg                                                                                                                               >  6:25486/1‑50 (MQ=255)                                                                                                                        aCTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                             <  5:81176/148‑1 (MQ=255)                                                                                                                                caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTcgc                                                                                     <  5:233884/84‑1 (MQ=255)                                                                                                                                caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTcgc                                                                                     >  6:233884/1‑84 (MQ=255)                                                                                                                                caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCtt                    >  4:218445/1‑149 (MQ=255)                                                                                                                                caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCtt                    >  5:329099/1‑149 (MQ=255)                                                                                                                                    aTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCttt                   <  7:88985/146‑1 (MQ=255)                                                                                                                                    aTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCttt                   >  8:88985/1‑146 (MQ=255)                                                                                                                                     tCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGGCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCg               <  3:279401/149‑1 (MQ=255)                                                                                                                                       aCCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCggg             <  6:275328/149‑1 (MQ=255)                                                                                                                                         ccTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGc            >  2:191968/1‑148 (MQ=38)                                                                                                                                          cTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCaa          <  8:350855/149‑1 (MQ=255)                                                                                                                                             ccGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCaatac       >  4:236585/1‑146 (MQ=255)                                                                                                                                                gTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGTCaatacgca    <  3:193477/149‑7 (MQ=255)                                                                                                                                                gTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCt    <  7:140133/149‑1 (MQ=255)                                                                                                                                                  aCTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGc  >  5:370008/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CTTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGC  >  NC_002947/3010824‑3011117

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