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

TCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGC  >  NC_002947/278598‑278884
                                                                                                                                                  |                                                                                                                                                 
tCGCTACAGCAGCTAGCGACAGACGTGTGCTGTTCCGGTGTTCCGGTATATATACCATAATTCTTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                              <  2:295664/149‑3 (MQ=255)
  gCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAaccc                                                                                                                                           >  7:28983/1‑148 (MQ=255)
   cTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAacccg                                                                                                                                          >  5:369558/1‑147 (MQ=255)
       aGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTcc                                                                                                                                       <  2:109730/149‑1 (MQ=255)
         cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                     <  8:204928/149‑1 (MQ=255)
         cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                     <  2:270431/149‑1 (MQ=255)
         cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                     <  3:263342/149‑1 (MQ=255)
                cgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACAt                                                                                                                              >  5:282674/1‑149 (MQ=255)
                  aCAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                                            <  8:151228/149‑1 (MQ=255)
                        gtgtGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACGCACATGCCTCaa                                                                                                                       <  6:244037/148‑1 (MQ=255)
                           tgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATc                                                                                                                   >  5:107556/1‑149 (MQ=255)
                             cTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                                  <  6:26637/148‑1 (MQ=255)
                              tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                                >  3:61395/1‑149 (MQ=255)
                              tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                                >  8:29597/1‑149 (MQ=255)
                                gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                                <  6:233888/147‑1 (MQ=255)
                                   cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCAc                                                                                                            <  6:107556/148‑1 (MQ=255)
                                    ggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTg                                                                                                          >  3:194738/1‑149 (MQ=255)
                                      tggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGCAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                                        <  6:307774/149‑1 (MQ=255)
                                        gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAg                                                                                                       <  7:122727/148‑1 (MQ=255)
                                             gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                                 <  4:194738/149‑1 (MQ=255)
                                                  tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTa                                                                                           <  2:378348/149‑1 (MQ=255)
                                                  tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTGCTACACACATGCCGCAAGATCTGTTGACTGTTGGAGCGGGTTTa                                                                                           >  5:34210/1‑149 (MQ=255)
                                                       ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                                      >  6:285234/1‑149 (MQ=255)
                                                           aTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGaaa                                                                                  <  6:282674/149‑1 (MQ=255)
                                                              aTTTATAAAATTTTTATGAACTGGCTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCGGCGAAACACCTGCCACAAGAAATGGGGACGGTAGGGCGGTGTTTAACCGGCAAAggg                                                                               >  7:59866/1‑149 (MQ=255)
                                                              aTTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTaca                                                                                                                                   <  7:116433/98‑1 (MQ=255)
                                                              aTTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTaca                                                                                                                                   >  8:116433/1‑98 (MQ=255)
                                                                 tataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                                            >  4:88749/1‑149 (MQ=255)
                                                                  ataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                                            >  4:229800/1‑148 (MQ=255)
                                                                   taATATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGt                                                                                                         <  2:153034/119‑1 (MQ=255)
                                                                   taAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGt                                                                                                         >  1:153034/1‑119 (MQ=255)
                                                                         ttttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCCCGCCTAAACACCTGCCTCAAGAACTGTTGACTGTAGGAGAGGGTTTAACCGCGAAAGGGCCGGAACAGGt                                                                    >  4:161208/1‑149 (MQ=255)
                                                                               gAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTc                                                              <  3:88749/149‑1 (MQ=255)
                                                                                  ctaactTTTTGTAATTTTTATATACGGGGTTTTTGTTACAGCGGGGCCGGCCCATTCGCGGGTAAACCCGCTCAAACACACTTGCCTCAACATCTCTCAACTGTAGGTGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                           <  8:28983/149‑1 (MQ=19)
                                                                                  ctaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCg                                                             >  7:301464/1‑147 (MQ=255)
                                                                                       tttttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAACATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCg                                                      <  7:29597/149‑1 (MQ=255)
                                                                                           tgaaaTATTAATACCAGGGGTTTTTGTTACCGCGGTGCCGGGCTCTGGCCGGTCAACCCCGCCCCAACCCAAATCCCCCACGATCTGTGCACTGTAGGCGGGGTTTAACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAg                                                   <  4:226859/145‑1 (MQ=255)
                                                                                                 atTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAg                                            >  5:165621/1‑149 (MQ=255)
                                                                                                     ataGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTATACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAAc                                        >  7:143059/1‑149 (MQ=255)
                                                                                                        gAAGGGGGATTTGTTACAGAGGGGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTTTGAACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACtt                                      <  8:174262/148‑1 (MQ=39)
                                                                                                         aaGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgc                                    <  2:243703/149‑1 (MQ=255)
                                                                                                           gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTTCCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                                   >  3:215353/1‑148 (MQ=255)
                                                                                                              ggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcgc                                  >  5:173977/1‑146 (MQ=255)
                                                                                                                aTTTGTTACCGAGGTGCCGGGCGCTTCGCGGGAAAACCCGCTCCTACACACATGCCTCACCATCTGTGCACTGTCGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGaa                             <  6:369558/149‑1 (MQ=255)
                                                                                                                      tACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTcc                                                >  8:347808/1‑124 (MQ=255)
                                                                                                                      tACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTcc                                                <  7:347808/124‑1 (MQ=255)
                                                                                                                         agagGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGt                    <  4:61395/149‑1 (MQ=255)
                                                                                                                                ccGGGCTCTTCGCGGGTAACCCCGCTCCTACACACAGGCCTCACGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGctct             <  6:322905/149‑1 (MQ=255)
                                                                                                                                 ggggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTctt            <  5:285234/148‑1 (MQ=255)
                                                                                                                                  gggCTCTTCGCGGGAAACCCCGCTCATACACACATGCCTCAAGATCTGTGCACTGAAGGTGCGGTTTAACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttc           <  3:168404/149‑1 (MQ=255)
                                                                                                                                        ttCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcgc                                  >  7:324340/1‑120 (MQ=255)
                                                                                                                                         tCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGgc    >  7:263344/1‑149 (MQ=255)
                                                                                                                                           gcgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGgcgc  >  1:204384/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGC  >  NC_002947/278598‑278884

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: