Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 310 (T)9→10 intergenic (+55/‑27) thrL → / → thrA thr operon leader peptide/fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133011.T93.8% 109.4 / ‑0.6 32intergenic (+46/‑36)thrL/thrAthr operon leader peptide/fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (22/8);  total (23/9)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGC  >  NC_000913/225‑371
                                                                             |                                                                       
caccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtttt                                                                     >  1:46/1‑81 (MQ=255)
 accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtttt                                                                     <  1:13/80‑1 (MQ=255)
            aGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGAcc                                                          >  2:47/1‑80 (MQ=255)
            aGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGAcc                                                          >  1:50/1‑80 (MQ=255)
             ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCaa                                                        <  2:14/81‑1 (MQ=255)
             ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCa                                                         >  1:62/1‑80 (MQ=255)
                aaCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAgg                                                     >  1:60/1‑81 (MQ=255)
                 aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGt                                                    <  2:15/81‑1 (MQ=255)
                   ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTaa                                                  >  2:67/1‑81 (MQ=255)
                   ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTaa                                                  <  2:16/81‑1 (MQ=255)
                   ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTaa                                                  >  1:66/1‑81 (MQ=255)
                                gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                      >  2:51/1‑80 (MQ=255)
                                 cgTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                    <  1:17/81‑1 (MQ=255)
                                 cgTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                    >  1:59/1‑81 (MQ=255)
                                     cAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                 >  1:92/1‑80 (MQ=255)
                                       ggAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                              >  1:73/1‑81 (MQ=255)
                                        gAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                             >  1:68/1‑81 (MQ=255)
                                         aaaCACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                            >  1:63/1‑81 (MQ=255)
                                              caGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                         >  2:56/1‑79 (MQ=255)
                                                 aaaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGtt                    <  1:19/81‑1 (MQ=255)
                                                  aaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                   >  1:61/1‑81 (MQ=255)
                                                    aaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                 >  1:57/1‑81 (MQ=255)
                                                      aGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg               >  2:64/1‑81 (MQ=255)
                                                      aGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg               >  2:58/1‑81 (MQ=255)
                                                       gCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcgg              >  2:83/1‑81 (MQ=255)
                                                        cccGCACCTGACAGTGCGGGCTTTTGTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt             <  2:20/81‑1 (MQ=255)
                                                         ccGCACCTGACAGTGCGGGCGTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa            >  1:55/1‑80 (MQ=255)
                                                          cGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa            >  1:76/1‑80 (MQ=255)
                                                          cGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc           <  1:21/81‑1 (MQ=255)
                                                             aCCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt         >  1:70/1‑80 (MQ=255)
                                                               cTGACAGTGCGGGCTTTTTTTTTTCGACCAAGGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg      <  2:22/81‑1 (MQ=255)
                                                               cTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg      >  1:80/1‑81 (MQ=255)
                                                                  aCAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg    <  1:23/80‑1 (MQ=255)
                                                                    aGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc  >  1:82/1‑80 (MQ=255)
                                                                             |                                                                       
CACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGC  >  NC_000913/225‑371

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