Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,139,356 C→T 7.5% intergenic (+145/‑91) PP_2754 → / → sfnB OprD family outer membrane porin/SfnB family sulfur acquisition oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,139,3560CT 7.5% 150.7 / 2.4 53intergenic (+145/‑91)PP_2754/sfnBOprD family outer membrane porin/SfnB family sulfur acquisition oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (18/31);  new base T (2/2);  total (20/33)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

GAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATCACGCTTGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAG  >  NC_002947/3139210‑3139494
                                                                                                                                                  |                                                                                                                                           
gAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGa                                                                                                                                           >  1:116706/1‑149 (MQ=255)
         cTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGaa                                                                                                                                  <  7:72486/149‑1 (MQ=255)
            aCGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTg                                                                                                                               <  6:116676/149‑1 (MQ=255)
              ggTCTCGCGAGATACCGCAATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCa                                                                                                                             <  4:207954/149‑1 (MQ=255)
               gTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCAc                                                                                                                            <  6:221569/149‑1 (MQ=255)
               gTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCAc                                                                                                                            <  8:88726/149‑1 (MQ=255)
               gTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCAc                                                                                                                            <  6:257946/149‑1 (MQ=255)
                    gcgAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCaa                                                                                                                       <  7:51436/149‑1 (MQ=255)
                    gcgAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCa                                                                                                                        <  2:116706/148‑1 (MQ=255)
                          tACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCAtt                                                                                                                 <  6:137749/149‑1 (MQ=255)
                          tACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCAtt                                                                                                                 <  7:99878/149‑1 (MQ=255)
                          tACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCAtt                                                                                                                 <  8:29842/149‑1 (MQ=255)
                             cGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGaa                                                                                                              >  3:23595/1‑149 (MQ=255)
                                     gtgttgcaTGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAgg                                                                                                      <  2:154472/142‑1 (MQ=255)
                                           gATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTcc                                                                                                <  6:132594/149‑1 (MQ=255)
                                                  gcCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTg                                                                                         >  5:233119/1‑149 (MQ=255)
                                                     tGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTa                                                                                      >  3:104334/1‑149 (MQ=255)
                                                      ggCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCTCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTAc                                                                                     >  7:199896/1‑149 (MQ=255)
                                                                 gCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCAc                                                                         <  8:19827/149‑1 (MQ=255)
                                                                 gCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCAc                                                                         <  8:176769/149‑1 (MQ=255)
                                                                 gCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCAc                                                                         >  5:33230/1‑149 (MQ=255)
                                                                  cGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCAc                                                                         >  1:240296/1‑148 (MQ=255)
                                                                   gaggagCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAAtt                                                                                                                                <  7:135161/93‑1 (MQ=255)
                                                                   gaggagCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAAtt                                                                                                                                >  8:135161/1‑93 (MQ=255)
                                                                        gCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCTCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCt                                                                  <  8:199896/149‑1 (MQ=255)
                                                                          ggTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTgg                                                                <  7:126914/149‑1 (MQ=255)
                                                                          ggTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTgg                                                                <  7:30794/149‑1 (MQ=255)
                                                                          ggTAGGCATGTTCTGAGGCAAGTGTACGCACTGCAACAGCAAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCTAGTGCAC‑AGcccgta                                                                <  7:90641/149‑4 (MQ=255)
                                                                            tCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTggg                                                               >  8:159847/1‑148 (MQ=255)
                                                                              ggCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGaa                                                            <  4:23595/149‑1 (MQ=255)
                                                                              ggCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGaa                                                            >  8:173224/1‑149 (MQ=255)
                                                                                   gTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTa                                                                       >  8:232456/1‑133 (MQ=255)
                                                                                   gTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTa                                                                       <  7:232456/133‑1 (MQ=255)
                                                                                    ttCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGTCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGt                                                                                                                          <  1:184288/82‑1 (MQ=255)
                                                                                    ttCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGt                                                                                                                          >  2:184288/1‑82 (MQ=255)
                                                                                     tCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATaagaa                                                                                                        >  5:52789/1‑99 (MQ=255)
                                                                                     tCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATaagaa                                                                                                        <  6:52789/99‑1 (MQ=255)
                                                                                               gtgtCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGtt                                                                                          <  6:248733/103‑1 (MQ=255)
                                                                                                  tCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATCACGCTTGGAGCaa                                         <  6:33230/148‑1 (MQ=255)
                                                                                                      cccGGCAACAGACAATCCACAGTTTGCCGACACAGGTGATTTTCCGATTATCTGAATTGCACGTCAAGGCATTGAATAAGAAGGATTTCCGGAGTTGGTACGGGCT‑GTTCCCTATCCCTGGGGAATCCATCACGCTTGGAGCAAAAGAt                                    >  2:101291/1‑149 (MQ=255)
                                                                                                                 ccAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATCACGCTTGGAGCAACAGATGACCCAGTca                          >  4:127906/1‑148 (MQ=255)
                                                                                                                     tCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGAGTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATCACGCTTGGAGCAACAGATGACCCAGTCACAAcc                     <  7:173224/149‑1 (MQ=255)
                                                                                                                       cacaGTTTGACGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATCACGCTTGGAGCAACAGATGACCCAGTCACAACccg                   >  4:103157/1‑149 (MQ=255)
                                                                                                                        acaGTTTGCCGCCACAGGTGCTTTTCTGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTggg                                                               >  2:147555/1‑104 (MQ=255)
                                                                                                                        acaGTTTGCCGCCACAGGTGCTTTTCTGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTggg                                                               >  4:11644/1‑104 (MQ=255)
                                                                                                                        acaGTTTGCCGCCACAGGTGCTTTTCTGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTggg                                                               <  3:11644/104‑1 (MQ=255)
                                                                                                                        acaGTTTGCCGCCACAGGTGCTTTTCTGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTggg                                                               <  1:147555/104‑1 (MQ=255)
                                                                                                                        acaGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATCACGCTTGGAGCAACAGATGACCCAGTCACAACccgc                  <  4:223397/149‑1 (MQ=255)
                                                                                                                         caGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATCACGCTTGGAGCAACAGATGACCCAGTCACAACccgcc                 <  6:233119/149‑1 (MQ=255)
                                                                                                                               gccgccACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATc                                                      <  1:45385/106‑1 (MQ=255)
                                                                                                                               gccgccACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATc                                                      >  2:45385/1‑106 (MQ=255)
                                                                                                                                         gTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATCACGCTTGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAg  >  7:17040/1‑148 (MQ=255)
                                                                                                                                         gTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATCACGCTTGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAg  <  8:17040/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GAGCACCAGCTCACGGTCTCGCGAGATACCGCCATGCAACACCGATGGCGGCCTGGCAAAAACAGGCGAGGAGCGGTCGGCAGGTTCTGAGGCAAGTGTCCGCCCGGCAACAGCCAATCCACAGTTTGCCGCCACAGGTGCTTTTCCGATTATCTGAATTGCACGTCAAGTCATTGAATAAGAAGGATTTCCTGAGTTGGTACGGGCT‑GTGCACTAGCCCTGGGGAATCCATCACGCTTGGAGCAACAGATGACCCAGTCACAACCCGCCGATCTGCGTGAGCAG  >  NC_002947/3139210‑3139494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑