Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,355 T→C T145T (ACA→ACG 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,3550TC92.6% 65.9 / ‑5.5 28T145T (ACA→ACGytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (12/13);  minor base G (2/0);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGG  >  NC_000913/4624119‑4624593
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
aCGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCCGGCGTAAAATGTTAGTGGGGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGGATGGGTGTCGAGGGTCAGGCGCCGAt                                                                                                                                                                                                                                    >  1:25843/1‑250 (MQ=255)
    tttCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTAACCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCCGAGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGGGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTATGGGGGTCGAGCGTCAGGCGCAGATAGCCg                                                                                                                                                                                                                               >  1:211911/1‑251 (MQ=255)
              ttttgtaATGGGTCCCTCTTTTTTCCCGCTTTCTGCTTCTTCCGGCACTGGTTTGCTGGCGGTTAACAGGAATGGAGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCTCCTGTTGGGCTTGGGTGGCGAGCGTATGGGTGTCTAGCGTCAGGCGCAGATAGCCGAGCGGTcc                                                                                                                                                                                                                       <  2:143663/245‑1 (MQ=255)
              tGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGCTTTCTTCCCGCGGGGACGTTTTCCCTTTAGCAGCGTGCGGGGCAGCACTAAGCCCATTGCCAGTGAAAGCCGCAACAACAGGCGTAAAATGTTAGGGGTGGTTTCCCCCTGTTTGGGTTTGGTGGGGAGCGTATGGGGGGCGAGGGTCAGgcgctgggaggcgggcgggccgg                                                                                                                                                                                                                     >  2:195674/1‑232 (MQ=255)
                 gttaaggtctcctcctTTTGTCTCCCTTTCTTCCTCTTCCTGCACCGGCTTGCTGGCGTTTCACAGGTATGGCTATTGCTTCCAGCGGGTACGTTTTCCCTGTATCAGCGTGCGGGTCAGCGCCACGCAGATGTCCAGTGAACGCAGCATCATCAGGCGTAAAATGTTAGTGGTGTTATCCTCCTGTTGGGCTTTGGTGGCGAGTGTATGGGTGTCGAGCGTCAGGTGCAGATTGCCGTGCGTTCCGtttt                                                                                                                                                                                                                 <  2:214738/237‑1 (MQ=255)
                        ggTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGAAGTTTACCCTGTAGCAGCGGGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACCTCCGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGGATGGGTGTCGGGCGTCAGGCGCAGATAGCCGGGCGGTCCGttttttccgac                                                                                                                                                                                                          >  1:245741/1‑245 (MQ=255)
                               ctgtgttcccGCTTTCTTCCTCTTCCGTCACCGGTTTGCTGGCGGTTAACAGGGAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTTCGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTGCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTg                                                                                                                                                                                                     <  2:120173/240‑1 (MQ=255)
                                    ttCTCGCTTTCTTCCTCTTCCGGCCCCGGTTTGCTGGCGGGTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccgg                                                                                                                                                                                                           <  2:296068/238‑6 (MQ=255)
                                    ttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGGAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGGATGGGGGTCGAGCGTCCGGCGCAGATAGCCGGGCGGGCCGtttttcccggggatcgggagagg                                                                                                                                                                                              >  1:296068/1‑233 (MQ=255)
                                      ccggCTTTCTCCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGAGTGCTGCCAGCGGGTTCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACaa                                                                                                                                                                                            <  2:129525/248‑1 (MQ=255)
                                       tcGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTTCCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGGTTGCCAGTGAAAGCAGCAACATCCGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGGGCGTCAGGCGCAGATAGCCGCGCGGGCCGGTTTTCCCGGCGATTGGTTCGACAAt                                                                                                                                                                                           >  1:1916/1‑251 (MQ=255)
                                              cttcctcttccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGGGCGGGTCCGCACTACGCCGATTGCCAGTGAAAGCAGCCACATCAGGCGTAAAATGTTAGGGGTGTTATCCACCTGTTGGGGTTCGGTGGCGAGCGTATGGGGGACGAGCGTCAGGCGCAGATAACCGAGCGGTCCGGTTTTCCCGGCAAATGGGTCGACAATCTGCTgg                                                                                                                                                                                    >  1:165534/1‑251 (MQ=255)
                                                      ccgggCATCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTTCGTTTACCCGGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTGCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTGAAAATa                                                                                                                                                                            <  2:165534/248‑1 (MQ=255)
                                                             ccGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAgccgccc                                                                                                                                                                     <  1:195674/251‑2 (MQ=255)
                                                                                              gctgcCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTcgcg                                                                                                                                     <  1:26192/250‑1 (MQ=255)
                                                                                                                                                          gTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCcgcg                                                                        >  1:141239/1‑251 (MQ=255)
                                                                                                                                                            gaaagcagcAACATCAGGCGTAAAATGTTAGTGGTGTTATCCCCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTTGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTc                                                                      <  1:13240/248‑1 (MQ=255)
                                                                                                                                                                  agcaACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGAt                                                                <  1:26266/251‑1 (MQ=255)
                                                                                                                                                                   gcaACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGCGGAt                                                                >  1:222245/1‑250 (MQ=255)
                                                                                                                                                                            ggCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTc                                                       >  1:336467/1‑250 (MQ=255)
                                                                                                                                                                                     tGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAAc                                             >  1:324571/1‑251 (MQ=255)
                                                                                                                                                                                      gtgtgtgggtgtctcctccttttGGGCGTCGGTGGCGGGCGTATGGGTGTCGAGCGGTAGGCGCAGATAGCCGAGCGGTCGTGTTTTCCCGGCAATGGGTGCGACAATCTGCTGGTGAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACGTCGACGCTTTCGCCAGAACGGGCGATAAGATCTCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAAc                                             <  2:136968/229‑1 (MQ=255)
                                                                                                                                                                                                 tATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGAGTAAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGGTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCATTCGACGCTTTCGCCAGAACGGGCGATAAGGTCGACTTGTTCGTCATACACGCCCGCCGCGAGGTTACGGATTTCATCCGTTTACTGCACGAGGagc                                 >  2:106966/1‑249 (MQ=255)
                                                                                                                                                                                                      cacgtttGGGCTTCGGGTGCGGGCGTATGGGTGGCGAGCGGCAGGTGCCGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCTCTTTTTTACCGTCGGGCGCCTGCCGGTCGCGCACTTCGACTCTTTCGCCATAACGGGCGATAAGATCGCCTTGGTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTg                            <  2:193796/246‑1 (MQ=255)
                                                                                                                                                                                                       ccTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGa                           <  1:246161/251‑1 (MQ=255)
                                                                                                                                                                                                             gggCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATCAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGt                     >  1:158659/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGttttt                 >  1:307536/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCgg  >  1:271081/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
ACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGG  >  NC_000913/4624119‑4624593

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