Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,243,755	T→C	7.5%	intergenic (+44/+55)	tsaA → / ← rnt	peroxiredoxin/ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,243,755	0	T	C	7.5%	147.7 / 3.0	53	intergenic (+44/+55)	tsaA/rnt	peroxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base T (26/23);  new base C (2/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGTGTCA  >  NC_002947/1243611‑1243894                                                                                                                                                 |                                                                                                                                               gCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGc                                                                                                                                            >  8:197194/1‑149 (MQ=255)                 cAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGcc                                                                                                                           >  3:298294/1‑148 (MQ=255)                       gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGcc                                                                                                                               <  8:35924/138‑1 (MQ=255)                         aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                              <  6:264822/123‑1 (MQ=255)                         aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                              >  5:264822/1‑123 (MQ=255)                         aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCgg                                                                                                                  <  3:17802/149‑1 (MQ=255)                                 gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGccc                                                                                                          >  4:285251/1‑149 (MQ=255)                                       ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAAt                                                                                                    >  1:158431/1‑149 (MQ=255)                                        gCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTATTCGCGGGCGCGCCCGCGaa                                                                                                     >  6:255544/1‑147 (MQ=255)                                         cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTATTCGCGGGCGCGCCCGCGaagag                                                                                                  <  5:255544/149‑4 (MQ=255)                                           tGAACGGTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTTTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCt                                                                                                <  7:299941/149‑1 (MQ=255)                                            gAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                               >  2:50103/1‑149 (MQ=255)                                             aaCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                               <  7:304054/148‑1 (MQ=255)                                              aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                               >  8:62287/1‑147 (MQ=255)                                               cGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                            <  4:298294/149‑1 (MQ=255)                                                gCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                            >  1:276112/1‑148 (MQ=255)                                                  ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTTCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                            >  5:246538/1‑146 (MQ=255)                                                        tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                                          >  6:219384/1‑142 (MQ=255)                                                        tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                                          <  5:219384/142‑1 (MQ=255)                                                          aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATc                                                                                  <  3:285251/148‑1 (MQ=255)                                                          aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCg                                                                                 <  7:197194/149‑1 (MQ=255)                                                                 tgctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACGGGCCTTCGTGGCCGTTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGaaa                                                                           <  1:50103/148‑1 (MQ=255)                                                                           cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgcc                                                                >  8:129113/1‑149 (MQ=255)                                                                            tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGACCGCGAATGCCTTTAGCAGATCAATCGTTGAAATAGCGCCAgccg                                                               >  7:285056/1‑149 (MQ=255)                                                                             gAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgc                                                              >  2:297274/1‑149 (MQ=255)                                                                                cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                                              >  1:60717/1‑68 (MQ=17)                                                                                cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                                              <  2:60717/68‑1 (MQ=17)                                                                                 gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAt                                                                                   >  7:253533/1‑124 (MQ=255)                                                                                 gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAt                                                                                   <  8:253533/124‑1 (MQ=255)                                                                                 gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCATAAGCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgcc                                                             >  1:108291/1‑148 (MQ=255)                                                                                      cGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctt                                                     >  8:75362/1‑149 (MQ=255)                                                                                        cTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                                                   <  2:276112/149‑1 (MQ=255)                                                                                                   tAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGa                                        <  7:62287/148‑1 (MQ=255)                                                                                                    aaTTGCTGCGTGCATGAAAAAAGCGGCCTTCGTGGCGGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGGAGCTCAATGGTTGAAATCGCGGCCGCCGCCCATTTGTTTCCA‑GCGGTTGAc                                       >  7:53452/1‑148 (MQ=255)                                                                                                      ttGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCGGATCAATCGtt                                                                               >  1:27697/1‑107 (MQ=255)                                                                                                      ttGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                               <  2:27697/107‑1 (MQ=255)                                                                                                           gcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCg                                >  2:215914/1‑148 (MQ=255)                                                                                                            cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGc                               <  2:158431/148‑1 (MQ=255)                                                                                                            cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCa                              <  5:56570/149‑1 (MQ=255)                                                                                                            cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCCAATGCCATTAGCAGATAAATAGTTTAAATACCGCCAGCCGCCCATTTCATTACA‑GCGGTTGACGATGccccc                              >  6:56570/1‑146 (MQ=255)                                                                                                               gCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGaa                           >  3:157922/1‑149 (MQ=38)                                                                                                                 ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGaa                           >  2:67314/4‑147 (MQ=255)                                                                                                                   gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                                                   >  4:195099/1‑122 (MQ=255)                                                                                                                   aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGaa                           <  3:324684/144‑1 (MQ=255)                                                                                                                   aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGaa                           >  4:324684/2‑145 (MQ=255)                                                                                                                     aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGa                                        >  3:10366/1‑130 (MQ=255)                                                                                                                     aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGa                                        <  4:10366/130‑1 (MQ=255)                                                                                                                                  cGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtctc         <  7:75362/148‑1 (MQ=255)                                                                                                                                   gTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCg        <  7:129113/148‑1 (MQ=255)                                                                                                                                     ggCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCATGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGgt      <  1:67314/149‑1 (MQ=255)                                                                                                                                      gCCGGTTTTTTCATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGgtgt    <  4:157922/149‑1 (MQ=255)                                                                                                                                        cGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGTGTCa  >  3:184144/1‑149 (MQ=255)                                                                                                                                        cGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGTGTCa  <  1:215914/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                               GCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCA‑GCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGTGTCA  >  NC_002947/1243611‑1243894

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: ‑