Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	278,744	+A	intergenic (‑239/+75)	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	95.7%	140.6 / 0.2	47	intergenic (‑237/+77)	PP_0224/sctC	DszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (24/21);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

TCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑T‑AACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCT  >  NC_002947/278598‑278885                                                                                                                                                    |                                                                                                                                                tCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑T‑aaac                                                                                                                                             >  1:306353/1‑147 (MQ=255)      aCAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCt                                                                                                                                         >  3:13584/1‑149 (MQ=255)       catcTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTc                                                                                                                                        <  3:213793/146‑1 (MQ=255)       cAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCt                                                                                                                                         <  6:174461/148‑1 (MQ=255)         gCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTcc                                                                                                                                       >  5:155336/1‑148 (MQ=255)          cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTa                                                                                                                                     >  8:291549/1‑149 (MQ=255)          cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTa                                                                                                                                     <  7:202529/149‑1 (MQ=255)             aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTaca                                                                                                                                   <  1:313514/148‑1 (MQ=255)             aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTaca                                                                                                                                   <  8:220352/148‑1 (MQ=255)             aTCGCCACAGAAGTGTGCAGGTCCGGTGGTCCCGTATAAATAACCTAATGATTTATAAAATTTTTATTAACTAACTTTTTTGAATATTTATAGAATGGGGATTTTTTACATAGGTGCCCGGGCTCTTCGCGGG‑TAAACCCGCTCCTa                                                                                                                                     <  1:23893/147‑1 (MQ=255)                       aCGTGTGCTGGTCCGGTGGTCCGGTATATATACACTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTATACACATGCCTc                                                                                                                         >  2:25332/1‑148 (MQ=255)                            tgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATc                                                                                                                   >  3:220843/1‑149 (MQ=255)                            tgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAAACCGCTCCTACACACATCCCACAAGATc                                                                                                                   >  2:221952/1‑149 (MQ=255)                              cTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                                  <  4:13584/148‑1 (MQ=255)                                       tggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGTGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                                        <  2:306353/149‑1 (MQ=255)                                             ggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCg                                                                                                  >  7:231259/1‑149 (MQ=255)                                              gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTATACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                                 <  1:25332/149‑1 (MQ=255)                                               tatatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCggg                                                                                                <  4:220843/149‑1 (MQ=255)                                                 tatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt                                                                                               <  2:370678/148‑1 (MQ=255)                                                     taCCCTTATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTTTTACAGAGGTG‑CCGGGCTCTTCGCTGG‑TAAACCCGCTCCTACACACATGCCTCAAAATCTGTGCACTGTAGGAGCGGGTTTAcc                                                                                          <  1:221952/149‑1 (MQ=255)                                                          tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGTAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGa                                                                                                     <  7:266018/133‑1 (MQ=255)                                                          tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGTAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGa                                                                                                     >  8:266018/1‑133 (MQ=255)                                                          tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACt                                                                                                           <  8:223694/127‑1 (MQ=255)                                                          tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgcg                                                                                      >  5:205065/1‑148 (MQ=255)                                                            aTGATTTATAAAATTTTTATGAACTAACTTTTTGTAATATTTATAGAAGGTGTATTTGTTACAGATGTG‑ACGGGCTATTAGCGGG‑TAAACCCGCTCCTAAACACATGCCTCAATATCTGTTCACTTTAGGAGCGGGTTTAACAGAGaaa                                                                                   >  7:359878/1‑149 (MQ=255)                                                                                gAACTAACTATTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTc                                                               >  7:101581/1‑149 (MQ=255)                                                                                 aactaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTc                                                               <  6:155336/148‑1 (MQ=255)                                                                                 aactaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTACTACACAAATGCCTAAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTaatag                                                              >  3:83862/1‑145 (MQ=39)                                                                                 aactaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAAACAGCTCCTACACACATGCCTCAAGCTCTGTGCACTGTAGGAGCGGGGTTACCCGCGAAAGTGCAGGAACAGGTTTACTCg                                                              >  3:89371/1‑149 (MQ=255)                                                                                   ctaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                            <  6:205065/149‑1 (MQ=255)                                                                                      actTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAg                                                           <  7:291549/147‑1 (MQ=255)                                                                                         ttttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACCCCTGCCTCAAGATCTGTGCACTGTAGGAGCGGGGTTAACCGCGAAAGGGCCGGAACAGGGTTACTCGGAGGCCg                                                       >  8:53441/1‑148 (MQ=255)                                                                                              taaTATTTATAGAAGTGGTATTTGTTACAGAGGGG‑CCGGGCTCTTCTCGTG‑TAAACCCGCTCATACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTc                                                  <  4:83862/147‑1 (MQ=38)                                                                                              gAATAATTATAGAAGGGGCATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGGTT‑AAACCGCAAGATCACAAATGCCTCAAGATGTATGAAATGTAGCATCCGGTTTACCCGCGAAAGGGACGGAAAAGGTTTAAACGGAGGcccgactac                                                 >  8:17768/1‑142 (MQ=255)                                                                                                    ttATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATATGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAgg                                            >  3:249149/1‑148 (MQ=255)                                                                                                            gggggATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                                    >  5:221125/1‑148 (MQ=255)                                                                                                              gggATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGGTT‑AACCCCCAACTAAACACATGGCGCAAAATATTTGGACGGTAGGAGCGGGGTTTAACGCGAAAAGGCCGGAAAAGGTTTAATAGGAGGCCCCCATCCCCAGGAACATTCGCa                                  >  4:367630/1‑148 (MQ=255)                                                                                                                        aCAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAAACCCCTCCTACACACAAGCCTCAATAACATTGCAATGTAGGAGCGGGGTTACCCGCCAAAAGGCCGGAACAAGTTTACACGGAGGCCGCAGTACCGAGGAACTTGCGCAAGAACTGCcgc                       >  4:82373/1‑149 (MQ=255)                                                                                                                         cAGAGGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgcg                      >  6:400694/1‑149 (MQ=255)                                                                                                                          agagGTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTGCAACACACAATCCTCAAGATCTGGGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAAAAATGCcccg                      >  2:287175/1‑148 (MQ=38)                                                                                                                              gTG‑CCGGGCTCTTCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCGCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGc                 >  1:39660/1‑149 (MQ=255)                                                                                                                                  ccGGGCTCTTCGCGGG‑TAAACCCGCTCCGACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGctct              >  7:32324/1‑149 (MQ=255)                                                                                                                                        tctTCGCGGG‑TAAAACCGCTCCGACACACATGCCAAAAGATAAGTGCACAGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGAACAGGGTTAAAAGGAGGCCGCAGTCCCGAGGAACATGCGCAAGAACTGCCGCGTTCGCActtcct          >  8:242254/1‑147 (MQ=255)                                                                                                                                          ttCGCGGG‑TAAACCCGCTCCTTCACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTTCCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTACCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCgg      >  2:213826/1‑149 (MQ=255)                                                                                                                                          ttCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGtt                                                                                              <  2:115148/61‑1 (MQ=255)                                                                                                                                          ttCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGtt                                                                                              >  1:115148/1‑61 (MQ=255)                                                                                                                                          ttCGCGGG‑TAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttccttc        <  8:231259/147‑1 (MQ=255)                                                                                                                                              cgGG‑TAAACCCGCTCCTACACACATCCCTCAAGATCTGTGCACTGTCGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCt  <  7:242254/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                TCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTG‑CCGGGCTCTTCGCGGG‑T‑AACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCT  >  NC_002947/278598‑278885

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