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.A96.0% 77.6 / ‑3.0 25intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (11/13);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCG  >  NC_002947/278597‑278880
                                                                                                                                                  |                                                                                                                                          
gTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTa                                                                                                                                           <  5:120785/147‑1 (MQ=255)
     tACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACGCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTAGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaaaccg                                                                                                                                     >  4:4561/1‑143 (MQ=255)
     tACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAacccg                                                                                                                                     >  7:49864/1‑146 (MQ=255)
                  gACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCGtcacacca                                                                                                                         >  5:119492/1‑146 (MQ=255)
                       aCGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTc                                                                                                                    <  7:46725/148‑1 (MQ=255)
                         gtgtGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                                  <  2:33964/148‑1 (MQ=255)
                             gCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGAt                                                                                                               <  4:49550/147‑1 (MQ=255)
                                  tccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTGTTTGGAATATTTATAGAAGGGGGAGTTGTTACAGAGGGGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgtg                                                                                                          <  7:5808/147‑1 (MQ=39)
                                   ccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCa                                                                                                        <  2:132949/148‑1 (MQ=255)
                                      gtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACt                                                                                                      >  3:50807/1‑147 (MQ=255)
                                        ggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                                   <  5:69120/148‑1 (MQ=255)
                                          tccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAg                                                                                                  <  3:103094/147‑1 (MQ=255)
                                            cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACAAATGCCTCAAGATCTGTGCACTGTAGGAg                                                                                               >  2:51417/1‑148 (MQ=255)
                                                        ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcg                                                                                   <  8:92232/148‑1 (MQ=255)
                                                        ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcg                                                                                   >  1:84233/1‑148 (MQ=255)
                                                         cTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTGCTACCCACATGCCGCCAGCTCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                                  >  5:118562/1‑148 (MQ=255)
                                                                   ataAATTTTTTATGAACTAACTTTTTGTAATATTTATAGACGGGGGATTTGTTACAGCGGGGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGGGCACTGTCGGCGCGGGTTTACCCGCGAAAGGGcc                                                                         <  6:118562/147‑1 (MQ=38)
                                                                         tttttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCGGTGCCCGGTAGGAGCGGGGTTACCCGCGAAAGGGCCGGAACAg                                                                  >  1:46377/1‑148 (MQ=255)
                                                                              aTGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCGCCTACACACATGCCGCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTa                                                             >  4:93209/1‑148 (MQ=255)
                                                                                  actaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCggg                                                                                           >  7:66816/1‑114 (MQ=255)
                                                                                  actaaatTTTTGTAATTTTTATCGACGGGGTTTTTTTTACCGTGGGGCGGGGCGGTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTCCAGTGTAGGAGCGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCg                                                         <  8:50667/148‑1 (MQ=255)
                                                                                   ctaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt                                                                                          <  8:66816/114‑1 (MQ=255)
                                                                                    taactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTCCTACACACATGCCTCAAGATCTGGGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                       >  6:25207/1‑148 (MQ=255)
                                                                                               aaTATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTc                                                          >  5:111293/1‑134 (MQ=255)
                                                                                                atatTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCg                                                         <  6:111293/134‑1 (MQ=255)
                                                                                                                          agagGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgc                  <  4:50807/147‑1 (MQ=255)
                                                                                                                                         ttCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGGTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGCGGCCGCCGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCg  >  7:31947/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: