Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,079 T→C I158V (ATC→GTC)  araD ← L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,0790TC92.0% 68.9 / ‑3.6 25I158V (ATC→GTC) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (8/15);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

AATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGT  >  NC_000913/65955‑66202
                                                                                                                            |                                                                                                                           
aaTATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGtt                                                                                                                     >  2:375957/1‑133 (MQ=255)
     ccccATATCAGCGACCTCTTCCAGCACGATCGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCa                                                                                                                <  1:131943/133‑1 (MQ=255)
     ccccATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCa                                                                                                                <  1:364952/133‑1 (MQ=255)
      cccATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCaa                                                                                                               <  1:398583/133‑1 (MQ=255)
              aGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCt                                                                                                       >  2:47235/1‑133 (MQ=255)
              aGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCt                                                                                                       >  1:417628/1‑133 (MQ=255)
                  aCCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCACCGACACCTTGTCTTTCAAAGGCTTCTACGa                                                                                                   >  1:216775/1‑133 (MQ=255)
                   ccTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGAt                                                                                                  <  2:481630/133‑1 (MQ=255)
                   ccTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGAt                                                                                                  >  1:473264/1‑133 (MQ=255)
                   ccTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGAt                                                                                                  >  1:271692/1‑133 (MQ=255)
                       ttCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACg                                                                                              <  2:390527/133‑1 (MQ=255)
                             cACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCg                                                                                        >  1:113357/1‑133 (MQ=255)
                             cACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCACATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCg                                                                                        <  2:384610/133‑1 (MQ=255)
                                gATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGtt                                                                                     <  2:113357/133‑1 (MQ=255)
                                  tGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTc                                                                                   >  1:57422/1‑133 (MQ=255)
                                            cACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCatat                                                                         <  2:422917/133‑1 (MQ=255)
                                                             ttttGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGc                                                        <  2:417628/133‑1 (MQ=255)
                                                               ttGCCCCATGAAAACGGGCCGTGGTAATCGACCAGAACGCCCGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGt                                                      <  2:216775/133‑1 (MQ=255)
                                                                            aCGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCg                                         <  1:413130/133‑1 (MQ=255)
                                                                             cGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCgg                                        <  1:339373/133‑1 (MQ=255)
                                                                                gCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTg                                     <  1:47235/133‑1 (MQ=255)
                                                                                 ccGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGc                                    >  2:355207/1‑133 (MQ=255)
                                                                                       gAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGaa                              <  2:271692/133‑1 (MQ=255)
                                                                                           ggACCAGAACGCCGGGCATTTGCGCTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGt                          <  1:375957/133‑1 (MQ=255)
                                                                                                                   cTGCATCGACACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggt  <  2:320377/133‑1 (MQ=255)
                                                                                                                            |                                                                                                                           
AATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGT  >  NC_000913/65955‑66202

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