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.A90.7% 160.5 / 10.6 54intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base A (19/30);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCT  >  NC_002947/278595‑278885
                                                                                                                                                    |                                                                                                                                                 
aTGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGt                                                                                                                                                    <  3:72597/148‑1 (MQ=255)
     gCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATCATTTATAAAATTTTTATCAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAaccc                                                                                                                                           >  7:240067/1‑148 (MQ=255)
      cTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAacccg                                                                                                                                          >  1:273885/1‑147 (MQ=255)
      cTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAacccg                                                                                                                                          >  3:164184/1‑147 (MQ=255)
         cAGCTGATCGCTACAGACGTGTGCTGGTCCGGTGTTCAAGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCt                                                                                                                                         <  1:124511/148‑1 (MQ=255)
         cAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTc                                                                                                                                        <  3:22769/149‑1 (MQ=255)
            cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                     >  6:71822/1‑149 (MQ=255)
                    gACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATg                                                                                                                             <  5:178525/149‑1 (MQ=255)
                    gACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATg                                                                                                                             <  4:247867/149‑1 (MQ=255)
                          cGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                                       <  4:38648/149‑1 (MQ=255)
                             gtgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGa                                                                                                                     >  8:296624/1‑148 (MQ=255)
                                cTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtg                                                                                                                 >  1:156009/1‑149 (MQ=255)
                                 tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                                >  8:141709/1‑149 (MQ=255)
                                 tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                                >  7:74042/1‑149 (MQ=255)
                                     ccgtttgtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAAGTTTTTTGAATATTTATAGAATGGTGATTTGTTACAGAGGTGCCGGGCTCTTCGCGTGTAAACCCGCTCCTACACACATTCCTCAAGATCTGTGCa                                                                                                             <  1:364730/142‑1 (MQ=255)
                                       ggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGAt                                                                                                                    <  7:207907/139‑1 (MQ=255)
                                       ggtggtccggtATAAATACACTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGAt                                                                                                                    >  8:207907/1‑139 (MQ=255)
                                                 tatatataCCATAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                                 <  7:296624/148‑1 (MQ=255)
                                                       taCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAc                                                                                           <  5:71822/148‑1 (MQ=255)
                                                         cccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTc                                                                                                                         <  6:29545/116‑1 (MQ=255)
                                                         cccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTc                                                                                                                         >  5:29545/1‑116 (MQ=255)
                                                         cccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAccc                                                                                         <  1:272226/148‑1 (MQ=255)
                                                         cccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcg                                                                                        <  2:156009/149‑1 (MQ=255)
                                                         cccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcg                                                                                        <  5:152568/149‑1 (MQ=255)
                                                          ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                                       <  7:141709/149‑1 (MQ=255)
                                                            tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACCCATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgcg                                                                                      >  2:264709/1‑148 (MQ=255)
                                                            tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGa                                                                                     <  6:321594/149‑1 (MQ=255)
                                                            tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAAACGCTCCTTAACACATGGCGCAAGATATGTGGACTGTAGGAGCGGGTTTTACcgcg                                                                                      >  3:204776/1‑148 (MQ=255)
                                                               tGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTATACCCCCTACTACACACATTCCTCAAGATCTGTTCACTGTAGGAGCGGGTTTAAACGCGaaat                                                                                  >  6:326647/1‑148 (MQ=255)
                                                                 aTTTAAAATATTTTTATTAACAAACTTTTTGGAATATTTATCGAAGGGGGATTTGTTACAGCGGGGCCGGGCTCTTCGGGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTCCACTGTAGGAGCGGGTTTACCCGCGAAAgg                                                                                 <  2:299439/148‑1 (MQ=255)
                                                                   ttATAATTTTTTTATTAACTATCTTTTTTTAATTTTTATAGACGGGGTATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAACATCTGTTCACTGTCGGAGCGGGTTAACCCGCGAAAGGGc                                                                               <  2:273885/148‑1 (MQ=255)
                                                                      taAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCa                                                                                                                        <  4:121376/104‑1 (MQ=255)
                                                                          aTTTTAATGAACTAACTTTTTGGAAAATTTAAAGACGGGGGATTTGTTCCAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCAACACACATCCCTCAAGATCTGTCCACTGTCGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAg                                                                       <  6:207395/149‑1 (MQ=255)
                                                                           tttttATGAACTAACTTTTTGTAATATTTATAGAAGGGGCATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAgg                                                                      <  1:256976/149‑1 (MQ=255)
                                                                              ttttGAACTAATTTTTTGTAATATTTATAGATGGGGGATTTGTTACAGCGGGGCCGGGCTCTTCGGGGTTAAACCCGCTCCTACACACTGGCCTCAATATCTGTGCACTGTCGGCGGGGGTTTACCCGCGAAAGGGCCGGAACAGGtt                                                                    <  4:164184/145‑1 (MQ=255)
                                                                                   aactaactTTTTGGAATATTTATAGAAGGGGGATTTGTTAGAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCg                                                              >  6:69338/1‑149 (MQ=255)
                                                                                      taactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGTGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAg                                                           >  1:81578/1‑149 (MQ=255)
                                                                                                     tttATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAg                                             >  3:275977/1‑148 (MQ=255)
                                                                                                           gAAGGGGTATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTg                                      <  2:115364/149‑1 (MQ=255)
                                                                                                             aGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCGCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                                    >  2:23327/1‑149 (MQ=255)
                                                                                                              gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGAt                                                                                                                    >  7:162133/1‑68 (MQ=255)
                                                                                                              gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGAt                                                                                                                    <  8:162133/68‑1 (MQ=255)
                                                                                                                   aTTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGaaa                                                                                   <  6:196636/96‑1 (MQ=255)
                                                                                                                        ttACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGcc                         >  6:357760/1‑149 (MQ=255)
                                                                                                                          aCAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTAc                                                                 <  1:278199/107‑1 (MQ=255)
                                                                                                                          aCAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTAc                                                                 >  2:278199/1‑107 (MQ=255)
                                                                                                                           cAGAGGTGCCGGGCTCTGCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgc                       <  8:213880/148‑1 (MQ=255)
                                                                                                                                gggcggggCTGTTCGCGGTTAACCCCGAACCTACACACTTGCCTAAATATCTGTCCACTGTAGGCGCGGGTTAACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACg                  <  1:169125/143‑1 (MQ=37)
                                                                                                                                  gCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                            <  7:188088/104‑1 (MQ=255)
                                                                                                                                  gCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                            >  8:188088/1‑104 (MQ=255)
                                                                                                                                          ctTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCg       >  1:343241/1‑149 (MQ=255)
                                                                                                                                           ttCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                            >  4:29856/1‑95 (MQ=255)
                                                                                                                                           ttCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                            <  3:29856/95‑1 (MQ=255)
                                                                                                                                               cgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCt  <  1:23327/149‑1 (MQ=255)
                                                                                                                                               cgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCt  <  4:275977/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCT  >  NC_002947/278595‑278885

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: