Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,346,688 G→A V213V (GTC→GTT fumB ← fumarase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,6880GA84.4% 75.2 / 2.0 34V213V (GTC→GTTfumBfumarase B
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (9/18);  total (13/21)
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

GAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGAC  >  NC_000913/4346445‑4346924
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
gAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGTACCGAGGGGACGCATTTTCTCAACGAggg                                                                                                                                                                                                                                        >  1:292974/1‑250 (MQ=255)
                    gcgcgcCAAGACCGAGTTTCTGGGCCTCTTCAAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGCCCGTGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGCGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAAGCACAATCGCGATATGGTACGTCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGa                                                                                                                                                                                                                    <  2:284447/251‑1 (MQ=255)
                          cAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagc                                                                                                                                                                                                              <  1:405083/251‑1 (MQ=255)
                                   ggttCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGTCATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTg                                                                                                                                                                                                     <  2:343143/249‑1 (MQ=255)
                                           gCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTg                                                                                                                                                                                             <  2:292974/251‑1 (MQ=255)
                                              tctTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGCTGCAGAACCGAGGGTACGCATTTTCTCAACGAGGAAGATTTTTAGGTTGGCGGGGGTCAGCAGGGGTTTGGGTTCCTGGtg                                                                                                                                                                                          >  1:386113/1‑250 (MQ=255)
                                                      cagTTCCTGTTCCAGCTGGACATCGCGGACCGCTTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTATCGTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGATTCAGCAGGGCTTTGGTTTCCTGGTAGAGATagg                                                                                                                                                                                  <  1:239188/251‑3 (MQ=255)
                                                       agTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGt                                                                                                                                                                                 <  1:247233/251‑1 (MQ=255)
                                                                 ccAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCCGGAACCGGGGGACGCGTTTTTTTTACGAGGGGGGTTTTTAGGTTGCCGGGAGTCAGCAGGGCTTTGGGTTCCTgggcgcggttcgttttgtttgcg                                                                                                                                                                       >  2:183366/1‑230 (MQ=255)
                                                                                         gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTcacaacgc                                                                                                                                                 <  1:111797/249‑1 (MQ=255)
                                                                                                          gccGTCGTCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATc                                                                                                                              <  2:77224/250‑1 (MQ=255)
                                                                                                           ccGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCg                                                                                                                             <  1:106525/251‑1 (MQ=255)
                                                                                                             ggcGGCAGTTCAGCGTAATAGGGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTGCTGCAGACGTACCGCCAAGCACAAACGCGATATGGTACGGCGGTCAGGCTGCAGTACCGAGGGTACGGATTTTCTCAACGAGGAAGGGTTTCAGTTTGCCGGGAGTCAGCAGGGTTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCAGCGcc                                                                                                                           <  1:184613/249‑1 (MQ=255)
                                                                                                             gTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTGTCGCAACGCAAAGGAATTTGTACTCATCGcc                                                                                                                           <  1:137426/251‑1 (MQ=255)
                                                                                                                ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGttt                                                                                                                                                                                                                           <  2:141594/152‑1 (MQ=255)
                                                                                                                ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCAGGACCCAGGGGGCGCCTTTTCTCCACGAGGGAGTTTTTTAGttt                                                                                                                                                                                                                           >  1:141594/1‑152 (MQ=255)
                                                                                                                    gTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAcc                                                                                                                    <  1:258033/251‑1 (MQ=255)
                                                                                                                           gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTCAGTTTGCCGGGGGTCAGCAGGGCTTTGGTTTCCTGGGAGAGATACGTTTTGTTGGGAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTTTACTCATCGCCCTCTACCGCGTAc                                                                                                              >  1:355187/1‑250 (MQ=255)
                                                                                                                              aTAGTGAGCGCTGGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGt                                                                                                          <  2:14044/251‑1 (MQ=255)
                                                                                                                                          tGCTAACTTGACGGTTTTCAGGTTGGTTTCTGCAGCCGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCGCGGGCAACGAGGGGCCGCTCTTTCTCAAAGCGGAAGGTTTTTAGTTTGCCGGGGGTAAGCAGGGGTTTGGGTTCTTGGTCGAGGCAAGTTTTTTTTGCCGAGCCGGCGCCTTTTGCAACGCAACGGACATTGTAATCATTGCCCTCTACCGCGTaacgggcggcttgcgccg                                                                                              >  1:270904/1‑234 (MQ=255)
                                                                                                                                                       gTTTTCAGGTGGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCGTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGTCAGGTTAGTGcc                                                                                <  2:415201/251‑1 (MQ=255)
                                                                                                                                                          ttCAGGTTGGTTTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTTCCGGGGGTCAGCAGGGCTTTTGGTTCCTGGGAGAGATAGGGTTTGTTTGCAGAGCCGCCGCCTTTCGCCACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGGTTTGTGCAGGCAGGTTAGTGCCggg                                                                             >  1:219852/1‑250 (MQ=255)
                                                                                                                                                                ttggttTCTGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCAGTACCGAGGGGAAGCATTTTTTTAACGAGGGAGTTTTTCAGTTTGGCGGGGGTCAGCAGGGGTTTGGTTTTCTGGTAGAGATAGGGTTTGTTTGGAGAGCCGCCGCCCTTCGCAACGCAAAGGCATTTTTACTCATCGCCACCTACCGCGTACAGGGCGCTTTGCGCTAGGCGGGttgtggccgtggttgc                                                                        >  2:93393/1‑237 (MQ=255)
                                                                                                                                                                                tACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGATGCCGTAACGAGGGGAAGCATTTTTTCAACGAGGAAGTTTTTTCGTTTGCCGGGAGAGATCAGGGGTTTGGTTTCCTGGGAGAGATAGGTTTTGTTTGGAGAGACGCCGCCTTTCGCAACGCAAAGGGATTTGTATTCCTCGCCATCTACCCCGGACAGGTAGGTTTGTGCAGGCCGGGTAGTGCCCGTGTTTACCTCTTTTTACATggcc                                                       >  1:143247/1‑248 (MQ=255)
                                                                                                                                                                                 aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCa                                                      <  1:183366/251‑1 (MQ=255)
                                                                                                                                                                                      cAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCCGCACCGAGGGGACGCATTTTTTCAACGAGGAAGGTTTTTAGGTTGCCGGGGGGCAGCAGGGCTTTGGTTTTCTGGTAGAGATAGGGTTTTTTGGGAGAGCCGCCGCCCTTTGCAAAGCCAAGGGATTTTTACTCACCGCCCTCTACCGCGGACAGGGCGATTTGCGCCGGCAGGGTAGTGCCCGGGGTGTCCTCTTTGTACATTTcccccgcc                                                 >  1:2391/1‑245 (MQ=255)
                                                                                                                                                                                              aaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGcgccgcgttttg                                         <  1:93393/251‑7 (MQ=255)
                                                                                                                                                                                                cgcgATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGGCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCGTTTtgtg                                       <  2:2391/251‑1 (MQ=255)
                                                                                                                                                                                                                   cAGGCTGCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCGTTTTGTGAATAGCGCAGGTTATCTTc                    >  2:356257/1‑251 (MQ=255)
                                                                                                                                                                                                                       cTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTAAATAGCGCAGGTTATCTTCGat                 >  2:316694/1‑250 (MQ=255)
                                                                                                                                                                                                                         gCAGTACCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCGTTTTGTGAATAGCGCAGGTTATCTTCGatat               <  1:348510/250‑1 (MQ=255)
                                                                                                                                                                                                                                aCCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTa                                                                                     <  1:4857/175‑1 (MQ=255)
                                                                                                                                                                                                                                aCCGAGGGTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTa                                                                                     >  2:4857/1‑175 (MQ=255)
                                                                                                                                                                                                                                       gTACGCATTTTCTCAACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAGGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCGTTTTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGAc  <  1:279745/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
GAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGAC  >  NC_000913/4346445‑4346924

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