Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,425,329 Δ1 bp intergenic (+90/‑48) trkG → / → ynaK K(+) transporter TrkG/ParB‑like nuclease domain‑containing protein YnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,425,3290C.93.5% 119.8 / 0.2 31intergenic (+90/‑48)trkG/ynaKK(+) transporter TrkG/ParB‑like nuclease domain‑containing protein YnaK
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (16/13);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATGGGATTAGGTTTTGGGGCTACTGCATCGTCATTCGGAGTGCTTAATGACATTGCAAAATGCTTAATGTGCATAGCTATGATTCTTGGTCGCCTTGAAATTTATCCTGTTATTATATTGTTTTCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGA  >  NC_000913/1425091‑1425570
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
atatGGGATTAGGTTTTGGGGCTACTGCATCGTCATTCGGAGTGCTTAATGACATTGCAAAATGCTTAATGTGCATAGCTATGATTCTTGGTCGCCTTGAAATTTATCCTGTTATTATATTGTTTTCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCAGTGgct                                                                                                                                                                                                                                      <  1:72477/251‑1 (MQ=255)
                                    tCGGAGTGCTTAATGACATTGCAAAATGCTTAATGTGCATAGCTATGATTCTTGGTCGCCTTGAAATTTATCCTGTTATTATATTGTTTTCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTa                                                                                                                                                                                                   <  1:179385/250‑1 (MQ=255)
                                     cGGAGTGCTTAATGACATTGCAAAATGCTTAATGTGCATAGCTATGATTCTTGGTCGCCTTGAAATTTATCCTGTTATTATATTGTTTTCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATg                                                                                                                                                                                                 >  1:125353/1‑251 (MQ=255)
                                                   aCATTGCAAAATGCTTAATGTGCATAGCTATGATTCTTGGTCGCCTTGAAATTTATCCTGTTATTATATTGTTTTCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATACAGCAAGGCTGTATGAGCGAGAAATTaa                                                                                                                                                                                    >  1:139514/1‑250 (MQ=255)
                                                        gCAAAATGCCTAATGTGCATAGTTATGAGTCTTGGTCGCCTTGAAAATTATCCTGTTATTATATTGTTTTCTGGTTTTTTTTGGCTCTTCTTATATATGTCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTTCCATCAACACCTTCGTCCCAGTAGCACTGGCCGCTTGGGTAATTTTCATTCATTAAGCGAGCCTGTAGGAGCGAGAAATTACAGATAg                                                                                                                                                                              <  2:346247/251‑1 (MQ=255)
                                                                 ttAATGTGAATAGCTATGATTCTTGGACGCATTGAAATTTATCCTGTTATTCTATTGTTATCAGGTTTTTTTTGGCGCTCCTAATATATGACAGATTTATAATTGTGAGTTTAAAATTATATGGACTCACTCATTGATCCATAACCTAAGTTTACAAGAAACCCATCCGTCCAAGTAGACATGGACGCTGGGGTGCGTTTTATTAATAATGACAGGATGTCTGAGAGGGAAATTAACGATAGTCTATctcc                                                                                                                                                                     >  2:210876/1‑248 (MQ=255)
                                                                        gCATAGCTATGATTCTTGGTCGCCTTGAAATTTATCCTGTTATTATATTGTTTTCAGGTTTTTTTTGTCGCTCCTAATATATTGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACa                                                                                                                                                              <  1:108708/251‑1 (MQ=255)
                                                                                 tGATTCTTGGTCGCCTTGAAATTTATCCTGTTATTATATTGTTTTCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAACACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTTTATCGCCCATTACAAGAATTGTc                                                                                                                                                     >  2:113300/1‑251 (MQ=255)
                                                                                                    aaTTTATCCTGTTATTATATTGTTTTCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGACACCGTATGCGCAAAACGcc                                                                                                                                  >  2:272060/1‑251 (MQ=255)
                                                                                                       ttATCCTGTTATTATATTGTTTTCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAACACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAg                                                                                                                                <  1:113300/250‑1 (MQ=255)
                                                                                                          tCCTGTTATTATATTGTTTTCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACg                                                                                                                            >  1:340702/1‑251 (MQ=255)
                                                                                                                             tCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTg                                                                                                          <  1:129488/250‑1 (MQ=255)
                                                                                                                              cAGGTTTTTTTTGGCGCACCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCAGTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGc                                                                                                        <  2:164883/251‑1 (MQ=255)
                                                                                                                                       tttGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTgg                                                                                                <  1:272060/250‑1 (MQ=255)
                                                                                                                                               cccttagatagggcataattATAATTGTGAATTTAATAGTATGTTGATCCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCACAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCTCCCATTACAAGAATTGTCACCGTATGCGCACAACACCAGGACGCACAGTACTTAGCAGGTGGCACAACTGGTAGAAAGTa                                                                                       <  2:5437/233‑1 (MQ=255)
                                                                                                                                                  tAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGc                                                                                                               >  2:83176/1‑224 (MQ=255)
                                                                                                                                                      atatTGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATTCGCACAACGCCATGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGc                                                                                 <  2:143115/250‑1 (MQ=255)
                                                                                                                                                        atGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAAt                                                                              >  1:66909/1‑251 (MQ=255)
                                                                                                                                                              gATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCt                                                                        <  1:91015/251‑1 (MQ=255)
                                                                                                                                                                              tttAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCt                                                        >  1:311776/1‑251 (MQ=255)
                                                                                                                                                                                tAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCt                                                        >  1:141996/1‑249 (MQ=255)
                                                                                                                                                                                tAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCGCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCt                                                        <  2:141996/249‑1 (MQ=255)
                                                                                                                                                                                    attatGTTTACCCACTCATTGTTCCAATACATACCTTTACCAGCAACTCAGCCGGCCCAGTAGCTCTGGCTGCTTGGGTGAGGTATATGCTGAAAGTAAGGCGGTATGAGCGTGAAATTAAAGAGAGGCTCTCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGTCAGGTCGCACAGTACTGAGCATGTGTCGCAACGGGTAGATAGTTATAAGCAATTCGGCTTGTTTCATCCGGTGCTGAtag                                                   <  2:310337/250‑3 (MQ=255)
                                                                                                                                                                                    attatGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCTGATTg                                                   >  1:286246/1‑250 (MQ=255)
                                                                                                                                                                                                cgctcaTTGATCCGATACCTAACGTTGCCAGCATCACCTCCGGCTCAGTAGCACGGGATGCGGGGGTGCGTTTTAGTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCGCAAAGCCAGGACGCACAGTACTGAGCAGGTGGTACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGa                                      <  2:125869/249‑1 (MQ=255)
                                                                                                                                                                                                  ctcaTTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACAGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGACATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGCGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGCCTCATCCGGTGCTGATTGACGGAAAGGACGa                                      >  1:125869/1‑249 (MQ=255)
                                                                                                                                                                                                             aaTACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCACATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGGAATTATTGCGGGt                         >  2:196227/1‑251 (MQ=255)
                                                                                                                                                                                                              aTACCTAACTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATACAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTc                        >  2:314001/1‑251 (MQ=255)
                                                                                                                                                                                                                    aaCTTTACCAGCAACACCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCAGAAAGCAAGGCGGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGGGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGCGCTGATTGACGAAAAGGGCGAAATTATTGCGGGGCACGGTc                  <  1:315059/251‑1 (MQ=255)
                                                                                                                                                                                                                                   acCTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGGGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCAACGTCTGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCTGATTTACGAAAAGGGAGAAATTATTGCGGGTCACGGTCGTGTTTTggcggggg   >  2:222449/1‑248 (MQ=255)
                                                                                                                                                                                                                                     cTCCGTCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGa  >  1:101955/1‑250 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
ATATGGGATTAGGTTTTGGGGCTACTGCATCGTCATTCGGAGTGCTTAATGACATTGCAAAATGCTTAATGTGCATAGCTATGATTCTTGGTCGCCTTGAAATTTATCCTGTTATTATATTGTTTTCAGGTTTTTTTTGGCGCTCCTAATATATGGCTGATTTATAATTGTGAGTTTAATATTATGTTGACTCACTCATTGATCCAATACCTAACTTTACCAGCAACACCTCCGTCCCCAGTAGCACTGGCTGCTGGGGTGCGTTTTATTCATAAAGCAAGGCTGTATGAGCGAGAAATTAAAGATAGTCTATCGCCCATTACAAGAATTGTCACCGTATGCGCACAACGCCAGGACGCACAGTACTGAGCAGGTGGCACAACTGGTAGAAAGTATTAAGCAATTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGA  >  NC_000913/1425091‑1425570

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: