Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,621,030	(A)10→9	coding (40/417 nt)	PP_5654 →	phosphotriesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,621,021	0	A	.	85.0%	119.1 / 22.2	40	coding (31/417 nt)	PP_5654	phosphotriesterase
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base . (17/17);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCA  >  NC_002947/4620875‑4621153                                                                                                                                                   |                                                                                                                                     aaCTCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCACAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACAAGCAAATGAACAAACCTAACATCTCGAGATGGTCCa                                                                                                                                      >  1:181566/1‑147 (MQ=255)        tctaCCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaa                                                                                                                              <  1:106979/145‑1 (MQ=255)          gACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaac                                                                                                                             >  2:102355/1‑146 (MQ=255)           aCCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGTAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaa                                                                                                                              <  5:88487/145‑1 (MQ=255)           aCCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaa                                                                                                                              <  1:110159/145‑1 (MQ=255)              ttCCAAGAAGCGAAGTTTGATCGAATTATTTAAAATAAGCAAACAATAAAGGCCAATCAATACTCACACAAACTCCCGTAAAACTCTGAAGTGGAACATGCAATTGAACAAAGCTAACATCTGGATATGGTCC‑AAAAAAAAACAAAGa                                                                                                                       <  2:79604/148‑1 (MQ=255)              ttCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGa                                                                                                                       >  2:133299/1‑148 (MQ=255)                      agCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTa                                                                                                               <  5:91457/148‑1 (MQ=255)                         gAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACAtt                                                                                                                >  2:16790/1‑145 (MQ=255)                            gTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTAAGAAGTGGAACATGCAAATGAAAAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAAc                                                                                                         <  1:35293/148‑1 (MQ=255)                                 agCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATaa                                                                                                     <  1:25256/147‑1 (MQ=255)                                  gCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAGGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATaaaa                                                                                                   >  4:146390/1‑148 (MQ=255)                                  gCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATaaaa                                                                                                   >  3:68212/1‑148 (MQ=255)                                     aaTTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCa                                                                                                <  8:146822/148‑1 (MQ=255)                                      attattTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCa                                                                                                >  4:149849/1‑147 (MQ=255)                                                aTAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATAtttt                                                                                     <  3:140903/148‑1 (MQ=255)                                                     cAAACAATAAAGCCCAAGCAAAACTCAAACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGc                                                                                 <  1:85437/148‑1 (MQ=255)                                                      aaaCAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAATATAAAAGCAGGAGATATTTTCTGCt                                                                                >  5:130453/1‑147 (MQ=255)                                                       aaCAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGACCATAAAAGCAGGAGATATTTTCTGCTTc                                                                              >  7:134660/1‑148 (MQ=255)                                                          aaTAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCAc                                                                            >  7:56983/1‑147 (MQ=255)                                                            tAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCa                                                                         <  3:149849/148‑1 (MQ=255)                                                            tAAAGCCCAAGCAAAACTCACACAAACTCCCCTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGCAGATATTTTCTGCTTCACCCa                                                                         >  6:19083/1‑148 (MQ=255)                                                             aaaGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAg                                                                        <  2:138851/148‑1 (MQ=255)                                                                  ccAAGCAAAACTCACACAAACTCCCGTAAACCTCATAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGT‑CAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGa                                                                   <  4:114920/148‑1 (MQ=255)                                                                             tCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACAtatt                                                        >  5:83122/1‑148 (MQ=255)                                                                                 acaAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGCCGGAACATattatt                                                     >  4:88757/1‑148 (MQ=255)                                                                                 acaAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATattatt                                                     >  4:90890/1‑148 (MQ=255)                                                                                 acaAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAA‑AAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTAttt                                                    >  2:175187/1‑148 (MQ=255)                                                                                   aaaCTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGACCATATTAGTTTg                                                  >  8:111305/1‑148 (MQ=255)                                                                                          cGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCaaa                                            >  1:132422/1‑147 (MQ=255)                                                                                               aaCTCAGTACTGGAACATGCAAATGAACAAAGCTAACATCTGGTGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAataata                                       <  8:5967/147‑1 (MQ=255)                                                                                               aaCTCAGAAGTGGAACATGCATATGAACAAAGCTAACATCTGGAGATTGTCC‑AAAAAAAAACAAAGAATACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAataata                                       <  8:8222/147‑1 (MQ=255)                                                                                               aaCTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAatat                                      <  2:27987/148‑1 (MQ=255)                                                                                               aaCTCAGAAGTGCAACATGCAAATGAACAAAGCTAACATATGGAGATTGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGATTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAatat                                      <  4:140809/148‑1 (MQ=255)                                                                                                 cTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATc                                     >  6:34534/1‑147 (MQ=255)                                                                                                  tCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCaa                                   <  2:54099/148‑1 (MQ=255)                                                                                                         tGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAaataat                            <  4:130870/148‑1 (MQ=255)                                                                                                           gAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCa                          >  8:133625/1‑148 (MQ=255)                                                                                                                         aCAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGc             >  8:155644/1‑147 (MQ=255)                                                                                                                          cAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCt            >  3:134820/1‑147 (MQ=255)                                                                                                                          cAAAGCTAACATCTGCAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTg           <  7:123555/148‑1 (MQ=255)                                                                                                                            aaGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCCCCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGaa         >  2:97370/1‑148 (MQ=255)                                                                                                                                 aaCATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAatat     >  6:107706/1‑147 (MQ=255)                                                                                                                                  aCATCTGGAGATGGGCC‑AAAACAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTc   <  6:107280/148‑1 (MQ=255)                                                                                                                                   cATCTGGATATTGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTTGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCa  <  3:32973/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                     AACTCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCA  >  NC_002947/4620875‑4621153

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