Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	802,230	C→T	8.9%	R239C (CGC→TGC)	obg →	GTPase Obg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	802,230	0	C	T	8.9%	126.2 / 2.9	45	R239C (CGC→TGC)	obg	GTPase Obg
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base T (2/2);  total (23/22)
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.33e-01

CCCGTTCACTACCCTGGTGCCAAACCTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGTGAGCGCTGGCTGGTGCTGAACAAGGCCGATAT  >  NC_002947/802082‑802375                                                                                                                                                     |                                                                                                                                                  cccGTTCACTACCCTGGTGCCAAACCTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTc                                                                                                                                                   >  7:116773/1‑149 (MQ=255)  ccGTTCACTACCCTGGTGCCAAACCTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTcg                                                                                                                                                  >  6:253271/1‑149 (MQ=255)   cGTTCACTACCCTGGTGCCAAACCTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTcgc                                                                                                                                                 <  3:222431/149‑1 (MQ=255)   cGTTCACTACCCTGGTGCCAAACCTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTcgc                                                                                                                                                 >  3:7736/1‑149 (MQ=255)            cccTGGTGCCAAACCTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGtgctgctg                                                                                                                                        >  6:238566/1‑149 (MQ=255)                     cAAACCTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTAGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCtggtg                                                                                                                               >  6:197827/1‑149 (MQ=255)                      aaaCCTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCtggtgg                                                                                                                              <  2:89189/149‑1 (MQ=255)                      aaaCCTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCtggtg                                                                                                                               >  6:45853/1‑148 (MQ=255)                          cTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTTGCGTGCTGCTGCACCTGGTGGACAt                                                                                                                          >  8:59143/1‑149 (MQ=255)                          cTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTTGCGTGCTGCTGCACCTGGTGGACAt                                                                                                                          >  8:193759/1‑149 (MQ=255)                                    tCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGCTCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTg                                                                                                                >  5:179883/1‑149 (MQ=255)                                      aGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGa                                                                                                              <  5:45853/149‑1 (MQ=255)                                            gACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGaaa                                                                                                         <  8:79132/148‑1 (MQ=255)                                            gACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAag                                                                                                        <  1:208688/149‑1 (MQ=255)                                            gACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAag                                                                                                        >  1:166294/1‑149 (MQ=255)                                             aCCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAag                                                                                                        >  8:153000/1‑148 (MQ=255)                                                  tGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGccc                                                                                                  >  7:115311/1‑149 (MQ=255)                                                  tGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGccc                                                                                                  >  6:151732/1‑149 (MQ=255)                                                    gAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGccc                                                                                                  <  5:151732/147‑1 (MQ=255)                                                    gAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCg                                                                                                 <  5:197827/148‑1 (MQ=255)                                                     aaGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGc                                                                                               <  8:116773/149‑1 (MQ=255)                                                         gCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGAt                                                                                           <  8:115311/149‑1 (MQ=255)                                                              gTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGAtatcc                                                                                       <  4:7736/148‑5 (MQ=255)                                                                        aCATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAagc                                                                                                       <  8:122056/122‑1 (MQ=255)                                                                        aCATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAagc                                                                                                       >  7:122056/1‑122 (MQ=255)                                                                          aTCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGt                                                                          <  5:253271/149‑1 (MQ=255)                                                                             ccAGGCCTGATCGAAGGCGCTTCCGAAGCTGCCGGCCTGGGTATCCGCTTCATCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGAAGCCCGGCCGATGCCGCCGAAGTCATCGTCaa                                                                       >  3:219663/1‑149 (MQ=255)                                                                                    tGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCt                                                                 >  6:176524/1‑148 (MQ=255)                                                                                       tCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTTGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGAcc                                                             <  7:59143/149‑1 (MQ=255)                                                                                       tCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTTGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGAcc                                                             <  7:193759/149‑1 (MQ=255)                                                                                              cgcTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCa                                                      <  7:145308/149‑1 (MQ=255)                                                                                                  tCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGccc                                                  >  3:144682/1‑149 (MQ=255)                                                                                                   ccGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGaa                                                                                  <  4:256850/116‑1 (MQ=255)                                                                                                             cGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAg                                                                   >  1:194838/1‑121 (MQ=255)                                                                                                             cGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAg                                                                   <  2:194838/121‑1 (MQ=255)                                                                                                                ccTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGc                                    <  7:153000/149‑1 (MQ=255)                                                                                                                   gggTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGc                                                    >  1:104996/1‑130 (MQ=255)                                                                                                                   gggTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGc                                                    <  2:104996/130‑1 (MQ=255)                                                                                                                    ggTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGCGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGt                                  >  6:140212/1‑147 (MQ=255)                                                                                                                    ggTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGt                                  >  6:136486/1‑147 (MQ=255)                                                                                                                    ggTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGTGa                                <  1:251871/149‑1 (MQ=255)                                                                                                                         ccGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGa                                                               >  4:97715/1‑113 (MQ=255)                                                                                                                         ccGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGa                                                               <  3:97715/113‑1 (MQ=255)                                                                                                                         ccGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGTGAGCgct                           >  3:83698/1‑149 (MQ=255)                                                                                                                          cGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGTGAGCgctg                          >  4:41865/1‑149 (MQ=255)                                                                                                                                 tCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGTGAGCGCTGGCTGGTg                   >  7:101949/1‑149 (MQ=255)                                                                                                                                 tCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGTGAGCGCTGGCTGGTg                   >  6:273040/1‑149 (MQ=255)                                                                                                                                   aaGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGTGAGCGCTGGCTGGTGc                  <  5:238566/148‑1 (MQ=255)                                                                                                                                                   cTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGTGAGCGCTGGCTGGTGCTGAACAAGGCCGatat  <  5:176524/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCCGTTCACTACCCTGGTGCCAAACCTGGGCGTGGTCAGCGTCGACCGCTGGAAGAGCTTCGTCATCGCCGACATCCCAGGCCTGATCGAAGGCGCTTCCGAAGGTGCCGGCCTGGGTATCCGCTTCCTCAAGCACCTGGCGCGTACTCGCGTGCTGCTGCACCTGGTGGACATTGCACCGCTGGACGAAAGCAGCCCGGCCGATGCCGCCGAAGTCATCGTCAACGAGCTGACCCGCTTCAGCCCGTCGCTGGCCGAGCGTGAGCGCTGGCTGGTGCTGAACAAGGCCGATAT  >  NC_002947/802082‑802375

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