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

AACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCA  >  NC_000913/2410207‑2410655
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
aaCCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCACTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGcc                                                                                                                                                                                                        >  1:184160/1‑251 (MQ=255)
aaCCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGACCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGcc                                                                                                                                                                                                        >  1:146917/1‑251 (MQ=255)
    aGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGgag                                                                                                                                                                                                    <  1:139409/251‑1 (MQ=255)
     ggATGATATCGCCCCGCCGCCGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGgaga                                                                                                                                                                                                   <  2:104381/251‑1 (MQ=255)
         gatAACGCCCGGCAGCATAGCCTCCACCGCCAGCGACCCTTCCAGCGCCTCGCAATTGCGCTTGAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGAGCCGAGTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCGCATCAAGGGCCTTGTCAGAAAAATACTCGCCGCGCAATACCAGCGGCTCCAGCAAGTGCTCGCTTCAAAATTGCCGGAGATCGa                                                                                                                                                                                               <  2:211801/251‑1 (MQ=255)
                       gcagAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGAc                                                                                                                                                                                 <  1:67459/251‑1 (MQ=255)
                                cagtcgcacgCGAACCTTCCAGCCTCACGCCATTGCGCTTGAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATCGACCAATCAGTTCTGACGCAGGAACTAATGAAATCTCTGCCATGCCCACATTAAGGGCCTGGTCAGAAAAATACTCGCCGCTCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCAGTAGATCGACATTATCCGCTGACATATCAATa                                                                                                                                                                        <  2:57817/242‑1 (MQ=255)
                                    cGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGGCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACAAATCAATAAGc                                                                                                                                                                     >  2:185753/1‑250 (MQ=255)
                                       cAGCGAACCTTCAAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTAAGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACAACAGACGCAGACATATAAATAAGaaaa                                                                                                                                                                  >  2:1382/1‑250 (MQ=255)
                                              ccTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCAcg                                                                                                                                                          <  1:185753/251‑1 (MQ=255)
                                                            gccaTTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTAGCGCACCGATTTACCAATCAGTTCCGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGAGAAATTc                                                                                                                                              >  1:238345/1‑249 (MQ=255)
                                                                                            cGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGcgacg                                                                                                             >  1:295254/1‑250 (MQ=255)
                                                                                                  gggTACGAAATCAAATTTCGCGCCCCTATTTACCAACCAGTTCTGACTCAGGAATTAATGAAACCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCTAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAAcgc                                                                                                      <  2:261812/251‑1 (MQ=255)
                                                                                                                           cGATTTACCAATCAGTTCTGAATCAGGAATTAATGAAATCTATGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGt                                                                             >  1:277629/1‑251 (MQ=255)
                                                                                                                              tttACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGtt                                                                            >  1:221061/1‑249 (MQ=255)
                                                                                                                                 aCCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATACCGCGACGAAAACGCCGCCAGCGTTCAACACCTATGACGTTAgcgc                                                                       >  2:102872/1‑251 (MQ=255)
                                                                                                                                                       agtaataaAATCTCCGTCCTGCCCCCATCAAGGGTCTGGTCATAAAAATGCTCGCCTCGACATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGAGCGACATCAGCCGCACACATATCAATAAGCAAAACGTCACTCGCGCGAAATTCAGAAAACCCATTAACGTCCACGTTAACGCGACGAAAACGCCGCCAGCGTTCAACGCCGATGACGTTAGCGCCATAACGCTCACGTAATTTgag                                                 <  2:82044/249‑1 (MQ=255)
                                                                                                                                                                                    aaGGGCCTGGTCAGAAAAAGACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCg                    <  2:132083/251‑1 (MQ=255)
                                                                                                                                                                                           tGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTAACTACAAAATTGCCGGAGATCGACATCAGCCGCAGCCATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCAAAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCgcatcggg             >  2:233773/1‑244 (MQ=255)
                                                                                                                                                                                                     aaTACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGcc   >  1:304732/1‑251 (MQ=255)
                                                                                                                                                                                                      aTACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCa  <  1:94870/251‑1 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
AACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCA  >  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: