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

GTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTT  >  NC_002947/278597‑278878
                                                                                                                                                  |                                                                                                                                         
gTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaa                                                                                                                                        >  2:334936/1‑148 (MQ=255)
 tCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                       <  1:128944/149‑3 (MQ=255)
       cAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACATAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTc                                                                                                                                 <  1:241114/149‑1 (MQ=255)
       cAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTc                                                                                                                                 <  3:103001/149‑1 (MQ=255)
       cAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTc                                                                                                                                 <  2:190318/149‑1 (MQ=255)
          cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                              <  1:334936/149‑1 (MQ=255)
                 cgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTTATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGGTACAGAGGTGCGGGGCTCTTCGCGGGTAAACCCGCTCCTACACACAt                                                                                                                       <  5:310158/149‑1 (MQ=255)
                   aCAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                                     >  2:158256/1‑149 (MQ=255)
                                 gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgtg                                                                                                        >  4:279397/1‑148 (MQ=255)
                                 gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgtg                                                                                                        >  3:145088/1‑148 (MQ=255)
                                 gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                         <  2:77092/147‑1 (MQ=255)
                                 ggccggtggtccggtCTAAATACCCTAATGATTTATAAAATTCTTATGAACTACCTTTTTGGAATATTTATAGAAGGGGCATTTGTTACAGAGCTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                         <  4:104301/145‑1 (MQ=255)
                                     ggtggtccgctATATATACACTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTg                                                                                                   >  8:76884/1‑149 (MQ=255)
                                              gtATATATAACCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                          <  7:364127/149‑1 (MQ=255)
                                                 tatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGtt                                                                                       <  1:158256/149‑1 (MQ=255)
                                                    ataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACACGATACTACACACATGCCTCAAGATATGTGAACTGTAGGAGCGGGTTTAc                                                                                    >  5:294305/1‑149 (MQ=255)
                                                         cTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCAGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                                >  1:288051/1‑148 (MQ=255)
                                                          tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgcg                                                                               >  5:212110/1‑148 (MQ=255)
                                                             tGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCAGCTCCTACACAAATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAg                                                                           >  3:272533/1‑149 (MQ=255)
                                                              gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAgg                                                                          >  2:54604/1‑149 (MQ=255)
                                                                         tttttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTTCACACATGCCTCAAGATCtgtg                                                                                                        <  2:138338/108‑1 (MQ=255)
                                                                         tttttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCAGCTCCTACACACATGCCTCAAGATCtgtg                                                                                                        >  1:138338/1‑108 (MQ=255)
                                                                            ttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGtt                                                             <  1:54604/148‑1 (MQ=255)
                                                                                gAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGAACAGGTTTACTc                                                        >  1:219307/1‑149 (MQ=255)
                                                                                  actaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgcg                                                                               <  6:101303/124‑1 (MQ=255)
                                                                                  actaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgcg                                                                               >  5:101303/1‑124 (MQ=255)
                                                                                  acaaactTTTTGTAATATTTTTACACGGGGGATTTGTAACAGAGGTGCCGGGCTCTTCGCGGTTAAACCCGCTCCTACACACATGCCTCAAGATCTGTCCACTGTAGGAGGGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCgg                                                      <  4:145088/146‑1 (MQ=255)
                                                                                   ctaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTCGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                     <  3:279397/149‑1 (MQ=255)
                                                                                    taactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTAAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAATGGCCGGAACAGGATTACTCGGAg                                                    >  6:132271/1‑149 (MQ=255)
                                                                                           tttgAATATTTATAGATGGGTTATTTGTTAACGTGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAATATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCa                                              <  6:294305/145‑1 (MQ=38)
                                                                                             ggAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTc                                           >  5:220321/1‑149 (MQ=255)
                                                                                                ttatTTATAGATGGGGTATTTGTTACTGTGGTGCGGGGGTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGCGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCg                                        <  2:250804/148‑1 (MQ=38)
                                                                                                atatTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTccc                                         <  6:212110/148‑1 (MQ=255)
                                                                                                      ataGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGaa                                   <  7:76884/148‑1 (MQ=255)
                                                                                                             ggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATCCCTCAAGATCTGTCCACTGTAGGCGCGGGTTAACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGcc                                                 <  5:30923/127‑1 (MQ=38)
                                                                                                             ggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTCCTACACACAATCCTCAAGATCAGTGCACTGTAGGAGCGGGTTTAACCCCGAAAGGGCCCGCACAAGTTTACTAGGAGGcc                                                 >  6:30923/1‑127 (MQ=255)
                                                                                                                  tttGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGTCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGaa                                   >  1:294097/1‑136 (MQ=255)
                                                                                                                  tttGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGaa                                   <  2:294097/136‑1 (MQ=255)
                                                                                                                          agcggTGCCGGGCTGTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGt              <  4:272533/146‑1 (MQ=255)
                                                                                                                            agGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGATAGGGCCGGATCAGGTTTACGCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTAc            >  7:292629/1‑149 (MQ=39)
                                                                                                                                 ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGctct       >  3:154917/1‑149 (MQ=255)
                                                                                                                                 ccGGGCTCTTCGCGGGTAAACCCGCTCATACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGctct       <  4:154917/149‑1 (MQ=255)
                                                                                                                                      ctctTCGCGGGTAACCCCGCTCCTAAACACATGCCTCAAGATCTGTCCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttcctt  <  6:220321/149‑1 (MQ=255)
                                                                                                                                      ctctGCGCGGTAAACCCCGCCCCAACACACAGCCCACAAGATCTGTCCACTGTCGGCGCGGTTTAACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttcctt  <  1:372705/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                         
GTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTT  >  NC_002947/278597‑278878

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: