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.A92.3% 116.7 / 3.8 39intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base A (18/18);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑A‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCC  >  NC_002947/278607‑278876
                                                                                                                                        |                                                                                                                                        
cTGATCGCGACAGACGTGTGCTGGTCCGGGGGTCCGGTTTATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTa                                                                                                                            <  7:190388/149‑1 (MQ=255)
       cgtcAGACGTGTGCTGGTCCTGTGGTCCGTTATATATACCCTAATGATTTATATAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTTCAGAGGTGCCGGGCTCTTCGCGGGT‑ATA‑CCCGCTCCTAGAAACAt                                                                                                                     <  2:345732/146‑1 (MQ=255)
                  tgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATc                                                                                                          <  1:48071/149‑1 (MQ=255)
                  tgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATc                                                                                                          >  3:148834/1‑149 (MQ=255)
                     tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTaca                                                                                                                          >  8:106472/1‑130 (MQ=255)
                     tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTaca                                                                                                                          <  7:106472/130‑1 (MQ=255)
                       gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGc                                                                                                     >  3:119137/1‑149 (MQ=255)
                        tccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCa                                                                                                    >  1:121676/1‑149 (MQ=255)
                              ggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTc                                                                                                                <  1:119325/131‑1 (MQ=255)
                              ggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTc                                                                                                                >  2:119325/1‑131 (MQ=255)
                              ggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAg                                                                                              >  2:262836/1‑149 (MQ=255)
                                  cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                       <  8:84522/136‑1 (MQ=255)
                                  cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                       >  7:84522/1‑136 (MQ=255)
                                    gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCaa                                                                                                              <  8:165112/127‑1 (MQ=255)
                                    gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCaa                                                                                                              >  7:165112/1‑127 (MQ=255)
                                              ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAATATATGTGCACTGTAGGAGCGGGTTTACCcgc                                                                              >  4:59674/1‑149 (MQ=255)
                                              ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                              <  3:47893/149‑1 (MQ=255)
                                               cTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                              >  2:61041/1‑148 (MQ=255)
                                                                  ttattaACTAACTTTTTGTAATATTTATAGACGGGGTATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTTTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGtt                                                           <  2:121676/148‑1 (MQ=39)
                                                                   tATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGttt                                                          <  3:59674/148‑1 (MQ=255)
                                                                     tGAACTAACTTTTTGGAAAATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCGCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTAc                                                        <  5:109252/148‑1 (MQ=255)
                                                                          taactTTTTGGAATTATTTTAGACGGGGGATTTGTTACAGTGGTGTCGGGCTCTTCGCGGGTAA‑A‑CCCGTTCTTACAAAAATTCCTAAAGATCTGTGCACTGTAGTAGAGGTTTAACCAGCGAAAGGGCGGGAACAGGTTTACTCGGAg                                                  <  2:271620/149‑1 (MQ=255)
                                                                              tttttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                                    >  1:109171/1‑127 (MQ=255)
                                                                              tttttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                                    <  2:109171/127‑1 (MQ=255)
                                                                                      atatTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCg                                      <  1:61041/149‑1 (MQ=255)
                                                                                               gAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTg                             >  5:67624/1‑149 (MQ=255)
                                                                                                  tctgtaTTTGTTACCCCGGTGCCGGGCTCTTCGCGGGT‑A‑ACCCCGCTCCTACACACATGCCTCAACAGCTGTGCACTGTCGCCGCGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                           <  4:63088/143‑1 (MQ=255)
                                                                                                      gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACCCATGCCTCAAGATCTGTGGACTGTAGGCGCGGGTTTTCCCGCGAAAGGGCCGGGACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCagct                      >  3:63088/1‑146 (MQ=255)
                                                                                                       aTTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGaa                     <  4:119137/149‑1 (MQ=255)
                                                                                                            ttACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGcc                >  1:73762/1‑149 (MQ=255)
                                                                                                             tACAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGTAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGcc                <  1:262836/148‑1 (MQ=255)
                                                                                                               cAGAGGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGGACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgcg             >  3:268596/1‑149 (MQ=255)
                                                                                                                  agGTGCCGTGCTCTTCGCTGGT‑A‑ACCCCGCTCCTACACACATGCCTAAAGATCTGTGCACTGTAGTAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgc              <  4:231865/145‑1 (MQ=255)
                                                                                                                  agGTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgc              >  3:231865/1‑145 (MQ=255)
                                                                                                                    gTGCCGGGCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGTAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGTGCAAGAACTGCCGCGTACGc        <  4:148834/149‑1 (MQ=255)
                                                                                                                         gggCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttc   >  6:877/1‑149 (MQ=255)
                                                                                                                         gggCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttc   <  2:73762/149‑1 (MQ=255)
                                                                                                                         gggCTCTTCGCGGGTAA‑A‑CCCGCACCTACACACAAGGCTCAAGATCTGTGCACTGTAGGATCGGGTTTACCCGCGAAAGGGCCGGAAACGGTTTACTCGGAGGCCCCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACCCTcttc   >  6:342932/1‑149 (MQ=255)
                                                                                                                          ggCTCTTCGCGGGTAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttcc  <  5:877/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
CTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑A‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCC  >  NC_002947/278607‑278876

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