Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,586 T→C S68S (TCA→TCG ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,5860TC91.2% 73.6 / ‑1.2 34S68S (TCA→TCGytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (16/15);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

ATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCA  >  NC_000913/4624343‑4624826
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
aTGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCgg                                                                                                                                                                                                                                            >  1:271081/1‑251 (MQ=255)
             gtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCgg                                                                                                                                                                                                                               >  1:150564/1‑251 (MQ=255)
                             cGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGCGGCTACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGACTGCGTTTTTCATCCGGCGAGCCGGTACGCATCAGCGGTGCAACGTTCAGCGtt                                                                                                                                                                                                               >  1:247334/1‑250 (MQ=255)
                                   gTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTgg                                                                                                                                                                                                         >  1:259552/1‑250 (MQ=255)
                                   gTCCGTTTTTCCCGGCAATTGGTTCGACAACCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATCCACGCCCGCGTCGAGGATACGGCTTTCATCCGTTATCTGATCGAGGATCGCCTGAATGCGTTGTTCATCCGTCGAGTCGGTACGCATCAGCGGTGCAACTTTCAGCGTTACCTGa                                                                                                                                                                                                         <  1:331406/251‑1 (MQ=255)
                                             cccggCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAgg                                                                                                                                                                                                <  1:130340/246‑1 (MQ=255)
                                                    aTTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGc                                                                                                                                                                                        >  1:96470/1‑251 (MQ=255)
                                                          tCGACAATCTGCGGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGtt                                                                                                                                                                                   <  1:58309/250‑1 (MQ=255)
                                                            gACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCAAGCCGGTCGCGCACTTCGACGCTTTCGCACGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGAAACGGCTTTCATACGTTAACTGCTCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGCTCAGCGTTACCTGACGCGCCAGGGTGAGGGCCAGttctt                                                                                                                                                                                >  1:182916/1‑251 (MQ=255)
                                                                          tAAAATAGCCGCCCGCTTTTTTACCGTCTAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGCGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGGGCAAAGTTCAGCGTTACCTGACGCGCCAGGGGGCGGGCCCGTTCTTCCGGCTGTgggta                                                                                                                                                                   >  2:284025/1‑247 (MQ=255)
                                                                           aaaaTAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTGTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTa                                                                                                                                                                  <  1:284025/250‑1 (MQ=255)
                                                                               tAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGCTAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACgttg                                                                                                                                                             >  1:279220/1‑248 (MQ=255)
                                                                                 gcctcccgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGACGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTc                                                                                                                                                           <  1:80294/244‑1 (MQ=255)
                                                                                           tttttACCGTCGGGCCCCAGCCGGTCGCGCACTTCGACGCTTTCGCCCGAACGTGCGATAAGACCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTggtggt                                                                                                                                                  <  2:96470/250‑1 (MQ=255)
                                                                                                     cGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAAc                                                                                                                                       <  1:295830/251‑1 (MQ=255)
                                                                                                         cgcCAGCCGGTCGCGAACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTCCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCAAGCTGTGGATTACGTTGTCGCTGGTGGTTTTAACTAAACAATg                                                                                                                                  >  1:187333/1‑251 (MQ=255)
                                                                                                                     gcgcACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGGTCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGc                                                                                                                       >  1:43477/1‑250 (MQ=255)
                                                                                                                                    ttCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAAc                                                                                                        <  1:64123/250‑1 (MQ=255)
                                                                                                                                       gCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGCGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGAGCCAGGGAGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCAGCAACGCCCCTAACaggg                                                                                                    >  1:55798/1‑248 (MQ=255)
                                                                                                                                             acgggCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGAGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACACGGCAAGAc                                                                                              >  1:58416/5‑251 (MQ=255)
                                                                                                                                              cgggCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGCTCGAGGATCGCCTGAATGCGTTTTTCATACGGCGAGTCGGTACGCATCAGCGGGGCAACGTTCAGCGTTACCTGACGCCCCAGGGGGCGGGCCCGTTCATCCAGCTGTGGATGACGTTTTCGCTGGTGGTTTTGCATAACCCATTGACCTCCCAGCATCCGCGCCCCAACCAAGGCCAGAc                                                                                              >  2:195673/4‑250 (MQ=255)
                                                                                                                                              cgggCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACa                                                                                             <  1:332829/248‑1 (MQ=255)
                                                                                                                                              cgggCGATAAGATAGACTTGTTCGTCCTACACGCCCGCGACGAGGATACGGCTTTCATACGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCCGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCAAGCTGTGGATTACGTTGTAGCTGGTGGTTTTGCCTAAACCATGACGCTCCCTGCATAAGCGACACTACTAAGGCAAgcc                                                                                              >  2:321267/4‑248 (MQ=255)
                                                                                                                                                    agaaGATCGCTTTTTTCGTCATACACGCCCGCGGCTATGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCACCCGGCGAGTCGGTACGCATCAGCGTTGCAACTTTCAGCGTTACCTGACGCGTCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAg                                                                                       <  2:324571/249‑1 (MQ=255)
                                                                                                                                                                   tCGTCATACTCTCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACTTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCGGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACACGGCAATACAGAACAGGACAATCTCTGccc                                                                         <  2:307536/250‑1 (MQ=255)
                                                                                                                                                                    cGTCATACACACCCGCGTCGAGGATACGTCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCTTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACCAACAAGGCAAGACAGAACAGGACGATCACTGCCCgg                                                                       >  1:41237/1‑250 (MQ=255)
                                                                                                                                                                                                   ttCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGcc                                         >  1:258307/1‑250 (MQ=255)
                                                                                                                                                                                                        ccGTTAACTGAGCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGTTTCAACGTTCAGCGTTACCTGACGCGGCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAGCAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCtt                                    <  1:31330/250‑1 (MQ=255)
                                                                                                                                                                                                               gcggtgaaGGATCGCTTGAATGCGTTTTACATCCGGCGATTCGTGAGTTGTCAGCGGTTCCACGTTGAGCTTTACCTGACGGGCCTGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTGTTGACTAACCCATGACGCTCCCTGCAGCAGCGCCACTAAAAAGGCAAGTCAGAACAGGACAATCACTGCCCGATGCATCCGGAATGTCAGTTTTGTGCGAGCTATCTTCCAACCtt                            <  2:245741/244‑1 (MQ=255)
                                                                                                                                                                                                                 gATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTg                          <  1:211432/251‑1 (MQ=255)
                                                                                                                                                                                                                       ggATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAAtt                    <  1:142685/251‑1 (MQ=255)
                                                                                                                                                                                                                               tGAATGCGTTTTTCATCCGGCGAGTCGGTACGCCTCAGCGGTGCAACGTTCAGCGTTACCTGAAGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGAGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACtt            >  1:324834/1‑251 (MQ=255)
                                                                                                                                                                                                                                gAATGCGTTTTTAATCCAACAAGTCGGTACGCATCAACGGTGCAACGTGCAGCGTTACCTGGCGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTACACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACCGGACAATAACTGCCCGCTGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTgagaattt           >  2:257117/1‑247 (MQ=255)
                                                                                                                                                                                                                                          ttCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCa  <  1:195673/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
ATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCA  >  NC_000913/4624343‑4624826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: