Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	278,742:1	+A	92.5%	intergenic (‑237/+77)	PP_0224 ← / ← sctC	DszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	278,742	1	.	A	92.5%	158.5 / 6.4	53	intergenic (‑237/+77)	PP_0224/sctC	DszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base A (26/23);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCT  >  NC_002947/278598‑278877                                                                                                                                                  |                                                                                                                                         tCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑aaac                                                                                                                                      >  7:166438/1‑147 (MQ=255) tCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑aaac                                                                                                                                      <  4:89232/149‑3 (MQ=255)  cGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AA‑acc                                                                                                                                    >  2:29231/1‑149 (MQ=255)     tACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTACGGTATATATACACTAATGATTTATAAAATTTTTATGAACTACCTTTTTGGAATATTTATAGAAGGTGGATATGTTACAGAGGTGAAGGGCTCTACGCGGGT‑AA‑acccg                                                                                                                                  >  8:290203/1‑146 (MQ=255)                   aCAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGc                                                                                                                    >  8:28135/1‑149 (MQ=255)                    cAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGcc                                                                                                                   >  3:224765/1‑149 (MQ=255)                     aGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTGGCGGGTAAA‑CCCGCTCCTACACACATGCCt                                                                                                                  <  1:29231/149‑1 (MQ=255)                              cttgtcaggtggtccggtATATATACACTAATTATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGACGGGGGATTTGTTACAGAGGTGCCGTGCTCTTCGCGTGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCt                                                                                                          <  1:226326/145‑1 (MQ=255)                              cTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCtg                                                                                                         >  5:80751/1‑149 (MQ=255)                              cTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACAAGCCTCAAGATcag                                                                                                         >  2:213312/1‑147 (MQ=255)                               tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGTGATTTGTTACAGAGTTGCCGGGCTCTTCGCGGGTAAA‑CACGCACCTACAGACATGGCTCAAGATCtgt                                                                                                        >  7:137982/1‑149 (MQ=255)                               tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                        >  1:122935/1‑149 (MQ=255)                                 gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                        <  7:28135/147‑1 (MQ=255)                                 gtacggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTTGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                        <  1:241512/144‑1 (MQ=255)                                     ggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTg                                                                                                  <  8:241487/149‑1 (MQ=255)                                              gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                         <  4:224765/149‑1 (MQ=255)                                               tatatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCggg                                                                                        <  2:122935/149‑1 (MQ=255)                                                atatataCGCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGTGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCggg                                                                                        <  7:290203/148‑1 (MQ=255)                                                 tatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGtt                                                                                      >  3:113639/1‑149 (MQ=255)                                                               aTTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCaa                                                                                                               >  2:257004/1‑110 (MQ=255)                                                               aTTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCaa                                                                                                               <  1:257004/110‑1 (MQ=255)                                                                  tataAAATTTTTATCAACAACCTTTTTGGAATTTTTATACACGGTGCAGTTGTTACCGCGGCGCCGGGCTCTTCGTGGGT‑AACCCCCCTCCTCCACACATGCCTCAAGATCTGTGCACGGTCGGCGCGGGTTTACCCGCGAAAGGGcc                                                                      <  8:137982/148‑1 (MQ=255)                                                                    taAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCAGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCATCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCgg                                                                    >  7:239301/1‑148 (MQ=255)                                                                         tttttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAgg                                                              <  6:80751/149‑1 (MQ=255)                                                                          ttttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTGAAGATCTGTGCAGTGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGt                                                             >  6:200159/1‑149 (MQ=39)                                                                             tATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAg                                                               <  2:267569/144‑1 (MQ=255)                                                                             tATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCAGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAg                                                               >  1:267569/1‑144 (MQ=255)                                                                             tATGAACTAACTTTTTGCAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGAGCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGttt                                                           <  5:200159/148‑1 (MQ=255)                                                                                tatctaattTTTTGTAATTTTTATTTATGGGTTTTTTTTTACTGTGGTGCCGGGCTCTTCGCGGTTAAA‑CCCGCTCCTAAAAATTTGCCTAATGATTTTTGCATTTTTGGTGGGTTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTc                                                       <  8:287253/146‑1 (MQ=255)                                                                                 aactaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCggg                                                                                        <  1:40840/115‑1 (MQ=255)                                                                                 aactaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCggg                                                                                        >  2:40840/1‑115 (MQ=255)                                                                                  actaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑ACAGCTACTAAAAAAATTCCTAAATATATGTTCACTGTATGATCGGGTTTTCCCGCGAAAGGGCCGGAACAAGTTTaa                                                         >  7:287253/1‑144 (MQ=255)                                                                                        tttttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑ACCCCTCCAACCCACCAGGCGCAACATGAGTGGACAGTAGGAGGCGGTTTACCCGCGAAAGGGGCGGGACACGGTTAATCGgggggcc                                               >  2:132511/1‑146 (MQ=255)                                                                                              gAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACTTGACTCAAGATCTGTGCACTGTAGGAGc                                                                                           <  8:253067/99‑1 (MQ=38)                                                                                                  atTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCGACACACAAGCCTCAAGATCTGTGGACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGACAGGGTTACTCGGAGGCCGCAGTCCCGAg                                     >  2:212949/1‑149 (MQ=39)                                                                                                     tataGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGaa                                  >  1:220988/1‑149 (MQ=255)                                                                                                      ataGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTc                                          <  3:29192/140‑1 (MQ=255)                                                                                                      ataGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCGCCGAAACACAAGCCGCAAGATCTGTGCACAGTAGGAGCGGGTTTTCCCCCGAAAGGGCCGGCACACGGTTAAACGGAGGGCGCAGGCCAGatcgg                                  >  4:29192/1‑144 (MQ=255)                                                                                                         gAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACtt                               <  8:166438/148‑1 (MQ=255)                                                                                                          aaGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGGCGGTACAGGGTTACTCGGAGGCCGCAGTCCCGAGGAACGTgc                             >  3:210506/1‑149 (MQ=255)                                                                                                            gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTg                              <  7:73566/146‑1 (MQ=255)                                                                                                            gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTg                              >  8:73566/1‑146 (MQ=255)                                                                                                            gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCGCCAACACCCAAACCACAAAAAAAGTGCCAAGTAGGGGCCGGTTTTCCCGCGCAAGCGCCGCGACAGCTTATCTACCAGGCCCCAGGCCCCAGGCAcaaccgc                           >  4:206049/1‑143 (MQ=255)                                                                                                               ggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AACCCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTCGGCGCGGGTTTACCCGCGAAAGGGcc                                                                      <  1:97578/103‑1 (MQ=39)                                                                                                                taTTTGTTACAGTGGTCCGGGGCTTTTCGCTGGTAAA‑CCCGATCCAACACACTTGCCTCAAGAGGTTTGCACTGTTGGTGGGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGa                       <  4:210506/148‑1 (MQ=255)                                                                                                                     gTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCGCCTACCCACATGCCTCAAGATCTGTGCACTGTAGGCGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGGCCCGAGGAACTTGCGCAAGAACTgg                  >  8:176957/1‑148 (MQ=38)                                                                                                                        aCAGAGGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgc               >  1:226719/1‑149 (MQ=255)                                                                                                                          agagGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGGGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgcg              >  5:297869/1‑148 (MQ=255)                                                                                                                           gagGTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTa            <  2:226719/149‑1 (MQ=255)                                                                                                                             ggTGCCGGGCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACg          <  2:220988/149‑1 (MQ=255)                                                                                                                                   gggCGCTTCGCGGTTAAA‑CCCGCTCCTACCCACATGCCTCATGATCTGTGCAGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttc    <  1:213312/149‑1 (MQ=38)                                                                                                                                    ggCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAg                                                   >  4:72706/1‑101 (MQ=255)                                                                                                                                     gCTCTTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTGAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttcct  >  1:231543/1‑149 (MQ=255)                                                                                                                                       tctTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGt                                                             >  6:68970/1‑88 (MQ=255)                                                                                                                                       tctTCGCGGGTAAA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGt                                                             <  5:68970/88‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                         TCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AA‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCT  >  NC_002947/278598‑278877

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