Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,510,609 C→T P327L (CCG→CTG)  yfeO → putative transport protein YfeO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,510,6090CT95.5% 59.1 / NA 22P327L (CCG→CTG) yfeOputative transport protein YfeO
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (12/9);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGG  >  NC_000913/2510366‑2510809
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                         
ggTTGCACGAGATGATGCATCAAATGAAAAATACGGTGCTCGTGCTGGGTATTGGAGGATTTATTCTCGGTAATCTGGGGGTTATTGGGGGACCAGATTCGCTGTTTAAAGGGCTGGATGCGATGTCGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCCCCCTCGTCGTTGCTTCCGGCATTGGGCTTCGCGGTGTGGGAGTCTTCCCGGCGgtg                                                                                                                                                                                                    >  2:23375/1‑251 (MQ=255)
     cACGCGATGATGCATCAAATGAAAAATCCGGTGCGCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGt                                                                                                                                                                                               <  1:46238/251‑1 (MQ=255)
             gatgCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATGGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGACGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTg                                                                                                                                                                                        <  1:381698/250‑1 (MQ=255)
             gatgCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTAGTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTg                                                                                                                                                                                        <  1:398488/250‑1 (MQ=255)
              atgCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCGGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGc                                                                                                                                                                                      <  2:81209/251‑1 (MQ=255)
                    cagaTGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAg                                                                                                                                                                                 <  1:13931/247‑1 (MQ=255)
                          aaaaaTCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGAt                                                                                                                                                                          <  2:335239/251‑1 (MQ=255)
                                        cgggCTGGGTAGTGGCGGATTGTTTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCGGCGAGTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGACAGTGGCTTTCGCGGTTGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCa                                                                                                                                                              <  2:299899/246‑1 (MQ=255)
                                                                                 ttATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTCGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTAtt                                                                                                                  >  1:348769/1‑251 (MQ=255)
                                                                                             gggtatCGCTGTGAAAAGGGCTGGAGGATATGAAGCAGATGGTGGCAACTATGTCTTCCAGCACCATCGTTTTCTTTTTGCGGCGGGTAATTAACCTTGCCGCGCGGGTCGTTGCGGCCTCCTGTGGCTTGCGCGTTGGCCAAATCTTCCTGGAATTGTTGGGGGGCGTGGCTTTGGGGTTGATTCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTg                                                                                                      <  1:40846/246‑1 (MQ=255)
                                                                                                         ttAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTACACTTGCCGCCCTGGGCGTTGCTGACGGCAGTGGCGTTTGCGGGGGGCGAATCTTTCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCGGCATGGGCACGTTCCCGCCTTACCAGCGGCAATAACCGTTTCTTTCGCGATTCTCGGCATCGTGCTGGTGGTAac                                                                                           >  1:35391/1‑250 (MQ=255)
                                                                                                           aaaGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGGCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGGGGCATTAGGGGTGATGCTGCATTAGCACGTTCCCGCCGTACCAGGGGCAATAACCGTTTCTTGCGGCTTTCTCGGCCTCGGGCTGGGGGTAACAcg                                                                                        >  2:81171/1‑251 (MQ=255)
                                                                                                                 cTGGATGAGATGCAGAAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGAAATTAAACTTTCCGCCCTGGTCGTTGCTGCCGCCAGTGGGTTTCGCGGTGGGAGTATCTTCCTGGCAGTGTTTTGCGGCGTGGCATTAGGGTTGATGCTGAAAGGGCACGTTCCCGCCGTACCAGCGGCAATAACCGATTCTTGCGCTATTCTCGGCCACGTGCTGGTGGTaaaaccgggtgg                                                                                  >  1:312619/1‑241 (MQ=255)
                                                                                                                                cagATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGCAATTAAACTTGCCCCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGAGGGCGAATCTTCCTGGCACAGTTTGTCGGCGGGGCATTTGGGTTGATGCTGCAAGAACACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCtttt                                                                   >  1:186404/1‑251 (MQ=255)
                                                                                                                                            aaTCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCCGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATggcggcggg                                                       >  2:4132/1‑248 (MQ=255)
                                                                                                                                             aTCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATggcggccgtg                                                      >  1:15497/1‑247 (MQ=255)
                                                                                                                                                aGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATggcggccggcgtt                                                   >  1:154621/1‑244 (MQ=255)
                                                                                                                                                aGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGAACGTTCCCGCCGGACCAGCGGCAATAACCGTTTCTTTCGATATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAggtggt                                                   >  1:40287/1‑246 (MQ=255)
                                                                                                                                                                                   gTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTTTAt                >  1:385033/1‑251 (MQ=255)
                                                                                                                                                                                             cTTGCCGCCCTGGTCGTTGCTGCCGCAAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCCTTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGCTGGCTGGTTAAGTCTTTTTa                                                                 >  2:376254/1‑193 (MQ=255)
                                                                                                                                                                                                  cgtttgggTCGTTGGTGCCGCCGGTGGCTTCCGGGGGGGGCGATTCTTCCTGGCGGTGTTGGCCGGGGGGGCATTGGGGTTGATGCTGCATGAGCCCGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCgg  <  1:5613/245‑1 (MQ=255)
                                                                                                                                                                                                        ggTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCCGGCGTTGTAACCCATAACCCATTGCTACCGCTGCTCTGTATCGTCATGCTTCCgg  >  2:5613/1‑245 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                         
GGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGG  >  NC_000913/2510366‑2510809

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