Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,996,944	A→C	19.1%	M88R (ATG→AGG)	PP_3525 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,996,944	0	A	C	19.1%	76.7 / 18.3	42	M88R (ATG→AGG)	PP_3525	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/14);  new base C (4/4);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

CCAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAGGCCTGG  >  NC_002947/3996795‑3997093                                                                                                                                                      |                                                                                                                                                      ccAGCACGCCCTGCTGCAGTTGCAACCACTGCACGTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCa                                                                                                                                                       >  1:661624/1‑150 (MQ=255)  cAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCAt                                                                                                                                                      >  1:445587/1‑150 (MQ=255)    gCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATgc                                                                                                                                                    >  2:382724/1‑150 (MQ=255)    gCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATgc                                                                                                                                                    >  2:172848/1‑150 (MQ=255)                        aaCCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAAt                                                                                                                                >  1:576088/1‑150 (MQ=255)                          ccACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTGa                                                                                                                                               <  2:469623/133‑1 (MQ=255)                          ccACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTGa                                                                                                                                               >  1:469623/1‑133 (MQ=255)                           cACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTGatcagca                                                                                                                                         <  2:472627/138‑7 (MQ=255)                           cACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTGatcagca                                                                                                                                         >  1:472627/1‑132 (MQ=255)                              tGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcg                                                                                                                          <  1:94600/150‑1 (MQ=255)                                 aCCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTGa                                                                                                                                               <  2:679103/126‑1 (MQ=255)                                 aCCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTGa                                                                                                                                               >  1:679103/1‑126 (MQ=255)                                            ctGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAg                                                                                                            <  2:493017/150‑1 (MQ=255)                                                   gggTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGccc                                                                                                     >  2:508701/1‑150 (MQ=255)                                                      tACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGATGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTg                                                                                                  <  1:382724/150‑1 (MQ=255)                                                             tGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCTGCGCGTGCTGCCa                                                                                                             >  2:329121/1‑132 (MQ=255)                                                             tGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCa                                                                                                             <  1:329121/132‑1 (MQ=255)                                                              ggCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTc                                                                                          <  2:538453/150‑1 (MQ=255)                                                                           tCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTgccgcc                                                                             >  1:356797/1‑150 (MQ=255)                                                                             aCGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCtt                                                                           >  1:367347/1‑150 (MQ=255)                                                                                          aCCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCa                                                              <  2:661624/150‑1 (MQ=255)                                                                                             cGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcggc                                                                                                                        <  1:455435/89‑1 (MQ=255)                                                                                             cGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcggc                                                                                                                        >  2:455435/1‑89 (MQ=255)                                                                                                   cTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGCGACCGTCCata                                                     >  1:361989/1‑150 (MQ=255)                                                                                                    tGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCGTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCatat                                                    >  1:698616/1‑150 (MQ=255)                                                                                                    tGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCatat                                                    >  1:125597/1‑150 (MQ=255)                                                                                                    tGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCatat                                                    >  1:224360/1‑150 (MQ=255)                                                                                                     gTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATc                                                   >  2:107704/1‑150 (MQ=255)                                                                                                     gTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCCCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATc                                                   <  1:508701/150‑1 (MQ=255)                                                                                                      tAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATcg                                                  >  2:614574/1‑150 (MQ=255)                                                                                                           gcgcATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGTCCTGAATGTCGGCCACGAATACATCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCgcgcgc                                             >  1:139140/1‑150 (MQ=255)                                                                                                           gcgcATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCgcgcgc                                             <  1:172848/150‑1 (MQ=255)                                                                                                                  gCCGGGACCGCGAACAATTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGt                                      >  1:151723/1‑150 (MQ=255)                                                                                                                                   tttGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGc                     <  2:224360/150‑1 (MQ=255)                                                                                                                                     tGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCAc                   <  2:139140/150‑1 (MQ=255)                                                                                                                                        tGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCaa                >  1:358151/1‑150 (MQ=255)                                                                                                                                        tGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGACACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCaa                >  2:456582/1‑150 (MQ=255)                                                                                                                                          ttCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCtt                                                                           <  1:176607/88‑1 (MQ=255)                                                                                                                                          ttCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCtt                                                                           >  2:176607/2‑89 (MQ=255)                                                                                                                                             aGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCgtgt                         <  1:163538/136‑1 (MQ=255)                                                                                                                                             aGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCgtgt                         >  2:163538/1‑136 (MQ=255)                                                                                                                                               cgcAGGCATTCGCTTAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTa         <  2:246579/150‑1 (MQ=255)                                                                                                                                                gcAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAg        <  2:151723/150‑1 (MQ=255)                                                                                                                                                      aTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAGGCCTgg  >  2:391623/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CCAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAGGCCTGG  >  NC_002947/3996795‑3997093

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