Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,837,836	A→G	9.3%	E463G (GAG→GGG)	PP_3388 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,837,836	0	A	G	9.3%	95.5 / 2.5	43	E463G (GAG→GGG)	PP_3388	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (12/27);  new base G (2/2);  total (14/29)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

CCAGCAAGCCGCCGAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCC  >  NC_002947/3837687‑3837975                                                                                                                                                      |                                                                                                                                            ccAGCAAGCCGCCGAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGa                                                                                                                                             >  1:1526862/1‑150 (MQ=37)     cAAGCCGCCGAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGa                                                                                                                                             <  2:953753/146‑1 (MQ=35)     cAAGCCGCCGAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGa                                                                                                                                             >  1:953753/1‑146 (MQ=35)      aaGCCGCCGAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGa                                                                                                                                        <  1:1920625/150‑1 (MQ=35)        gccgccGAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGTCCAGTGCGTGATCGGCATGGAGTTGAGt                                                                                                                                      >  1:1588438/1‑150 (MQ=25)          cgccgAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGGGTTGAGTGa                                                                                                                                    >  1:484860/1‑150 (MQ=25)          cgccgAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGGGTTGAGTGa                                                                                                                                    >  1:484435/1‑150 (MQ=25)           gccgAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGaa                                                                                                                                   <  1:1628860/150‑1 (MQ=35)             cgAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAgg                                                                                                                                 <  2:1094585/150‑1 (MQ=35)               aaGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAgggg                                                                                                                               <  2:899013/150‑1 (MQ=35)                aGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGTCCAGTGCGTGATCGGCATGGAGTTGAGTGAAggggg                                                                                                                              <  2:1588438/150‑1 (MQ=25)                     gcgcCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGa                                                                                                                         >  2:1567783/1‑150 (MQ=35)                      cgcCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGaa                                                                                                                        <  1:279036/150‑1 (MQ=35)                         cAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCAc                                                                                                                           <  2:1143395/144‑1 (MQ=35)                         cAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCAc                                                                                                                           >  1:1143395/1‑144 (MQ=35)                         cAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGt                                                                                                                     <  1:672952/150‑1 (MQ=35)                         cAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGt                                                                                                                     >  2:1901084/1‑150 (MQ=35)                         cAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGAACGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGt                                                                                                                     >  2:1490319/1‑150 (MQ=25)                         cAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCGGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGt                                                                                                                     >  1:1861348/1‑150 (MQ=25)                           aCGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGGGTTGAGTGAAGGGGGCACGAAGGTAc                                                                                                                   <  2:484435/150‑1 (MQ=21)                           aCGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGGGTTGAGTGAAGGGGGCACGAAGGTAc                                                                                                                   <  2:484860/150‑1 (MQ=21)                            cGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGAAAGGAGTTGAGTGAAGGGGGCACGAAGGTACt                                                                                                                  <  1:51162/150‑1 (MQ=18)                             gACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGCCCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTg                                                                                                                 <  2:1856324/150‑1 (MQ=17)                             gACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTg                                                                                                                 <  2:177560/150‑1 (MQ=32)                              aCCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTgg                                                                                                                <  1:1901084/150‑1 (MQ=32)                               ccATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGa                                                                                                                                        <  2:967738/125‑1 (MQ=33)                               ccATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGa                                                                                                                                        >  1:967738/1‑125 (MQ=33)                                cATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGAc                                                                                                              <  2:932040/150‑1 (MQ=32)                                  tACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAggggg                                                                                                                              <  1:1136424/132‑1 (MQ=33)                                  tACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAggggg                                                                                                                              >  2:1136424/1‑132 (MQ=33)                                      cGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACtgg                                                                                                        <  1:673900/150‑1 (MQ=32)                                        cTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACtggtg                                                                                                      <  1:106188/150‑1 (MQ=32)                                        cTGGGCGGCCAGCGAACGCGTATTGGAAGCGCGGGACCGCTGCGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACtggtg                                                                                                      <  2:874575/150‑1 (MQ=11)                                          gggTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGgc                                                                                                    <  1:1188129/150‑1 (MQ=32)                                              ggtCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAg                                                                                                <  2:1370419/150‑1 (MQ=32)                                                    cGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGt                                                                                          <  2:454435/150‑1 (MQ=32)                                                    cGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGt                                                                                          <  2:1374909/150‑1 (MQ=32)                                                     gAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGTg                                                                                         >  2:238260/1‑150 (MQ=32)                                                     gAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGATGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGAGCAGTGACGTg                                                                                         >  2:1001642/1‑150 (MQ=17)                                                       aCGCTTATTGCAAGCGCTGGACCGTTACGACAACGCCGACCACATCAACGGCCTGGCCTGTGCCGAACAGACCGGCCAGGCCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGTGGc                                                                                       <  1:949991/150‑1 (MQ=11)                                                       aCGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGTGGc                                                                                       <  1:876125/150‑1 (MQ=32)                                                       aCGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGTGGc                                                                                       <  2:770616/150‑1 (MQ=32)                                                           ttATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGTGGCAGAc                                                                                   >  1:1231400/1‑150 (MQ=32)                                                            tATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGAGACGTGGCAGAcc                                                                                  >  2:307932/1‑150 (MQ=14)                                                                                                     atcgtcctggcctTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCaggag                                         <  2:1151320/145‑1 (MQ=1)                                                                                                                                            aTCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAgccgcc  <  2:100177/150‑1 (MQ=32)                                                                                                                                                      |                                                                                                                                            CCAGCAAGCCGCCGAAGACCGCGCCAACGACCATACCCGCTGGGTGGTCAGCGAACGCTTATTGCAAGCGCTGGACCGCTACGACAACGCCGACCTCATCAACGGCCTGGCCTTTGCCGAACAGACCGGCCAGTGCGTGATCGGCATGGAGTTGAGTGAAGGGGGCACGAAGGTACTGGACCACTGGTGGCGCAGTGACGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCC  >  NC_002947/3837687‑3837975

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑