Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,183,628	G→A	9.5%	G469E (GGG→GAG)	PP_3680 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,183,628	0	G	A	9.5%	114.7 / 3.1	42	G469E (GGG→GAG)	PP_3680	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/16);  new base A (2/2);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

AGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTCGACCCTACACG  >  NC_002947/4183482‑4183775                                                                                                                                                   |                                                                                                                                                    aGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGgc                                                                                                                                                   <  2:260230/149‑1 (MQ=255)   ccGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGgcgc                                                                                                                                                 >  3:111075/1‑149 (MQ=255)   ccGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGgcgc                                                                                                                                                 >  3:106494/1‑149 (MQ=255)     gTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCCCCCCTAGTGGAgcgatc                                                                                                                                               >  1:313678/1‑146 (MQ=255)     gTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCCCCCCTAGTGGAgcgatc                                                                                                                                               <  2:313678/149‑4 (MQ=255)           gCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCaa                                                                                                                                         >  4:250736/1‑149 (MQ=255)                   cGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCt                                                                                                                                 <  3:250736/149‑1 (MQ=255)                       ctctACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTggc                                                                                                                              <  8:62271/148‑1 (MQ=255)                                gCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGaa                                                                                                                    >  3:88572/1‑149 (MQ=255)                                          tggtACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGt                                                                                                          <  1:163321/149‑1 (MQ=255)                                             tACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTc                                                                                                                                           <  1:145440/113‑1 (MQ=255)                                             tACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTc                                                                                                                                           >  2:145440/1‑113 (MQ=255)                                                        ggggATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGgctg                                                                                                                         >  4:82905/1‑120 (MQ=255)                                                        ggggATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGgctg                                                                                                                         <  3:82905/120‑1 (MQ=255)                                                          ggATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGAccc                                                                                          >  1:89639/1‑149 (MQ=255)                                                                tGAACAAAACCTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATAGCTACAATAGGCGGGCGCCATGTCAATCACATCTGTCGGCTGCTGAACGCGCTGGGTAAACCTCATATCACCCTGCTCg                                                                                    >  4:264386/1‑149 (MQ=255)                                                                 gaagaaATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTGat                                                                                                 >  8:37657/1‑135 (MQ=255)                                                                  aagaaATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGAc                                                                                  >  5:85467/1‑149 (MQ=255)                                                                   agaaATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGAcc                                                                                 <  5:275836/149‑1 (MQ=255)                                                                                       tCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                             <  4:106494/149‑1 (MQ=255)                                                                                       tCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                             <  4:111075/149‑1 (MQ=255)                                                                                          ttCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                             <  3:119649/146‑1 (MQ=255)                                                                                          ttCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCaa                                                             >  4:119649/1‑146 (MQ=255)                                                                                          ttCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGt                                                          >  1:139210/1‑149 (MQ=255)                                                                                               gCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTg                                                     >  2:268010/1‑149 (MQ=255)                                                                                                     ggAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCg                                               >  7:149298/1‑149 (MQ=255)                                                                                                          tCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACCCCTCATATGACCCTGTTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGtgtg                                           >  8:189026/1‑148 (MQ=255)                                                                                                          tCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGCCCCTGCTCGACCTAGACGTCGGGCGCCATCAAGGTGGCTGGGGGCGtgtg                                           >  4:10908/1‑148 (MQ=255)                                                                                                           cccGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGGCAATCACTTCTGGCGGCTGCGGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGGGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGaa                                         <  2:188066/149‑1 (MQ=255)                                                                                                            ccGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAAc                                        <  4:88572/149‑1 (MQ=255)                                                                                                                  aCGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGtt                                  <  3:264386/149‑1 (MQ=255)                                                                                                                   cGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTc                                 <  1:268010/149‑1 (MQ=255)                                                                                                                          gCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACt                            >  2:315628/1‑147 (MQ=255)                                                                                                                             gTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCt                       <  2:139210/149‑1 (MQ=255)                                                                                                                                catggcGCCCCCCTAGTGGAGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTc                                 >  5:36847/6‑136 (MQ=255)                                                                                                                                catggcGCCCCCCTAGTGGAGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTc                                 <  6:36847/131‑1 (MQ=255)                                                                                                                                  gATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATg                  >  7:99220/1‑149 (MQ=255)                                                                                                                                      gCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTAtt              <  6:85467/149‑1 (MQ=255)                                                                                                                                      gCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTAtt              >  3:155777/1‑149 (MQ=255)                                                                                                                                       cTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACATCATATGACCCTGCTCGACCTAGACGTGGGGCGACATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTc             >  4:282598/1‑149 (MQ=255)                                                                                                                                        tCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTc             >  6:153661/1‑148 (MQ=255)                                                                                                                                          cACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTCGAc          >  5:232156/1‑149 (MQ=255)                                                                                                                                                  gggCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGTCCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTGGATGTATTCGACCCTACACg  >  5:188658/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AGCCGTCCAGGCATTTCCCGAGCTCTACTTCGCACGCTTGGTGGTACTTGGTGAGGGGGATAGTGAAGAAATCGTGCTACCTCGAATCCTTCAAGCTAAGGGAGCTCCCGTTGACGAATACGCAGTGGCGATCGCTCCACTAGGGGGGCGCCATGTCAATCACTTCTGGCGGCTGCTGAACGCGCTGGGTACACCTCATATGACCCTGCTCGACCTAGACGTGGGGCGCCATCAAGGTGGCTGGGGGCGTGTGAACTATGTTCGTACTCAGCTTGATGTATTCGACCCTACACG  >  NC_002947/4183482‑4183775

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