Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,227,339 C→T T369I (ACC→ATC)  ygjI → putative transporter YgjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,227,3390CT89.7% 77.0 / ‑5.6 29T369I (ACC→ATC) ygjIputative transporter YgjI
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (16/10);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 7.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAAT‑TTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTGC  >  NC_000913/3227095‑3227582
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
ggcggtcgtgtgtCCTTTACCGCGATTGTCGGTTCCCTGCTGATGTGGACCGCAACGCCG‑TGGAAATCTTTCTTCTCCGAAATCCCGGAAGGCAACTTTGGTAAGAAAACCGTCGAACGGAACGAACACGGCGTTCCGGCGCGCGCATCGGGGATCCAGTTCCTGAGCGGCATCCCGCTTATGATTATCCGGATGCTCGGTTCCAATACCTTGCAGGATCTGATAAATACTATTAATAATATGAtcgcct                                                                                                                                                                                                                                                <  2:175737/241‑2 (MQ=255)
    gTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAAT‑TTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCg                                                                                                                                                                                                                                           >  1:116182/1‑251 (MQ=255)
                            tCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAAT‑TTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATtcat                                                                                                                                                                                                                    <  1:14136/250‑1 (MQ=255)
                                  cccTGCTGATGTGGACCGCAACGCCGGTGAAAAT‑TTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTg                                                                                                                                                                                                             <  2:116182/251‑1 (MQ=255)
                                        tgATGTGGACCGCAACGCCGGTGAAAAT‑TTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTAc                                                                                                                                                                                                       >  1:52611/1‑251 (MQ=255)
                                             tgGACCGCAACGCCGGTGAAAAT‑TTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGa                                                                                                                                                                                                   >  2:87875/1‑250 (MQ=255)
                                                    cAACGCCGGTGAAAAT‑TTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTAcgc                                                                                                                                                                                           >  1:137628/1‑251 (MQ=255)
                                                               aaaaT‑TTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAAtt                                                                                                                                                                                   >  1:239281/1‑248 (MQ=255)
                                                                              cGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACGGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATCCCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCAcgcg                                                                                                                                                                  <  2:13328/251‑1 (MQ=255)
                                                                                  aTCCCGGAAGGCCTCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGAttt                                                                                                                                                              >  1:105412/1‑251 (MQ=255)
                                                                                    cccGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTcc                                                                                                                                                            <  1:150994/251‑1 (MQ=255)
                                                                                                  ttGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTccc                                                                                                                                                <  1:375277/249‑1 (MQ=255)
                                                                                                        aGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCAcc                                                                                                                                        >  2:218000/1‑251 (MQ=255)
                                                                                                               cGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGtatt                                                                                                                                  >  2:377675/1‑250 (MQ=255)
                                                                                                                    aacaaaaCGAAACCGGCGTCCCGGCGCGCGCAGCGTGGATCAAGTCCCGAATCGTCACCCCGGTGATGATTACCCCGATGCTGGGTTCAAATCCCGTGCAGGATCTGATGAATATTATTATAAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTCCCTGAATTTCCGCGCCAAATTAGATCCCCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCgtt                                                                                                                            <  2:105412/246‑1 (MQ=255)
                                                                                                                       tGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCaaa                                                                                                                                                                                     >  1:61924/1‑191 (MQ=255)
                                                                                                                       tGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCaaa                                                                                                                                                                                     <  2:61924/191‑1 (MQ=255)
                                                                                                                       agaaagAAAGCGGCGTTCCGGCGCGCGCAGGGGGGATCCAGTCCCTGATCGTCCCCCCGCGGATGTTTACCCCGATGCCCGTTTCAAATCCCTGGCGGTTTCTGATGAATCTTATTATAAATTTGATCGCCGCGGCGTCCATGCTCCCGCGTTTATTCATCATGCTGCCTTACCTGAATTTCCGCGCCAAATTAGTTCACCTGCCACTCGATTTCCTTGTGGGCTCCCGCCGCACCGGTATTATCgttgtt                                                                                                                         <  2:180445/246‑1 (MQ=255)
                                                                                                                             aaaaCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCACCCCGCTGATGATTACCCCGATGCTCGGTTCAAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTCCCGCCGTTATTCATCATGCTGGCTTCCCTGAATTTACGCGCCAAATTAGACCCCCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATg                                                                                                                   <  1:87875/251‑1 (MQ=255)
                                                                                                                              aaaCGGCGTTCCGGCGCGCGCAGCGTGGACCCAGTCCCTGATCGTCATCCCGCTGATGATTACCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCCCGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGc                                                                                                                  <  1:316645/251‑1 (MQ=255)
                                                                                                                                        ccGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGat                                                                                                        >  1:82818/1‑251 (MQ=255)
                                                                                                                                                      gTGGATCCCGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGCTATTTGCCGTAGGGt                                                                                          >  2:54404/1‑251 (MQ=255)
                                                                                                                                                            ccAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCg                                                                                    >  1:131399/1‑251 (MQ=255)
                                                                                                                                                              aGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCt                                                                                  >  1:423751/1‑251 (MQ=255)
                                                                                                                                                                          tCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATCCCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCAAAATTAGATCACCTGCCCCGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTccc                                                                       <  1:377675/250‑1 (MQ=255)
                                                                                                                                                                                                                       gggatctaataaatttttttataaATAGGACCGCCCCAGCGACACTGCTCCCGCCGTTATTCATCAGGCTGGCTTACCTGAATTTCCGCACCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGt                         <  2:304774/233‑1 (MQ=255)
                                                                                                                                                                                                                            cTGATGAATAATATTATTAATATGATCGCCGCAGCGTCCATGCTTCAGCCGTTATTCATCATGATGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGCCTGGCGCGAATATTCTGACCATCATTTTTTATAACGTAGGCGGTATTGt                    >  2:188338/1‑251 (MQ=255)
                                                                                                                                                                                                                              gtttacgactataataattatGACCGCCGCAGCCCCCTTGTCTCCGCATTTATTCATCATGCTGTCTTACCTTAATTTACGCGCCAAATTATACCACCTCCCCTGTGTTTTCCGTATGTGCCCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGtt                   <  1:238857/232‑1 (MQ=255)
                                                                                                                                                                                                                                tgaATACTATTTTTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTTTAACGTCGGCGGTATTGttttt                >  2:187907/1‑248 (MQ=255)
                                                                                                                                                                                                                                              aaTATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTGc  >  1:228642/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
GTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAAT‑TTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTGC  >  NC_000913/3227095‑3227582

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