Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,339,216 C→T 9.7% intergenic (‑218/+25) PP_2935 ← / ← ybiT membrane protein/ABC‑F family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,339,2160CT 9.7% 108.1 / 3.0 41intergenic (‑218/+25)PP_2935/ybiTmembrane protein/ABC‑F family ATPase
Reads supporting (aligned to +/- strand):  ref base C (24/13);  new base T (2/2);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

CGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTGGCTG  >  NC_002947/3339073‑3339364
                                                                                                                                               |                                                                                                                                                       
cGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCCTGCGAGAt                                                                                                                                                    >  6:94260/1‑149 (MQ=255)
       ttCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCggg                                                                                                                                             >  7:307759/1‑149 (MQ=39)
         cgtcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGc                                                                                                                                            <  8:307759/148‑1 (MQ=39)
            cgGGCCGTGCCGGTCGGTGGGAGGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑c                                                                                                                                     <  1:302210/149‑1 (MQ=255)
               gCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCCTCCGACATGGGCCGGGC‑‑GG‑ccc                                                                                                                                   >  2:119765/1‑148 (MQ=38)
                  gTGCCGGGCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCACAAGCGTCCAATACCTGTGCGGGCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCatca                                                                                                                               <  2:120293/149‑1 (MQ=255)
                          cggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGa                                                                                                                         <  4:229223/147‑1 (MQ=255)
                          cggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGa                                                                                                                         >  3:229223/1‑147 (MQ=255)
                            gtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAg                                                                                                                     <  6:138192/149‑1 (MQ=255)
                                      gatgTACTCTTCGCGGGCTTGCCCGCTCTCCAAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTgg                                                                                                           <  3:229100/149‑1 (MQ=255)
                                               ttCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCa                                                                                                    <  6:144954/147‑1 (MQ=255)
                                                 cgcgGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAgg                                                                                                  >  6:99664/1‑147 (MQ=255)
                                                           cccGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACAAGCACACCCTGGCTAAGCAGGTAGTCATCGTAg                                                                                      >  2:55887/1‑149 (MQ=255)
                                                                    ccAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGTTGGGGTCTGTCCCAGTGAACAAAAATCAGACCAGCCCACCCAGGCTACGCAGGTAGTCATCGTAGCTa                                                                                   >  1:296741/1‑146 (MQ=19)
                                                                    ccAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGt                                                                                                 <  8:207108/129‑1 (MQ=255)
                                                                    ccAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGt                                                                                                 >  7:207108/1‑129 (MQ=255)
                                                                     cAAGGCCCCACAAGCGTCCAATACCTGTGCGCTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCACGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTg                                                                             >  3:62155/1‑148 (MQ=255)
                                                                             cacaAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccac                                                                     >  1:232116/1‑148 (MQ=255)
                                                                                       cAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCgg                                                          >  3:181898/1‑149 (MQ=255)
                                                                                        aaTACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGc                                                         >  1:230634/1‑149 (MQ=255)
                                                                                            ccTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGa                                                      >  3:9755/1‑148 (MQ=255)
                                                                                             cTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGa                                                      >  6:79433/1‑147 (MQ=255)
                                                                                                    gTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCgat                                             >  3:305136/1‑149 (MQ=255)
                                                                                                     tCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCgat                                             >  1:331188/1‑148 (MQ=255)
                                                                                                             ggAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGt                                    <  2:296741/149‑1 (MQ=255)
                                                                                                             gcccggggAAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGt                                    <  1:166177/144‑1 (MQ=255)
                                                                                                                gCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCCACAGCTCGATGATGCGCGTGGc                                 >  6:208748/1‑149 (MQ=255)
                                                                                                                    gCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGc                             >  5:258942/1‑149 (MQ=255)
                                                                                                                      aaGCCCGCGAAGAGGCCAGGCCATCCGAGATGGCCAGGCC‑‑GT‑CACATCAATCAGACCAGCACACACTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCAACAACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGcga                           >  1:66141/1‑149 (MQ=255)
                                                                                                                          ccGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGcgacga                        >  2:233482/1‑148 (MQ=255)
                                                                                                                                    gCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTcgc             >  2:87883/1‑149 (MQ=255)
                                                                                                                                    gCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTcgc             >  5:108890/1‑149 (MQ=255)
                                                                                                                                      ctggCCATCTCGGATGGCCGGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGacac                     <  7:115424/137‑1 (MQ=255)
                                                                                                                                      ctggCCATCTCGGATGGCCGGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGacac                     >  8:115424/3‑139 (MQ=255)
                                                                                                                                      ccGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACg                    >  8:139964/1‑140 (MQ=255)
                                                                                                                                      ccGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACg                    <  7:139964/140‑1 (MQ=255)
                                                                                                                                           cATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTg      <  2:232116/149‑1 (MQ=255)
                                                                                                                                            aTCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTgg     <  6:108890/149‑1 (MQ=255)
                                                                                                                                             tCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTc                                                                        <  4:100667/81‑1 (MQ=255)
                                                                                                                                             tCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTc                                                                        >  3:100667/1‑81 (MQ=255)
                                                                                                                                             tCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTGGc    >  2:257899/1‑149 (MQ=255)
                                                                                                                                               cGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTGGCt   >  8:35032/1‑148 (MQ=255)
                                                                                                                                               cGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTGGCTg  <  4:9755/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
CGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCC‑‑GT‑CCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTGGCTG  >  NC_002947/3339073‑3339364

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