Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	160,311	A→G	9.6%	L26L (CTT→CTC)	PP_0149 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	160,311	0	A	G	9.6%	142.0 / 5.6	52	L26L (CTT→CTC)	PP_0149	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (3/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

AGGCATGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGT  >  NC_002947/160170‑160459                                                                                                                                              |                                                                                                                                                     aggcaTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAg                                                                                                                                              >  1:649158/1‑150 (MQ=255)       gttgtCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACACCTGCAGCTTGAGCTGCAggcgtcc                                                                                                                                        <  2:892402/150‑6 (MQ=255)                 cAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAggcgggc                                                                                                                              >  2:1094320/1‑150 (MQ=255)                     tGACCCGCACCCGGCGGTCAAGCTGTTGCTTCGCTCTTGGGGTCTGTTTGCCCATCGCCTCGTCTAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGc                                                                                                                          <  2:11820/150‑1 (MQ=255)                        cccGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCTTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGa                                                                                                                       <  2:998268/150‑1 (MQ=255)                         ccGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTGGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAggc                                                                                                                                  >  2:781795/1‑138 (MQ=255)                         ccGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAggc                                                                                                                                  <  1:781795/138‑1 (MQ=255)                          cGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGagc                                                                                                                     >  1:736109/1‑150 (MQ=255)                           gCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGNCGCCTGNAGCTCAAGCTGCNGGTCGGCCAAGGCGGGCGTGCCGagca                                                                                                                    <  1:1000046/150‑1 (MQ=255)                            cACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGagcag                                                                                                                   <  1:1094320/150‑1 (MQ=255)                                 ggcggTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCa                                                                                                              >  2:894045/1‑150 (MQ=255)                                   cggTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAgg                                                                                                            >  2:831856/1‑150 (MQ=255)                                         aGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGCCGTGCCGCGCAGCGTAAGGGCGCaa                                                                                                      >  2:913270/1‑149 (MQ=255)                                          gCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGc                                                                                                     >  2:111012/1‑150 (MQ=255)                                          gCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGc                                                                                                     >  2:397663/1‑150 (MQ=255)                                          gCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGc                                                                                                     <  1:111012/150‑1 (MQ=255)                                           cTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGcc                                                                                                    <  2:477266/150‑1 (MQ=255)                                             gTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCaa                                                                                                  <  1:705546/150‑1 (MQ=255)                                                gCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGAGTGAGGGGGCAGGGAagca                                                                                               >  2:936600/1‑150 (MQ=255)                                                 cTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGc                                                                                                                          <  2:516846/122‑1 (MQ=255)                                                 cTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGc                                                                                                                          >  1:516846/1‑122 (MQ=255)                                                    gAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCa                                                                                           <  1:894045/150‑1 (MQ=255)                                                    gAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCa                                                                                           <  1:913270/150‑1 (MQ=255)                                                     aaGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAg                                                                                          <  1:397663/150‑1 (MQ=255)                                                           ggTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCCGCCAAGCAGCCAGgcacgg                                                                                    <  2:649158/150‑2 (MQ=255)                                                                           tGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGc                                                                    <  1:583753/150‑1 (MQ=255)                                                                             gCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCTCGCGCTTCACAGTGCGa                                                                  <  1:831856/150‑1 (MQ=255)                                                                               cTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGt                                                                >  1:783119/1‑150 (MQ=255)                                                                               cTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGt                                                                >  2:121625/1‑150 (MQ=255)                                                                                tcgtcgAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTa                                                               >  1:160826/1‑150 (MQ=255)                                                                                     gAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCAAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGGGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCt                                                          >  1:699127/1‑150 (MQ=255)                                                                                       gCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTc                                                        >  2:227398/1‑150 (MQ=255)                                                                                          gggCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGc                                                     <  2:783119/150‑1 (MQ=255)                                                                                                  gAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCATCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCt                                             >  2:9277/1‑150 (MQ=255)                                                                                                       tttAGTTTGTTGTGGGTCGAGGCCGGACGCCCGCACTTGAATCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCg                                        <  2:699127/150‑1 (MQ=255)                                                                                                         tAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCt                                      <  1:187825/150‑1 (MQ=255)                                                                                                            tttGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGgc                                   >  1:899059/1‑150 (MQ=255)                                                                                                                 ggtgGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGcgccgc                              <  2:515879/149‑1 (MQ=255)                                                                                                                   tgGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCcgcgc                            <  2:899059/150‑1 (MQ=255)                                                                                                                   ggggTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCcgcgc                            <  2:554436/149‑1 (MQ=255)                                                                                                                   ctgggcGATGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCcgcgc                            >  2:1079892/6‑150 (MQ=255)                                                                                                                              cGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGag                                                                                                                      <  2:458430/49‑1 (MQ=255)                                                                                                                              cGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGag                                                                                                                      >  1:458430/1‑49 (MQ=255)                                                                                                                              cGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTAc                 <  1:121625/150‑1 (MQ=255)                                                                                                                              cGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTAc                 >  2:389629/1‑150 (MQ=255)                                                                                                                                 cgaccTGCAGCTTGAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcg              >  1:717128/4‑150 (MQ=255)                                                                                                                                 cgaccTGCAGCTTGAGCTGCAGGCCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcg              >  2:1032999/4‑150 (MQ=255)                                                                                                                                    ccTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGTCTCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctc           <  2:938435/150‑1 (MQ=255)                                                                                                                                      tGCAGCTTGAGCTGCAGGCCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCa                                                       <  1:1095701/104‑1 (MQ=255)                                                                                                                                      tGCAGCTTGAGCTGCAGGCCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCa                                                       >  2:1095701/1‑104 (MQ=255)                                                                                                                                           cTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAgg    <  1:227398/150‑1 (MQ=255)                                                                                                                                             cAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGt  >  2:1073910/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                     AGGCATGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGT  >  NC_002947/160170‑160459

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