Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	403,271	G→C	7.8%	intergenic (‑74/+95)	msrA ← / ← PP_0337	peptide methionine sulfoxide reductase MsrA/sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	403,271	0	G	C	7.8%	155.0 / 2.8	50	intergenic (‑74/+95)	msrA/PP_0337	peptide methionine sulfoxide reductase MsrA/sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base C (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

TCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTC  >  NC_002947/403124‑403413                                                                                                                                                        |                                                                                                                                               tCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCGCCGGGCATGCCCGCTCCCACAGGTACAg                                                                                                                                                <  1:187785/149‑1 (MQ=255)    cgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGGACAGca                                                                                                                                              >  4:115608/1‑147 (MQ=255)     ggCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑TGGGCATGCCCGCTCCCACAGGTACAGCACAg                                                                                                                                           <  1:18420/149‑1 (MQ=255)         aggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGttt                                                                                                                                       <  2:212520/149‑1 (MQ=255)             aggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCa                                                                                                                                     >  8:279855/1‑147 (MQ=255)              gggCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAg                                                                                                                                  >  6:294989/1‑149 (MQ=255)               ggCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAgg                                                                                                                                 <  5:265501/149‑1 (MQ=255)                   tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgt                                                                                                                             >  7:256612/1‑149 (MQ=255)                           aGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTcc                                                                                                                     <  7:7502/149‑1 (MQ=255)                                    gACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTACCACAGGTACAGCAAAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGc                                                                                                            >  3:239044/1‑149 (MQ=255)                                     aCGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCCTT‑TCG‑CGGGCACGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGa                                                                                                              <  4:212476/148‑1 (MQ=17)                                        tttttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                        >  5:329665/1‑149 (MQ=255)                                               tCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCcg                                                                                                 >  3:216295/1‑149 (MQ=255)                                                 aCCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCaaa                                                                                                                                   <  2:318720/113‑1 (MQ=255)                                                 aCCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCaaa                                                                                                                                   >  1:318720/1‑113 (MQ=255)                                                 aCCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCcgc                                                                                                >  3:303000/1‑148 (MQ=255)                                                    aGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                            <  3:262694/149‑1 (MQ=255)                                                      gATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAgg                                                                                          <  3:117475/149‑1 (MQ=255)                                                           cctaaCGCAGGACCATGGCCGTTCCTCTGTTTTTATTCAGTGGGGTTGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAACGGTGGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGccc                                                                                      <  4:303000/145‑1 (MQ=38)                                                           ccGACCGCCGGCCCCTGGCGGTTCCTCTGTGTTTTTGCTGGGGGGTCGGGGTTTGGGCC‑‑T‑TCG‑CGGCCAGGCCCGCACCCACTGGTACCCCACAGGTTTCAATGGTGGGGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGccc                                                                                      <  6:203594/148‑1 (MQ=255)                                                                 gCAGGTCCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                        >  7:44308/1‑124 (MQ=255)                                                                 gCAGGTCCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                        <  8:44308/124‑1 (MQ=255)                                                                     gACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTgg                                                                                                                <  7:215217/112‑1 (MQ=255)                                                                     gACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTgg                                                                                                                >  8:215217/1‑112 (MQ=255)                                                                      acaaTGGCCGTTCGTTTTTGTTTTTTCCGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCta                                                                          <  4:216295/146‑1 (MQ=37)                                                                        cATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCt                                                                              <  7:257067/143‑1 (MQ=255)                                                                        cATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCt                                                                              >  8:257067/1‑143 (MQ=255)                                                                          tGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TGTCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAg                                                                       >  1:257419/1‑149 (MQ=255)                                                                           ggCCTTTCCTCTTTGTTTATTCTGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACTGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGt                                                                      <  4:239044/148‑1 (MQ=255)                                                                                   ctctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACATGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGa                                                             >  5:25728/1‑149 (MQ=255)                                                                                          tttATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGc                                                       <  7:198077/148‑1 (MQ=255)                                                                                                 aGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGt                                                 >  2:72449/1‑147 (MQ=255)                                                                                                 aGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGAGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGAGTAAGCCGCTGt                                                 >  7:151859/1‑147 (MQ=255)                                                                                                           gggTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCggg                                     >  6:85466/1‑149 (MQ=255)                                                                                                           gggTTTGGGCC‑‑T‑TCG‑CCGGCATGCCCCCTACCACAGGTACAGCACATGTTTCAAATGTTGTGGTGTTACTGTGGGAACTGGCGTGCCCCCCAAACGGCCCGAAAGGCTTCTATAGTGCCCGAGGATTACGCAGCGGTCAGGCCACCcgg                                     >  8:312683/1‑149 (MQ=255)                                                                                                                tGGGCC‑‑T‑TCG‑CGGGGATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCAcg                                >  7:311147/1‑149 (MQ=255)                                                                                                                   gCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGCACCGCACAGGTTTCAACGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCgg                             <  4:209539/149‑1 (MQ=255)                                                                                                                                   aTGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAg                 >  6:220716/1‑149 (MQ=255)                                                                                                                                       ccGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAg                                                     <  2:305911/109‑1 (MQ=255)                                                                                                                                       ccGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGc             >  3:163813/1‑149 (MQ=255)                                                                                                                                       ccGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGc             >  3:244866/1‑149 (MQ=18)                                                                                                                                        cGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGa                            <  4:245746/133‑1 (MQ=255)                                                                                                                                        cGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGa                            >  3:245746/1‑133 (MQ=255)                                                                                                                                         gCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTc           >  6:142164/1‑149 (MQ=21)                                                                                                                                          cTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCg          <  5:220716/149‑1 (MQ=255)                                                                                                                                             ccACAGTTACAGCACTGGTTTCAATGGTTGTGTTGTTCCTGTGGTAGCGGGGGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGc        <  6:25728/148‑1 (MQ=255)                                                                                                                                              cacaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGc        >  2:144/1‑147 (MQ=255)                                                                                                                                              cacaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa      <  8:311147/149‑1 (MQ=255)                                                                                                                                              cacaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTACTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGCCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGc        >  5:291595/1‑147 (MQ=255)                                                                                                                                                 aGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGtt   <  5:85466/149‑1 (MQ=255)                                                                                                                                                  ggTACAGCACAGGTTTCAAAGGTTGTTGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCTCGCGGATAGCGCTTCAGCCGCTCGGCCAGTTc  <  1:144/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                               TCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑T‑TCG‑CGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTC  >  NC_002947/403124‑403413

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