Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,447 T→C E35G (GAG→GGG)  araD ← L‑ribulose‑5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,4470TC97.8% 154.6 / ‑5.8 46E35G (GAG→GGG) araDL‑ribulose‑5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (24/21);  total (25/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

GGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTT  >  NC_000913/66333‑66582
                                                                                                                  |                                                                                                                                       
ggCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGAGCCCGCGATCAACGGCGCTGACGTTGc                                                                                                                 >  1:98561/1‑139 (MQ=255)
 gCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGcc                                                                                                                <  1:1068553/139‑1 (MQ=255)
 gCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGcc                                                                                                                <  1:126891/139‑1 (MQ=255)
 gCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGcc                                                                                                                <  2:97958/139‑1 (MQ=255)
     tttCGTACCTTCAACCACTTCACCGGTTCCGATGATAACCACGACCATATCGTCAGCGGTCCTGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCAt                                                                                                            >  2:694577/1‑139 (MQ=255)
          tACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAg                                                                                                       <  2:966740/139‑1 (MQ=255)
                 aaCCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGAcc                                                                                                >  1:516202/1‑139 (MQ=255)
                      cTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACg                                                                                                                     <  2:456968/113‑1 (MQ=255)
                      cTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACg                                                                                                                     >  1:456968/1‑113 (MQ=255)
                       ttCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGAc                                                                                                                      <  1:607406/111‑1 (MQ=255)
                       ttCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGAc                                                                                                                      >  2:607406/1‑111 (MQ=255)
                         cACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAAc                                                                                                                               <  2:506692/100‑1 (MQ=255)
                         cACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAAc                                                                                                                               >  1:506692/1‑100 (MQ=255)
                                        tAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAg                                                                         >  1:762270/1‑139 (MQ=255)
                                         aaCCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAgg                                                                        >  1:973453/1‑139 (MQ=255)
                                             aCGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTgg                                                                    <  2:762270/139‑1 (MQ=255)
                                               gACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGTAGCGCCAGGTTGGCt                                                                  >  2:181969/1‑139 (MQ=255)
                                               gACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTg                                                                     >  1:409816/1‑136 (MQ=255)
                                               gACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTg                                                                     <  2:409816/136‑1 (MQ=255)
                                               gACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCt                                                                  >  2:781934/1‑139 (MQ=255)
                                                aCCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTg                                                                     <  2:237875/135‑1 (MQ=255)
                                                aCCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTg                                                                     >  1:237875/1‑135 (MQ=255)
                                                           aGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGgcgc                                                                                                                          >  1:714603/1‑71 (MQ=255)
                                                           aGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGgcgc                                                                                                                          <  2:714603/71‑1 (MQ=255)
                                                           aGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTgg                                                      >  1:688044/1‑139 (MQ=255)
                                                             cGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGt                                                                                                                    >  1:783536/1‑75 (MQ=255)
                                                             cGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGt                                                                                                                    <  2:783536/75‑1 (MQ=255)
                                                                 cATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACAAGGTTGTGTTTTGGCAGCGCCAGGTTGGCT‑‑TAATACCTGGCGTTTgag                                              <  2:98561/139‑1 (MQ=255)
                                                                  aTGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTga                                               <  2:973453/139‑1 (MQ=255)
                                                                   tGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTgag                                              <  2:688044/139‑1 (MQ=255)
                                                                          gTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGAtcttct                                       <  2:868721/139‑1 (MQ=255)
                                                                             aTCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACga                                                                                                                    <  2:671373/59‑2 (MQ=255)
                                                                             aTCGACGCCGGAAGGTTGGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGt                                                                                                                    >  1:671373/1‑59 (MQ=39)
                                                                                     cGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGgcgc                                                                                                                          <  1:331717/45‑1 (MQ=255)
                                                                                     cGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGgcgc                                                                                                                          >  2:331717/1‑45 (MQ=255)
                                                                                        aaGGTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTcc                         >  2:445819/1‑139 (MQ=255)
                                                                                          ggTTTGATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCtt                       <  2:467231/139‑1 (MQ=255)
                                                                                               gATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAg                                                                                                       >  1:644040/1‑54 (MQ=255)
                                                                                               gATCACAAAGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAg                                                                                                       <  2:644040/54‑1 (MQ=255)
                                                                                                     aaaGACGCCGCGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATg            >  2:884466/1‑139 (MQ=255)
                                                                                                          cgccgcGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAg                                                                                                       >  1:944811/1‑43 (MQ=255)
                                                                                                          cgccgcGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAg                                                                                                       <  2:944811/43‑1 (MQ=255)
                                                                                                            ccgcGCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTg     >  2:417869/1‑139 (MQ=255)
                                                                                                             cgcgCCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTgt    >  1:27839/1‑138 (MQ=255)
                                                                                                               cgcCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTaa                                     >  2:27123/1‑104 (MQ=255)
                                                                                                               cgcCCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCtt  >  1:100152/1‑139 (MQ=255)
                                                                                                                  |                                                                                                                                       
GGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTT  >  NC_000913/66333‑66582

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