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	.	87.5%	161.9 / 19.1	48	coding (31/417 nt)	PP_5654	phosphotriesterase
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base . (22/20);  total (25/23)
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.76e-01

GACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTAAA  >  NC_002947/4620884‑4621163                                                                                                                                          |                                                                                                                                               gACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaacaa                                                                                                                                     <  8:187820/149‑4 (MQ=255) gACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaa                                                                                                                                        <  3:160943/146‑1 (MQ=255)  aCCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCaaaaaaaaacaaa                                                                                                                                    <  2:306308/149‑5 (MQ=255)  aCCTTCCAACAAGCGAAGTTTTAGCGAATTATTTAAAATAAGCAAACAATAAACCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAATTGTAACATGCAAATTAAAAAAGCTAACATCTGGAGATGGTCCaaaaaaaaacaaa                                                                                                                                    <  1:237424/149‑5 (MQ=255)     ttCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGa                                                                                                                                 >  4:214290/1‑148 (MQ=255)     ttCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCAAACAAACTCCCGTAAAACTCAGAAGTGTAACATGCAAATGAAAAAACCTAACATCAGGAGATGGTCC‑AAAAAAAAACAAAGa                                                                                                                                 >  4:239018/1‑148 (MQ=255)        cAAGAAGCGAATTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTa                                                                                                                              <  7:265736/148‑1 (MQ=255)        cAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTa                                                                                                                              >  7:369972/1‑148 (MQ=255)         aagaagCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTa                                                                                                                              <  6:66083/147‑1 (MQ=255)                     tAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACAt                                                                                                                 <  2:61383/148‑1 (MQ=255)                         gCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAg                                                                                                            <  2:200619/149‑1 (MQ=255)                         gCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTAACACAAACTCACGTAAAAATCAGAAGTGGAAAATGCAAATGAACAAAGCTAAAATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAAAATaa                                                                                                               >  1:233063/1‑146 (MQ=255)                          cGAATTATTTAAAATTAGCAAACAATATAGCCCAACCAAAACTGACACACACTCCGGTAAAACTCAGAAGGGGAACATGCAAATGAACAAAGCTATCATCTGGAGATGGTCC‑ATAAAAAAACAAAGACTAc                                                                                                                             <  2:245110/131‑1 (MQ=255)                          cGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACACACACAAACTCCCGTAAAACTCAGAAGTCGAACATGCAAATGAACAAACCAAACAACTGGAGATGGTCC‑AAAAAAAAACAAAGACAAc                                                                                                                             >  1:245110/1‑131 (MQ=255)                            aaTTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAATCTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGAATACATTAAAGAACATAAAACCAg                                                                                                         >  7:176023/1‑149 (MQ=255)                            aaTTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAg                                                                                                         >  1:268987/1‑149 (MQ=255)                            aaTTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAg                                                                                                         >  7:32045/1‑149 (MQ=255)                             attattTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGCCTACATTAAATAACAAAAAAGCa                                                                                                          >  6:133425/1‑147 (MQ=255)                               tattAAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGTGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGga                                                                                                       <  6:117571/148‑1 (MQ=255)                                      aaTAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATAtttt                                                                                               <  8:369972/149‑1 (MQ=255)                                                  aTAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCa                                                                                   <  4:264844/149‑1 (MQ=255)                                                   tAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCa                                                                                   >  6:355555/1‑148 (MQ=255)                                                          cAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACg                                                                           >  8:280122/1‑149 (MQ=255)                                                             gCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑‑AAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCt                                                                                             >  2:180025/1‑127 (MQ=255)                                                             gCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑‑AAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCt                                                                                             <  1:180025/127‑1 (MQ=255)                                                              cAAAAATTACACAAACTCCCGTAAAACTAATAATTTTAACATGCAATTTAACAAAGATAACATCTGGATATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGaa                                                                        <  5:133425/148‑1 (MQ=255)                                                                aaaCTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACa                                                                      >  4:27326/1‑148 (MQ=255)                                                                 aaCTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACata                                                                    <  3:214290/149‑1 (MQ=255)                                                                 aaCTCACACAAACTCACGTAAAACTCATAAGTGTAACATGCAAATTAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACata                                                                    <  3:239018/149‑1 (MQ=255)                                                                      acacaAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATattatt                                                               >  4:30665/1‑149 (MQ=255)                                                                           aaCTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTggg                                                          <  8:32045/149‑1 (MQ=255)                                                                              tCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACatat                                                                   >  5:114444/1‑137 (MQ=255)                                                                               cccGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACAtatt                                                                  <  6:114444/137‑1 (MQ=255)                                                                                 cttAAAACTCAGAAGTGGAACATTCAAATGAACAAAGCTAACATATTGATATGGTCC‑AAAAAAAAACAAATACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAaata                                                    <  2:233063/147‑1 (MQ=255)                                                                                       aCTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAatat                                                >  5:327340/1‑147 (MQ=255)                                                                                        cTCAGAAGTTGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCaa                                             <  3:351337/149‑1 (MQ=255)                                                                                         tCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCaaa                                            <  7:280122/149‑1 (MQ=255)                                                                                           aGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCaaaaa                                          >  8:111099/1‑149 (MQ=255)                                                                                                  gAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACAtattat                                                                >  1:219508/1‑121 (MQ=255)                                                                                                  gAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACAtattat                                                                <  2:219508/121‑1 (MQ=255)                                                                                                  gAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTgg                                                           <  1:122514/126‑1 (MQ=255)                                                                                                  gAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTgg                                                           >  2:122514/1‑126 (MQ=255)                                                                                                  gAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCa                                    >  6:11628/1‑148 (MQ=255)                                                                                                       tGCAAATGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGc                               <  5:11628/148‑1 (MQ=255)                                                                                                            aTGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATg                           >  2:117235/1‑148 (MQ=255)                                                                                                            aTGAACAAAGCTAACATCTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATg                           >  2:90150/1‑147 (MQ=255)                                                                                                                             cTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCaaaaa                                          >  1:64283/1‑115 (MQ=255)                                                                                                                             cTGGAGATGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCaaaaa                                          <  2:64283/115‑1 (MQ=255)                                                                                                                                gagaTGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAaat                                                     <  1:232575/101‑1 (MQ=255)                                                                                                                                gagaTGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAaat                                                     >  2:232575/1‑101 (MQ=255)                                                                                                                                gagaTGGTCC‑AAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATgg      >  2:242898/1‑148 (MQ=255)                                                                                                                                      tttcAAAAAAAAAACAAAGACTACATTAAAGAACATAAAATCAGTAGATATTTTCTGCTTCACCAAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTaaa  <  5:121975/144‑1 (MQ=255)                                                                                                                                          |                                                                                                                                               GACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGGAGATGGTCCAAAAAAAAAACAAAGACTACATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTAAA  >  NC_002947/4620884‑4621163

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