Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,343 T→C T70A (ACG→GCG)  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,3430TC91.7% 31.5 / NA 12T70A (ACG→GCG) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (3/8);  minor base A (1/0);  total (4/8)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCT  >  NC_000913/66127‑66591
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        
cATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATAGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGACTGCGCCCAGATAGTGGCGTGGCGCGTGTGGGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGTAGTCGGTGGGTTGGCGGATCGGGGGAGGGGTTTTTCGACCCTTCAACCACTTCACCGGTTTCGATGCCAAcc                                                                                                                                                                                                                        >  1:254727/1‑251 (MQ=255)
              aTTTCTGCGTCGGTCATTTTGCGGGGGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATATAGCAGTCGGTGAGTTGGCTTATCGGAGGAGGGCTTGTTCGCACCTTCAACCACtt                                                                                                                                                                                                                                            <  1:83302/217‑1 (MQ=255)
                    gCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGGGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGTGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAgn                                                                                                                                                                                                       <  1:172975/248‑2 (MQ=255)
                                                                               gCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGAGAGTGTGCACAATGCCGCCAATGGACGGGGATGCCTGATAGAGCAGGCGGGGGGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTccg                                                                                                                                                                                                        >  1:114732/1‑186 (MQ=255)
                                                                                   nttGCTGGATTCGACTGGCACGCCTGCGCCCAGATAGTGGCGTGGTGCGCGTGCGTGTGAACAATGCCGACAAAGGACGGGAATGCCTGATAGAGCAGGCGGTGAGTTGGCGTATCGGAGGAGGGGTTTTTCGCACCTTCAACCAATTCACCGGTTTCGATGCTAACAACGACCATATCGTCGGAGGTCATGATGATGCAATAGACGACGGCAGGTTTGATCACAAAGACGCCGCGCTCGAGATCCACggg                                                                                                                                     >  1:51665/2‑250 (MQ=255)
                                                                                                         ccTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCGTCAACCACTTCACCGGTTTCGAGGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGc                                                                                                               <  1:224965/251‑1 (MQ=255)
                                                                                                                   gagagtggcgtggcGCGATTGGGTGGGCACAATGCCGCCAATGGACTCGAATGCCGGATAGAGCAGGCGGTGAGTTGTCGTATCGGAGGAAGGCTTTTTCGCGCCTTCAACCACTTCTCCGGTTTCGATGCTAACCACGACCATATCGTCAGCTGTCATGATGCTTTAATTGACGCCTGAAGTTTGGATGACAAACACCCCTCGCTCGCCATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGt                                                                                                     <  1:271431/247‑1 (MQ=255)
                                                                                                                            gtggcgCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGGGAGTTGTCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGTTCATTATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGgc                                                                                                                                     <  1:144455/210‑1 (MQ=255)
                                                                                                                                                               acgggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTg                                                         <  1:183184/249‑1 (MQ=255)
                                                                                                                                                                                            gTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCGCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAAGACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTg                            <  1:205704/251‑1 (MQ=255)
                                                                                                                                                                                                         anaangGCTTNTTCGNACCTTCAACCACTTCNCCGGTTTCGATGCTAACCANNANCATATNGTCAg                                                                                                                                                                                                        <  2:114732/63‑1 (MQ=255)
                                                                                                                                                                                                          ggaggGCTTTTTCGCACCTTCAACCANTTCACCNNTTNCGATGCTAACCACGACCATNTCNNNANNNGNCATGANGCNNTACTCGACGNCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGNCCCATGTGANNGTGCCCAGGTTATGTTTTGGCNGCGNNANNTNGGcncc                                                                         >  2:61915/1‑189 (MQ=255)
                                                                                                                                                                                                          ggaggGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGGTGGCCTCTAATAc                                                                   <  1:61915/198‑1 (MQ=255)
                                                                                                                                                                                                                       gNACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGGTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCt  >  1:215660/1‑250 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        
CATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCT  >  NC_000913/66127‑66591

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