Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	336,125	+T	coding (210/237 nt)	PP_0278 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	336,124	1	.	T	97.6%	136.4 / ‑3.0	42	coding (209/237 nt)	PP_0278	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (22/19);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

CCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑T‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCA  >  NC_002947/335987‑336261                                                                                                                                           |                                                                                                                                           ccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGc                                                                                                                                  <  1:124937/148‑1 (MQ=255)      aTCTCCTACGCGTAGAGCCCTGCGTCCCCGTCGACGCGGCCTACGAGAATGTCTTGATCCTGATTGGCTACATCAAGCACCTGGTGCGCGATGGCGATATGGAGGACTACCATAAATTCCTGGTGTCAGCCGA‑TCACCTGTGCGACATg                                                                                                                            <  7:298494/149‑1 (MQ=255)           cTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGcc                                                                                                                         >  7:167008/1‑147 (MQ=255)             aCGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGa                                                                                                                                             <  5:50079/126‑1 (MQ=255)             aCGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGa                                                                                                                                             >  6:50079/1‑126 (MQ=255)              cgcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGgc                                                                                                                    <  8:193813/149‑1 (MQ=255)                 gTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCt                                                                                                                 >  1:214644/1‑149 (MQ=255)                  tAGAGCCCGGCGTCCCGGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCtg                                                                                                                <  4:340473/149‑1 (MQ=255)                   agagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCtga                                                                                                               >  8:202906/1‑149 (MQ=255)                    gagCACGGCGTCCCCGTCGACGAGGCCTACGAGCATGTCTCGATCCTTATGGGCTACATCAAGCACCTTGTGCGCGAAGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCtga                                                                                                               <  5:240335/148‑1 (MQ=255)                        ccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTAAATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGTGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAAc                                                                                                          <  5:285987/149‑1 (MQ=255)                           gCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGat                                                                                                       <  7:202906/149‑1 (MQ=255)                           gCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGa                                                                                                        >  3:149152/1‑148 (MQ=255)                               ccccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGa                                                                                                                                             <  6:219546/108‑1 (MQ=255)                               ccccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGa                                                                                                                                             >  5:219546/1‑108 (MQ=255)                                 ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                    >  1:165135/1‑146 (MQ=255)                                 ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                    >  5:27943/1‑146 (MQ=255)                                 ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCgag                                                                                                 <  5:72941/149‑1 (MQ=255)                                       aCGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGc                                                                                            >  8:258940/1‑148 (MQ=255)                                              cTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGAc                                                                                                                               >  5:205719/1‑106 (MQ=255)                                              cTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGAc                                                                                                                               <  6:205719/106‑1 (MQ=255)                                                 cGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAgg                                                                                 <  7:258940/149‑1 (MQ=255)                                                         tctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAg                                                                         <  4:149152/149‑1 (MQ=255)                                                          ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                        >  8:18912/1‑149 (MQ=255)                                                               tCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAtt                                                                    >  1:978/1‑148 (MQ=255)                                                                        gCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCa                                                          >  6:101764/1‑149 (MQ=255)                                                                         cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAt                                                         >  7:193399/1‑149 (MQ=255)                                                                         cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAt                                                         <  6:27943/149‑1 (MQ=255)                                                                            cATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAttt                                                       <  5:101764/148‑1 (MQ=255)                                                                               cAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                   >  2:340543/1‑149 (MQ=255)                                                                                       tGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGcc                                           >  6:135186/1‑149 (MQ=255)                                                                                        ggTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCt                                          >  2:90271/1‑149 (MQ=255)                                                                                            cgcgAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGc                                       <  8:167008/148‑1 (MQ=255)                                                                                                    cGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAAt                              <  1:90271/149‑1 (MQ=255)                                                                                                     gATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGACGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGaaag                             >  6:38825/1‑147 (MQ=255)                                                                                                       tatGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                           >  6:100507/1‑149 (MQ=255)                                                                                                           gaggaCGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTg                        >  6:271087/1‑148 (MQ=255)                                                                                                             ggaCGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCaa                     <  7:18912/149‑1 (MQ=255)                                                                                                               acgacCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAag                   <  5:271087/149‑1 (MQ=255)                                                                                                                      tAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCTTTGCAGccca            <  5:135186/149‑2 (MQ=255)                                                                                                                       aaaTTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCg           >  5:31960/1‑149 (MQ=255)                                                                                                                         aTTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGcccccgt         >  6:42899/1‑146 (MQ=38)                                                                                                                          ttCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTg                        >  3:63288/1‑133 (MQ=255)                                                                                                                          ttCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTg                        <  4:63288/133‑1 (MQ=255)                                                                                                                          ttCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCg        <  6:31960/149‑1 (MQ=255)                                                                                                                                   gCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCa  >  1:91489/1‑146 (MQ=255)                                                                                                                                           |                                                                                                                                           CCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑T‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCA  >  NC_002947/335987‑336261

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