Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,410,419 G→A L312L (CTG→TTG)  yfbS ← putative transporter YfbS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,410,4190GA95.0% 43.5 / NA 21L312L (CTG→TTG) yfbSputative transporter YfbS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/9);  minor base C (1/0);  total (11/10)
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

AACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCA  >  NC_000913/2410207‑2410655
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
aaCCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CACTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGcc                                                                                                                                                                                                        >  1:184160/1‑251 (MQ=255)
aaCCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGACCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGcc                                                                                                                                                                                                        >  1:146917/1‑251 (MQ=255)
    aGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGgag                                                                                                                                                                                                    <  1:139409/251‑1 (MQ=255)
     ggATGATATCGCCCCGCCGCCGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGgaga                                                                                                                                                                                                   <  2:104381/251‑1 (MQ=255)
         gatAACGCCCGGCAGCATAGCCTCCACCGCCAGCGACCCTTCCAGCGCCTCGCAATTGCGCTTGAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGAGCCGAGTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCGCATCAAGGGCCTTGTCAGAAAAATACTCGCCGCG‑CAATACCAGCGGCTCCAGCAAGTGCTCGCTTCAAAATTGCCGGAGATCGa                                                                                                                                                                                               <  2:211801/251‑1 (MQ=255)
                       gcagAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGAc                                                                                                                                                                                 <  1:67459/251‑1 (MQ=255)
                                cagtcgcacgCGAACCTTCCAGCCTCACGCCATTGCGCTTGAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATCGACCAATCAGTTCTGACGCAGGAACTAATGAAATCTCTGCCATGCCCACATTAAGGGCCTGGTCAGAAAAATACTCGCCGCT‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCAGTAGATCGACATTATCCGCTGACATATCAATa                                                                                                                                                                        <  2:57817/242‑1 (MQ=255)
                                    cGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGGCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACAAATCAATAAGc                                                                                                                                                                     >  2:185753/1‑250 (MQ=255)
                                       cAGCGAACCTTCAAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTAAGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACAACAGACGCAGACATATAAATAAGaaaa                                                                                                                                                                  >  2:1382/1‑250 (MQ=255)
                                              ccTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCAcg                                                                                                                                                          <  1:185753/251‑1 (MQ=255)
                                                            gccaTTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTAGCGCACCGATTTACCAATCAGTTCCGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGAGAAATTc                                                                                                                                              >  1:238345/1‑249 (MQ=255)
                                                                                            cGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGcgacg                                                                                                             >  1:295254/1‑250 (MQ=255)
                                                                                                  gggTACGAAATCAAATTTCGCGCCCCTATTTACCAACCAGTTCTGACTCAGGAATTAATGAAACCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCTAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAAcgc                                                                                                      <  2:261812/251‑1 (MQ=255)
                                                                                                                           cGATTTACCAATCAGTTCTGAATCAGGAATTAATGAAATCTATGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGt                                                                             >  1:277629/1‑251 (MQ=255)
                                                                                                                              tttACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGtt                                                                            >  1:221061/1‑249 (MQ=255)
                                                                                                                                 aCCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATACCGCGACGAAAACGCCGCCAGCGTTCAACACCTATGACGTTAgcgc                                                                       >  2:102872/1‑251 (MQ=255)
                                                                                                                                                       agtaataaAATCTCCGTCCTGCCCCCATCAAGGGTCTGGTCATAAAAATGCTCGCCTCGACA‑TACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGAGCGACATCAGCCGCACACATATCAATAAGCAAAACGTCACTCGCGCGAAATTCAGAAAACCCATTAACGTCCACGTTAACGCGACGAAAACGCCGCCAGCGTTCAACGCCGATGACGTTAGCGCCATAACGCTCACGTAATTTgag                                                 <  2:82044/249‑1 (MQ=255)
                                                                                                                                                                                    aaGGGCCTGGTCAGAAAAAGACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCg                    <  2:132083/251‑1 (MQ=255)
                                                                                                                                                                                           tGGTCAGAAAAATACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTAACTACAAAATTGCCGGAGATCGACATCAGCCGCAGCCATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCAAAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCgcatcggg             >  2:233773/1‑244 (MQ=255)
                                                                                                                                                                                                     aaTACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGcc   >  1:304732/1‑251 (MQ=255)
                                                                                                                                                                                                      aTACTCGCCGCG‑CAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCa  <  1:94870/251‑1 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
AACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCG‑CAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCA  >  NC_000913/2410207‑2410655

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