Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,640,106	C→A	14.8%	Q400K (CAG→AAG)	PP_4950 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,640,106	0	C	A	14.8%	116.6 / 6.1	44	Q400K (CAG→AAG)	PP_4950	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (17/20);  new base A (5/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCTGAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGCCGGCCCTTTCGCGG  >  NC_002947/5639965‑5640252                                                                                                                                              |                                                                                                                                                   tACCTGAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTg                                                                                                                                               >  2:304562/1‑147 (MQ=255)      gAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGAACAGCTTGCTGGCCCCCCACCTGGAGAAcc                                                                                                                                        >  7:201167/1‑149 (MQ=255)       aaCTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCt                                                                                                                                       <  6:66848/149‑1 (MQ=255)               aGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAg                                                                                                                                >  6:296630/1‑148 (MQ=255)                 aCCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAg                                                                                                                             >  5:147162/1‑149 (MQ=255)                 aCCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAg                                                                                                                             <  7:121474/149‑1 (MQ=255)                             cgcgcATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATc                                                                                                                 <  2:97329/149‑1 (MQ=255)                             cgcgcATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATc                                                                                                                 <  7:234711/149‑1 (MQ=255)                              gcgcATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACTTGACCCAGGAGCGCGAGATGATcc                                                                                                                >  1:86318/1‑149 (MQ=255)                              gcgcATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATcc                                                                                                                <  8:24914/149‑1 (MQ=255)                              gcgcATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATcc                                                                                                                >  2:71627/1‑149 (MQ=255)                                  ctgcccGGGCGGCCCCGCGTCGCCACTTTCTGGGTTGAAGCCGGGCCGATCCCGGGGCCGGGGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTc                                                                                                            <  8:201167/145‑1 (MQ=255)                                      ccGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAg                                                                                                                                >  8:281608/1‑125 (MQ=255)                                      ccGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAg                                                                                                                                <  7:281608/125‑1 (MQ=255)                                        gggCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGAc                                                                                                                                         >  5:187220/1‑114 (MQ=255)                                        gggCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGAc                                                                                                                                         <  6:187220/114‑1 (MQ=255)                                             gACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACt                                                                                                  >  3:136802/1‑148 (MQ=255)                                              aCCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACTTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACt                                                                                                  <  2:86318/147‑1 (MQ=255)                                                            tttCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTTCCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGAAAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATCGACAGAGCAGTa                                                                                  >  7:118432/1‑149 (MQ=255)                                                                    tAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTc                                                                           >  5:1783/1‑148 (MQ=255)                                                                      gaagaCGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGa                                                                         <  3:92758/148‑1 (MQ=255)                                                                         gaCGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCa                                                                     >  3:267207/1‑149 (MQ=255)                                                                          aCGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCAt                                                                    >  2:91473/1‑149 (MQ=255)                                                                            gggCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCAttt                                                                  <  7:33702/149‑1 (MQ=255)                                                                                 taTCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCg                                                             <  7:107848/148‑1 (MQ=255)                                                                                   tCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCgg                                                            <  8:77200/148‑1 (MQ=255)                                                                                       ggggCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCAc                                                       >  7:228612/1‑149 (MQ=255)                                                                                         ggCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTg                                                     <  6:1783/149‑1 (MQ=255)                                                                                          gCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTgg                                                    <  1:304562/149‑1 (MQ=255)                                                                                                 aGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAg                                              <  4:267207/148‑1 (MQ=255)                                                                                                      cATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGcctga                                        <  6:147162/149‑1 (MQ=255)                                                                                                      cATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGcctga                                        >  1:324595/1‑149 (MQ=255)                                                                                                             ttCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGcctgaccctga                                  >  6:292850/1‑148 (MQ=255)                                                                                                               cGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACa                                <  2:324595/148‑1 (MQ=255)                                                                                                                    aCAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATtgtg                           <  8:228612/148‑1 (MQ=255)                                                                                                                       gCTTGTACAGCTTGCTGGGCAGCAAGCTGGAGGACCTGACCCAGTAGCGCGAGATTATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCTAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTg                       <  1:71627/149‑1 (MQ=255)                                                                                                                           gTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACt                   >  3:19666/1‑149 (MQ=255)                                                                                                                            tccAGCTGGCTGCCCAGCAAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGc                                            >  5:169739/3‑123 (MQ=255)                                                                                                                            tccAGCTGGCTGCCCAGCAAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGc                                            <  6:169739/121‑1 (MQ=255)                                                                                                                            tccAGCTGGCTGCCCAGCAAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACtt                  >  4:301428/3‑149 (MQ=255)                                                                                                                            tcaaGCTGGCTGCACAGAAAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACACCGCAGTACCGGGTCGAGTCAATTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGa                          >  7:138832/4‑141 (MQ=255)                                                                                                                              cAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTg                 <  4:136802/148‑1 (MQ=255)                                                                                                                              cAGCTGGCTGCCCAGCAAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGc                >  5:158809/1‑149 (MQ=255)                                                                                                                                gCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGCCg              >  1:238395/1‑149 (MQ=255)                                                                                                                                           agccagcTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGCCGGCCCTTTcgcg   <  1:91473/149‑1 (MQ=255)                                                                                                                                             ccagcTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGCCGGCCCTTTCGCgg  >  3:273558/1‑148 (MQ=255)                                                                                                                                             ccagcTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGCCGGCCCTTTCGCgg  <  4:273558/148‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   TACCTGAACTACTCAGACCTGCCAGCGGCGCGCATGACCGGGCTGACCCGCTTCGCCACTTTCTGGGTAGAAGACGGGCAGATCCAGGGGCCGGTGAGTACCATGCGCTTCGATGACAGCTTGTACAGCTTGCTGGGCAGCCAGCTGGAGGACCTGACCCAGGAGCGCGAGATGATCCTGTCGACCAGTACTTATGGACAGCGCAGTACCGGGTCGAGTCATTTGCCCGGGGCACTGGTCAAAGGCCTGACCCTGACATTGTGACTGGACTTGCCGGCCCTTTCGCGG  >  NC_002947/5639965‑5640252

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