Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,409,528	Δ1 bp	8.0%	coding (786/1398 nt)	PP_3893 ←	DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,409,528	0	G	.	8.0%	156.4 / 6.9	50	coding (786/1398 nt)	PP_3893	DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base . (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

AGATCGAGCGGCTTGGCCTTGTGCTGGAAGCGGGCGATGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGACTTGCTCTCGCGAATGGTCAGATATTCG  >  NC_002947/4409385‑4409666                                                                                                                                                |                                                                                                                                             aGAACGAGCGGCTTGGCCTTGTGCTGGAAGCGGGCGATGTTGCGGATTCGGCTCCGCGGCAGGCCGCCCTTCTCGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCGGTCGGCCag                                                                                                                                         >  8:337781/1‑148 (MQ=255)       aGCGGCTTGGCCTTGTGCTGGAAGCGGGCGATGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTTGGTCGGCATCAAGCATGCGCCCGACCGGAGCGGGTATTTTCTGGCTCTTCTCGCCGTCGGACATGGCttt                                                                                                                                   >  6:265656/1‑149 (MQ=255)        gCGGCTTGGCCTTGTGCTGGAAGCGGGCGATGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTg                                                                                                                                  >  6:59480/1‑149 (MQ=255)           gCTTGGCCTTGTGCTGGTAGCGGGCGATGTTGCGGATTCGGCTCAGCGGCAGGCCGGCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGGCCGACCGCAGCGGTTATTTTCTGGTTCTTCTCGCCGTCGGCCATGGCTTTGCCg                                                                                                                               <  2:256593/149‑1 (MQ=255)           gCTTGGCCTTGTGCTGGAAGCGGGCGATGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATAGCTTTGCCg                                                                                                                               >  4:299009/1‑149 (MQ=255)                     gtgCTGGAAGCGGGCGATGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGAt                                                                                                                       <  6:244138/147‑1 (MQ=255)                        cTGGAAGCGGGCGATGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGc                                                                                                                >  2:136970/1‑149 (MQ=255)                          ggAAGCGGGCGATGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCg                                                                                                              >  8:107796/1‑149 (MQ=255)                             aGCGGGCGATGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCggtg                                                                                                           <  5:59480/149‑1 (MQ=255)                                   cGATGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACt                                                                                                     >  8:19540/1‑149 (MQ=255)                                     aTGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATGTTCACGCCGGCGGCCATGGGTTTGCGGGTTACGTTG‑‑TTGCGGGTGGTAATTgc                                                                                                   >  2:344928/1‑148 (MQ=255)                                        ttGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTgacga                                                                                                >  1:213313/1‑149 (MQ=255)                                        ttGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGCTCGGCATCACCCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTgacga                                                                                                >  6:128896/1‑149 (MQ=255)                                          gCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGacgacg                                                                                              >  4:147151/1‑149 (MQ=255)                                              aTTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGGCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGa                                                                                          <  1:251898/149‑1 (MQ=255)                                                 cGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAgg                                                                                        <  7:19540/148‑1 (MQ=255)                                                   gCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCa                                                                                     <  5:128896/149‑1 (MQ=255)                                                   gCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCa                                                                                     <  1:136970/149‑1 (MQ=255)                                                   gCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCa                                                                                     <  2:213313/149‑1 (MQ=255)                                                             aGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACTCCGAGGCCAGGCTGCGCtt                                                                           >  2:300197/1‑149 (MQ=255)                                                               gccgccCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTC‑GCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGc                                                                        >  8:235289/1‑149 (MQ=255)                                                               gccgccCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTC‑GCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGc                                                                        >  8:235287/1‑149 (MQ=255)                                                                         tGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCCATG‑‑TAGCCGGTGGTAACTGACGCCGCCGAGGCCAGGCTGCGCTTGGCCAAcactat                                                               >  8:318776/1‑147 (MQ=255)                                                                           gCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCttt                                                               <  2:103502/147‑1 (MQ=255)                                                                           gCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑ATGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTATgg                                                             >  2:119979/1‑148 (MQ=255)                                                                              gATGCGTAGGTCGGCATCACGCATGCGCCCGACCGAAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTgccag                                                          >  4:145111/1‑149 (MQ=255)                                                                                      ggTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTC‑GCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCAtctcaa                                                 <  7:235289/149‑3 (MQ=255)                                                                                      ggTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTC‑GCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCAtctcaa                                                 <  7:235287/149‑3 (MQ=255)                                                                                       ggccgCATGAAGCATGCGGCCGACCGGATCGGTAATTTTCTGGCTCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCa                                                 <  1:344928/145‑1 (MQ=255)                                                                                            cATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTTGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCaga                                               >  1:321245/1‑146 (MQ=255)                                                                                              tCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGaa                                          >  5:224328/1‑149 (MQ=255)                                                                                                 cGCATGCGCCGGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATGGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAacac                                       <  7:318776/149‑1 (MQ=255)                                                                                                  gCATGCGCCCGACCACAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAacac                                       >  5:277161/1‑148 (MQ=255)                                                                                                       cgcCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAggg                                  <  7:107796/148‑1 (MQ=255)                                                                                                       cgcCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAggg                                  >  5:87789/1‑148 (MQ=255)                                                                                                          ccGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGa                              >  6:9795/1‑149 (MQ=255)                                                                                                             aCCGCAGAGGTAATTTTTTGGTTATTATCGCCGTCGGGCATGGCTTTGCGGGTTTCTATGCCTTGCCGGTGGTTTTTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGAc                             <  1:300197/149‑1 (MQ=255)                                                                                                                     ggTAATTTTCTGGATCTTCTCGCCGTCGGCCATAGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGACTTGCTCTcgc                   <  3:299009/149‑1 (MQ=255)                                                                                                                         aTTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGACTTGCTCTCGCGaa                <  3:145111/148‑1 (MQ=255)                                                                                                                          ttttCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTgccagcc                                                        >  5:363708/1‑107 (MQ=255)                                                                                                                          ttttCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTgccagcc                                                        <  6:363708/107‑1 (MQ=255)                                                                                                                             tCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGACTTGCTCTCGCGAATGGt            <  3:147151/148‑1 (MQ=255)                                                                                                                              cTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGACTTGCTCTCGCGAATGGTc           <  3:149438/148‑1 (MQ=255)                                                                                                                               tGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGACTTGCTCTCGCGAATGGt            >  2:310819/1‑146 (MQ=255)                                                                                                                                 gATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGaa                                          >  2:207223/1‑114 (MQ=255)                                                                                                                                 gATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGaa                                          <  1:207223/114‑1 (MQ=255)                                                                                                                                    cttctCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGAc                             <  8:212770/124‑1 (MQ=255)                                                                                                                                    cttctCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGAc                             >  7:212770/1‑124 (MQ=255)                                                                                                                                       ctcGCCGTTGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGACTTGCTCTCGCGAATGGTCAGatat     <  2:321245/145‑1 (MQ=255)                                                                                                                                       ctcGCCGTCGGCCATGGCTTGGCCGGTATCGATG‑‑TTGCCGGCGGTAACTGACGACGCCGAGGCCAGGCTGCGGTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGACTTGCTCTCGCGAATGGTCAGATATTCg  <  7:337781/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                             AGATCGAGCGGCTTGGCCTTGTGCTGGAAGCGGGCGATGTTGCGGATTCGGCTCAGCGGCAGGCCGCCCTTCTGGCAGATGCGTAGGTCGGCATCACGCATGCGCCCGACCGCAGCGGTAATTTTCTGGATCTTCTCGCCGTCGGCCATGGCTTTGCCGGTATCGATG‑‑TTGCCGGTGGTAACTGACGACGCCGAGGCCAGGCTGCGCTTGGCCAACTCTTTGCCAGCCATCTCCAGAGAGAACACCAGGGCAGACTTGCTCTCGCGAATGGTCAGATATTCG  >  NC_002947/4409385‑4409666

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