Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,244,369	C→G	10.6%	G39A (GGT→GCT)	rnt ←	ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,244,369	0	C	G	10.6%	155.5 / 6.0	55	G39A (GGT→GCT)	rnt	ribonuclease T
Reads supporting (aligned to +/- strand):  ref base C (20/29);  new base G (4/2);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

TCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTT  >  NC_002947/1244222‑1244516                                                                                                                                                    |                                                                                                                                                     tCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCTATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAaccacc                                                                                                                                                     >  4:147373/1‑149 (MQ=255)   aaGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCgg                                                                                                                                                   >  5:190682/1‑149 (MQ=255)   aaGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCgg                                                                                                                                                   >  1:62384/1‑149 (MQ=255)   aaGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCg                                                                                                                                                    <  1:231027/148‑1 (MQ=255)   aaGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCg                                                                                                                                                    <  6:198397/148‑1 (MQ=255)    aGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGt                                                                                                                                                  <  2:115874/149‑1 (MQ=255)     gTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGtt                                                                                                                                                 <  6:190682/149‑1 (MQ=255)           ccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCAGTGGTTTCAACa                                                                                                                                           <  7:344304/149‑1 (MQ=255)                  gATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATccacc                                                                                                                                     <  2:171234/148‑1 (MQ=255)                     gTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGa                                                                                                                                 >  8:242617/1‑149 (MQ=255)                     gTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGa                                                                                                                                 >  3:49389/1‑149 (MQ=255)                         gcgcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCg                                                                                                                              <  5:89025/148‑1 (MQ=255)                          cgcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGc                                                                                                                            <  5:178731/149‑1 (MQ=255)                                 gAACGGCTCCACTTGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGc                                                                                                                     <  7:102305/149‑1 (MQ=255)                                 gAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATa                                                                                                                       <  7:280142/147‑1 (MQ=255)                                 gAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGc                                                                                                                     >  4:368729/1‑149 (MQ=255)                                        tCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGGGGCAATTTCCATCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGaa                                                                                                              <  6:311631/149‑1 (MQ=255)                                                   gAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCa                                                                                                  >  1:56344/1‑149 (MQ=255)                                                      gTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCg                                                                                               <  1:9659/149‑1 (MQ=255)                                                           tATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATc                                                                                                <  2:45931/143‑1 (MQ=255)                                                           tATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATc                                                                                                >  1:45931/1‑143 (MQ=255)                                                            aTGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGc                                                                                         >  6:27271/1‑149 (MQ=255)                                                             tGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGc                                                                                         >  5:282208/1‑148 (MQ=255)                                                             tGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCg                                                                                        <  5:27271/149‑1 (MQ=255)                                                              gTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGc                                                                                         >  4:209171/1‑147 (MQ=255)                                                              gTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGc                                                                                         >  8:65798/1‑147 (MQ=255)                                                               ttCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGa                                                                                       <  7:342062/148‑1 (MQ=255)                                                               ttCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAACCGCTCG‑GCCATCGGGTGGCGAg                                                                                      >  6:209683/1‑149 (MQ=255)                                                               ttCGGGGAACAGGAAACCTTTCTCGTCAATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAg                                                                                      <  4:302279/149‑1 (MQ=255)                                                                     gAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCgc                                                                                >  2:16712/1‑149 (MQ=255)                                                                         aGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCgctgcc                                                                            <  7:65798/149‑2 (MQ=255)                                                                               ccTTTCTCGTCTATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCACTGCTGCCTTg                                                                      <  3:147373/149‑1 (MQ=255)                                                                                cTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAAGCACCTGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                     <  7:242617/149‑1 (MQ=255)                                                                                cTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                     <  4:49389/149‑1 (MQ=255)                                                                                cTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                     <  1:188548/149‑1 (MQ=255)                                                                                   tctcGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGgt                                                                   <  7:2235/148‑1 (MQ=255)                                                                                       gTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTc                                                                >  2:153497/1‑147 (MQ=255)                                                                                            tGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtcg                                                         >  2:102913/1‑149 (MQ=255)                                                                                                 aTGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                    <  2:62384/149‑1 (MQ=255)                                                                                                     tGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTaa                                                <  1:188130/149‑1 (MQ=255)                                                                                                       aCGGCGGCAATTTCCAGCAATGCAGCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAg                                              <  8:25687/149‑1 (MQ=255)                                                                                                          gcggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                    >  4:70165/1‑140 (MQ=255)                                                                                                          gcggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTc                                                    <  3:70165/140‑1 (MQ=255)                                                                                                          gcggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAggat                                           >  2:18351/1‑147 (MQ=255)                                                                                                                 tttCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGCCtacctgc                                    >  4:11854/1‑143 (MQ=255)                                                                                                                       gCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGt                              <  3:132060/149‑1 (MQ=255)                                                                                                                              aTCGGTAGCGCTGTTGAAACGAGCGGTTTCAACATCCACCAGGACCGGCAGATAGCCGCTGAAGCTCTCG‑TCCATCGGGTGGCGAGGACGGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTtctc                       <  5:209683/149‑1 (MQ=255)                                                                                                                              aTCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGaa                                                                                                              >  6:314371/1‑63 (MQ=255)                                                                                                                              aTCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGaa                                                                                                              <  5:314371/63‑1 (MQ=255)                                                                                                                                   tAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCcag                    >  3:210901/1‑147 (MQ=255)                                                                                                                                         tGTTGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGgcgc             >  1:128894/1‑148 (MQ=255)                                                                                                                                          gTTGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGgcgc             >  6:276235/1‑147 (MQ=255)                                                                                                                                           ttGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGCCCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGgcgc             >  1:109025/1‑147 (MQ=255)                                                                                                                                           ttGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTGTTCGTAAAGGTCTTAGGTGACGGGTTTTTCTCCAGCAGGCGCCAg          >  6:58727/1‑149 (MQ=255)                                                                                                                                           gtcaaaCCCCCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCGCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTc                                        <  6:349228/116‑1 (MQ=255)                                                                                                                                               aaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCa      <  4:210901/149‑1 (MQ=255)                                                                                                                                                   accGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCTGCAGGAGCCAGAGCAGttt  >  6:259994/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCG‑GCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTT  >  NC_002947/1244222‑1244516

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