Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,912,429 C→T 23.9% A135T (GCG→ACG)  pyrI ← aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,912,4290CT23.9% 91.1 / 23.8 46A135T (GCG→ACG) pyrIaspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base C (13/22);  new base T (6/5);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTG  >  minE/2912363‑2912499
                                                                  |                                                                      
aCCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa                                                                      >  1:261017/1‑69 (MQ=255)
  cGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGatat                                                                   >  1:738073/1‑70 (MQ=255)
    cAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAt                                                                <  1:135103/71‑1 (MQ=255)
         aattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgcg                                                           >  1:432723/1‑71 (MQ=255)
           ttGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgc                                                            <  1:754081/68‑1 (MQ=255)
              gCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGatat                                                                   >  1:718728/1‑58 (MQ=255)
              gCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGatat                                                                   >  1:4802/1‑58 (MQ=255)
              gCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGttt                                                      <  1:1744947/71‑1 (MQ=255)
                 aGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCg                                                   >  1:1792864/1‑71 (MQ=255)
                 aGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCCTTTTCg                                                   <  1:284599/71‑1 (MQ=255)
                  gCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGa                                                                      <  1:1260931/51‑1 (MQ=255)
                   caccacATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCa                                                 >  1:606914/1‑71 (MQ=255)
                   caccacATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa                                                                      <  1:894435/50‑1 (MQ=255)
                    accacATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCAc                                                >  1:8927/1‑71 (MQ=255)
                       acaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                              >  1:321195/1‑70 (MQ=255)
                        caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCTCGATATc                                                                  <  1:1607534/49‑1 (MQ=255)
                        caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGc                                             >  1:1862747/1‑70 (MQ=255)
                          ttATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGatat                                                                   <  1:2037611/46‑1 (MQ=255)
                          ttATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATc                                                                  <  1:564482/47‑1 (MQ=255)
                           tATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCa                                                                 >  1:1546264/1‑47 (MQ=255)
                           tATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCaaac                                         >  1:433307/1‑70 (MQ=255)
                           tATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAg                                         >  1:563178/1‑71 (MQ=255)
                            aTGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGc                                        >  1:1460679/1‑71 (MQ=255)
                             tGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGtttt                                                     >  1:1719454/1‑57 (MQ=255)
                               ggAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATc                                                                  <  1:1947752/42‑1 (MQ=255)
                               ggAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATc                                                                  <  1:1695175/42‑1 (MQ=255)
                               ggAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATc                                                                  <  1:448748/42‑1 (MQ=255)
                               ggAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTgg                                     <  1:673316/71‑1 (MQ=255)
                                  aaaCTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGCCAAAGCTGGATg                                  <  1:593037/71‑1 (MQ=255)
                                   aaCTCTTTTTCACAGTATTTGCATTTGAGAGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGa                                 <  1:2046234/71‑1 (MQ=255)
                                    aCTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGaa                                <  1:1964004/71‑1 (MQ=255)
                                    aCTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGa                                 <  1:1271266/70‑1 (MQ=255)
                                     ctctTTTTCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGaaa                               >  1:906837/1‑71 (MQ=255)
                                        tttttCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCg                            >  1:1880122/1‑71 (MQ=255)
                                        tttttCACAGTATTTGCATTTGAGCGCGATATCAtt                                                               <  1:39912/36‑1 (MQ=255)
                                          tttCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAAcc                             <  1:759058/68‑1 (MQ=255)
                                          tttCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCg                            <  1:1164606/69‑1 (MQ=255)
                                            tCACAGTATTTGCATTTGAGCGCGATATCATTTGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGtt                         <  1:1122150/70‑1 (MQ=255)
                                              acaGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCgg                      <  1:88153/71‑1 (MQ=255)
                                                  tATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTg                                      <  1:1616031/51‑1 (MQ=255)
                                                    tttGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAt                   <  1:714266/68‑1 (MQ=255)
                                                      tGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATg                  >  1:1616690/1‑67 (MQ=255)
                                                        cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa             <  1:288496/70‑1 (MQ=255)
                                                         aTTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATa           <  1:322211/71‑1 (MQ=255)
                                                         aTTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGATGATa           <  1:1742333/71‑1 (MQ=255)
                                                           ttGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATa           >  1:367367/1‑69 (MQ=255)
                                                                  cgATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTg  >  1:2141700/1‑71 (MQ=255)
                                                                  |                                                                      
ACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTG  >  minE/2912363‑2912499

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