Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,072,800 C→T 11.1% D298N (GAT→AAT)  epd ← D‑erythrose‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,072,8000CT11.1% 81.6 / 4.3 36D298N (GAT→AAT) epdD‑erythrose‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (16/16);  new base T (2/2);  total (18/18)
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.69e-01

GCTGCTTGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTTTCTTCAC  >  NC_000913/3072660‑3072941
                                                                                                                                            |                                                                                                                                             
gttgctTGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATGCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGa                                                                                                                                             <  2:160273/141‑1 (MQ=255)
                cTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGCGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTCACCCGGGTTTCGGTGCCATCGACAATGGCACTGTGCgg                                                                                                                             >  2:82688/1‑143 (MQ=255)
                cTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAACGCCCCATTCGTTATCGCACCAGACCAACGTGTTGATCAGGTGTACGCCACTGACCCGGGTTTGGGTGCCATCGACACTGGCGCTgcgctc                                                                                                                             >  1:244192/1‑141 (MQ=255)
                      gCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGGTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTg                                                                                                                       >  1:317325/1‑143 (MQ=255)
                            gTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGATGCAATCGACAATGGCACTATGCGCAGCGtcgaaaag                                                                                                                 >  2:472898/1‑136 (MQ=255)
                             tAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTaaaa                                                                                                                <  1:427160/143‑1 (MQ=255)
                              ggCCATAGCTAACGTCGCGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAAt                                                                                                               <  2:382536/142‑1 (MQ=255)
                               gCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATc                                                                                                              <  1:414710/143‑1 (MQ=255)
                                           gtcgtGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCaa                                                                                                  >  2:404310/1‑143 (MQ=255)
                                           gtcgtGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCaa                                                                                                  <  1:390228/143‑1 (MQ=255)
                                                gtCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCa                                                                                             <  2:289390/143‑1 (MQ=255)
                                                    aGCATTCGGCTAGCAAAGCCCCCTCCGTTATCGCACCAGCTCAACGTTTTCATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTc                                                                                         <  1:227139/143‑1 (MQ=255)
                                                      cATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCg                                                                                       <  1:357950/143‑1 (MQ=255)
                                                       aTTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGATGTGCGCCACTGACCCGGGTATGGGTGCCATCGACAATGGCACTGTGCGGCTCGTAGGTACAATCTACAGCGACCAACGGCACTTCCGt                                                                                      >  2:223617/1‑143 (MQ=255)
                                                          cGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATTGTCGATGGCAc                                                                                                                                    >  1:182265/1‑94 (MQ=255)
                                                          cGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATTGTCGATGGCAc                                                                                                                                    <  2:182265/94‑1 (MQ=255)
                                                          cGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATTGTCGATGGCAc                                                                                                                                    <  2:383443/94‑1 (MQ=255)
                                                          cGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATTGTCGATGGCAc                                                                                                                                    >  1:383443/1‑94 (MQ=255)
                                                            tttgtcAAAGCCCCACTCGTTATCCCACCTTATCAACTTTTTGATCATGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCTTGTTTAAAATCTACATAGACCAACGGCAATTCCGTATAGt                                                                                 <  2:448072/138‑1 (MQ=255)
                                                               aGCAATGCCGCATTCGTTCTCGTTCCAGACTAGTGTTTTTATCAGGTGTTCGCCACTGACCCGGGTTTGGGTGCCATCGAGAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAATGGCAATTCCGTATAGTCaa                                                                              <  2:228805/143‑1 (MQ=255)
                                                                                  tCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGGGGGTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCAcc                                                           >  1:37951/1‑143 (MQ=255)
                                                                                      aCCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCCCCTtgt                                                       >  2:328836/1‑143 (MQ=255)
                                                                                       ccAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACACtgtg                                                      >  1:211727/1‑143 (MQ=255)
                                                                                                  ttgtGAGCAGGTGTGCGCCACGGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTgc                                           <  1:325770/140‑1 (MQ=255)
                                                                                                       aTCAGGGGTCCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCaa                                      <  1:458916/143‑1 (MQ=255)
                                                                                                             tgtgCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGtt                                <  1:191179/143‑1 (MQ=255)
                                                                                                                 cgcCACTGGCCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACt                            <  1:530252/143‑1 (MQ=255)
                                                                                                                         aCCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGc                    >  1:260348/1‑143 (MQ=255)
                                                                                                                            cGGGTTTGGGTGGCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCttt                 <  2:194340/143‑1 (MQ=255)
                                                                                                                                      tGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCGTTCGCTTTTACAGGTttc       >  1:498515/1‑143 (MQ=255)
                                                                                                                                      tGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTttc       >  1:41803/1‑143 (MQ=255)
                                                                                                                                      tGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTCCAGCAACCGGTTCACTTCATTGGCATTTACAGCTttc       >  2:411300/1‑143 (MQ=255)
                                                                                                                                         cATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACGCGTtgcttc    >  1:475674/1‑143 (MQ=255)
                                                                                                                                         cATCGACAATGGCAATGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTttcttc    >  1:471998/1‑143 (MQ=255)
                                                                                                                                          aTCGACAATGGCGCTGTGCGGATCGTGGTTATCGGGTGCCGAGGTCAGCGGTAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTTTCTTCa   <  2:402950/143‑1 (MQ=255)
                                                                                                                                           tCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTCGCCTTCATTGGCTTTTACAGGTTTCTTCAc  >  2:173460/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                             
GCTGCTTGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTTTCTTCAC  >  NC_000913/3072660‑3072941

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