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_002947278,7421.A87.9% 95.4 / 7.2 33intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base A (16/13);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCG  >  NC_002947/278597‑278880
                                                                                                                                                   |                                                                                                                                           
gTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGtt                                                                                                                                             >  5:123983/1‑146 (MQ=255)
gTCGATACAGCTGAGCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGCGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                           <  8:41899/148‑1 (MQ=255)
  cGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                         <  2:128233/148‑3 (MQ=255)
    cTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAacc                                                                                                                                       <  4:31529/147‑1 (MQ=255)
     tACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAacccg                                                                                                                                     >  1:137479/1‑146 (MQ=255)
           tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                <  6:9859/148‑1 (MQ=255)
                       aCGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCt                                                                                                                     >  8:163902/1‑147 (MQ=255)
                        cGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCa                                                                                                                   <  4:178481/148‑1 (MQ=255)
                               tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                             <  1:44668/147‑1 (MQ=255)
                                  tccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgtg                                                                                                          <  3:67889/147‑1 (MQ=255)
                                  tccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                           <  2:98669/146‑1 (MQ=255)
                                      gtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTg                                                                                                     <  5:12771/148‑1 (MQ=255)
                                            cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCCAGCCTAAACACCAGCCTCAAGATGTGTTCACTGTAGGAg                                                                                               >  2:177271/1‑148 (MQ=255)
                                              gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTACTACACACATGCCTCAAGATCTGTGCACTGTAGGcgcg                                                                                             >  6:71441/1‑148 (MQ=255)
                                               tatatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                            <  2:161975/148‑1 (MQ=255)
                                                   tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTAcacaca                                                                                                                          >  1:83499/1‑114 (MQ=255)
                                                    ataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACAt                                                                                                                         <  2:83499/114‑1 (MQ=255)
                                                    ataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGttt                                                                                        <  5:2300/147‑1 (MQ=255)
                                                     taCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCGCAAGCTCTGTGCACTGTAGGAGCGGGTTTAc                                                                                      >  4:172295/1‑148 (MQ=255)
                                                                tttATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTGCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAggg                                                                           >  7:86170/1‑148 (MQ=255)
                                                                  tataAAATTTTTTTGAACTAACTTTTTGTAATATTTATAGACGGGGGATTTGTTACAGTGGTGCGGGGCTCTTCGCGGGTAAACCCGCTCCTACACACTTGCCTCAAGATCTGTGCACTGTAGGAGGGGGTTTACCCGCGAAAGGGc                                                                          <  8:47374/147‑1 (MQ=255)
                                                                  tataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                         >  4:128041/1‑148 (MQ=255)
                                                                     aaaaTTTTTATTAAATAACTTTTTTTAATATTTATAGAAGGGGTATTTTTTCCAGCGGTGCCGGGCTCTGCGCGGTTAAACCCGCTCCAACACACATCCCTAAATATCTGTGCAGTGTCGGAGCGGGTTTACCCGCGAAAGGGCCgg                                                                       <  2:50769/147‑1 (MQ=255)
                                                                       aaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCGCAAGATGTGTGCACTGTAGGAGCGGGGTTTCCCGCGAAAGGGCCGGa                                                                      >  5:157705/1‑146 (MQ=255)
                                                                         tttttATGAACTAACTTTTTGGAATATTTATAGACGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAg                                                                  <  6:96354/148‑1 (MQ=255)
                                                                         tttttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAg                                                                  >  1:122064/1‑148 (MQ=255)
                                                                         tttttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCGACACACATGCCTCAACATATGTGCACAGTAGGAGCGGGTATACCCGCGAAAAGGCCCGAACAg                                                                  >  8:94659/1‑148 (MQ=255)
                                                                                    taactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACCCCTGCCTCAAGATGTGTGCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTGTACTgggc                                                       >  6:55861/1‑147 (MQ=38)
                                                                                      actTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCa                                                                                                                   >  7:79143/1‑86 (MQ=255)
                                                                                       ctTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                                  <  8:79143/86‑1 (MQ=255)
                                                                                            tGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAg                                               >  4:38449/1‑148 (MQ=255)
                                                                                                atatTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCCCACCGACACACATGCCTCAAGATCTGTTCACGGTTGGAGCGGGTTTACCCGCGCAAGGGCCGGAACAGGGTTACACGGAGGCCGCAGGccc                                           >  5:3339/1‑148 (MQ=255)
                                                                                                        aGACGGGGTATTTGTAACTGCGGGGCCGGGCTCTGCGCGGGAAACCCCGCCCAAACACACCGCCCTCAATATCTGGCCACTGTAGCCGGGGGTTTCCCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAAc                                    <  6:123983/147‑1 (MQ=255)
                                                                                                               ggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACCCACATGCCTCAAGATCTGTGCACGGTAGGCGCGGGTTTAACCGCGAAAGGGCCGGAACAGGGTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                               >  1:109142/1‑145 (MQ=255)
                                                                                                                                  cGGGCTCTTCGCGGGAAAACCCGCTCCTACACACATGCCGCAAGATCTGTGCACTGTAGGAGCGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGctct         <  7:94659/148‑1 (MQ=255)
                                                                                                                                         ttCGCGGGTAAACCCGCTCCTACACACCTGCCTAAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGAACAGGTTTACTCGGAGGCCCCAGTCCCGAGGCACTTGCGCAAGAACTGCCGCGTACGCACTTCCTTCg  >  8:97777/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCG  >  NC_002947/278597‑278880

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: