Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,194,694 Δ1 bp intergenic (+40/‑229) metA → / → aceB homoserine O‑succinyltransferase/malate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,194,6940A.96.5% 131.9 / ‑4.0 29intergenic (+40/‑229)metA/aceBhomoserine O‑succinyltransferase/malate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (5/23);  minor base G (0/1);  total (5/24)
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

TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCAT  >  REL606/4194624‑4194763
                                                                      |                                                                       
tCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAg                                                                   <  2:1120034/76‑1 (MQ=255)
    cGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCAc                                                              <  2:836940/76‑1 (MQ=255)
       cacaTGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCtt                                                           <  2:596012/76‑1 (MQ=255)
         caTGAATCCAACGCTGGATTAATCTTCTGTGTTAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCg                                                         <  2:598788/76‑1 (MQ=255)
           tGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCggg                                                       <  1:988062/76‑1 (MQ=255)
                   atgCTGGATTAAGCTGTTGTGGTAGGCGATCGTTAAGCCATTCAGCCCCTTCCTCAGGCACCTTCGGGTGCCtttt                                              <  2:943792/74‑1 (MQ=255)
                   aCGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCtttt                                              <  2:943824/76‑1 (MQ=255)
                     gCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCtttttt                                            <  1:1050306/76‑1 (MQ=255)
                     gCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCtttttt                                            <  2:87364/76‑1 (MQ=255)
                     gCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCtttttt                                            <  1:1050291/76‑1 (MQ=255)
                       tGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTAt                                          <  2:660577/76‑1 (MQ=255)
                         gATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTAttt                                        <  2:88445/76‑1 (MQ=255)
                         gATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTAttt                                        <  2:31826/76‑1 (MQ=255)
                         gATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTAttt                                        <  2:568054/76‑1 (MQ=255)
                          aTTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTATTTc                                       <  1:981241/76‑1 (MQ=255)
                             aaTCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCCCGCACCCTCCGGGGGCGCtttttttttttt                                      >  1:243474/1‑74 (MQ=255)
                              aTCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGaa                                   <  2:382222/76‑1 (MQ=255)
                                   ctGTGATAGTCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTa                              <  2:666499/76‑1 (MQ=255)
                                            tCGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTcagcag                     <  2:70044/76‑1 (MQ=255)
                                             ggATCGTTAGGCGATTCAGCGCCTTCCTCAGGCAACTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAgg                    <  1:341654/75‑1 (MQ=255)
                                             cGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAgg                    <  1:341638/76‑1 (MQ=255)
                                             cGATCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGGGCCTTATTTTTTTTCGAAACGTACCTCAGCAgg                    >  1:49510/1‑76 (MQ=255)
                                               aTCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTg                  <  2:296436/76‑1 (MQ=255)
                                               aTCGTTAAGCGATTCAGCACCTGTCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTg                  <  2:274041/76‑1 (MQ=255)
                                                tCGTTAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGa                 <  1:916181/76‑1 (MQ=255)
                                                 cGTTAAGCGATTCAGCACCTTCCTCAGGCACCTNCGGGGGCCTTTTTTATTTTCGANACGTTCCTGCNCAGGTGaa                >  2:1015324/1‑76 (MQ=255)
                                                    tAAGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATaa             <  1:184034/76‑1 (MQ=255)
                                                      aGCGATTCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTTTTTTCCAAAACTACCCCAGCAGGGGAATAAAt           >  2:825467/1‑76 (MQ=255)
                                                           ttCAGCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTAt      <  2:606265/76‑1 (MQ=255)
                                                               gCACCTTCCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCat  >  1:207509/1‑76 (MQ=255)
                                                                      |                                                                       
TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCAT  >  REL606/4194624‑4194763

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