Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,912,429 C→T 35.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,4290CT35.9% 97.1 / 92.3 92A135T (GCG→ACG) pyrIaspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base C (34/25);  new base T (23/10);  total (57/35)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

CCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAG  >  minE/2912360‑2912493
                                                                     |                                                                 
ccAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                   <  1:2064041/70‑1 (MQ=255)
   aCCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                               >  1:964516/1‑71 (MQ=255)
   aCCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                               >  1:3590656/1‑71 (MQ=255)
          ttaattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACATTATTTGCATTTGAGCGCGatat                                                              <  1:1260040/65‑1 (MQ=255)
          ttaattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCAt                                                           <  1:1616209/68‑1 (MQ=255)
          ttaattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCAt                                                           <  1:2290086/68‑1 (MQ=255)
          ttaattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCAt                                                           <  1:2191697/68‑1 (MQ=255)
          ttaattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATc                                                             <  1:1139186/66‑1 (MQ=255)
          ttaattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg                                                        >  1:329704/1‑71 (MQ=255)
          ttaattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg                                                        >  1:324886/1‑71 (MQ=255)
          ttaattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg                                                        >  1:2904254/1‑71 (MQ=255)
          ttaattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg                                                        >  1:1385635/1‑71 (MQ=255)
           taattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgc                                                       >  1:4117360/1‑71 (MQ=255)
            aattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCATTg                                                         >  1:522943/1‑68 (MQ=255)
            aattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCATTg                                                         >  1:2707443/1‑68 (MQ=255)
            aattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCATTg                                                         >  1:2551459/1‑68 (MQ=255)
            aattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAt                                                           >  1:2678439/1‑66 (MQ=255)
            aattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAt                                                           >  1:1157822/1‑66 (MQ=255)
            aattGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTAGCATTTGAGCGTGATATCATTg                                                         >  1:1061596/1‑68 (MQ=255)
                ggCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCTATATCATTGgcgc                                                     >  1:2017634/1‑68 (MQ=255)
                ggCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtt                                                  >  1:2437108/1‑71 (MQ=255)
                ggCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtt                                                  >  1:2879982/1‑71 (MQ=255)
                ggCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtt                                                  >  1:2445943/1‑71 (MQ=255)
                 gCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGata                                                               <  1:3152953/57‑1 (MQ=255)
                 gCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGata                                                               <  1:3675213/57‑1 (MQ=255)
                 gCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGata                                                               <  1:632041/57‑1 (MQ=255)
                   cAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCAt                                                           <  1:854705/59‑1 (MQ=255)
                     gCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCg                                              >  1:1304614/1‑70 (MQ=255)
                       accacATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATc                                                             <  1:1988123/53‑1 (MQ=255)
                        ccacATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                   <  1:1808550/46‑1 (MQ=255)
                         cacaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                   <  1:3590284/45‑1 (MQ=255)
                         cacaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                   <  1:383145/45‑1 (MQ=255)
                         cacaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                   <  1:2709613/45‑1 (MQ=255)
                         cacaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                   <  1:400074/45‑1 (MQ=255)
                         cacaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                   <  1:3322964/45‑1 (MQ=255)
                         cacaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                   <  1:1368700/45‑1 (MQ=255)
                         cacaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                   <  1:970115/45‑1 (MQ=255)
                          acaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGc                                        <  1:138714/71‑1 (MQ=255)
                          acaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                         >  1:784327/1‑70 (MQ=255)
                          acaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                         >  1:2588901/1‑70 (MQ=255)
                          acaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                         >  1:1129937/1‑70 (MQ=255)
                          acaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                         >  1:1138894/1‑70 (MQ=255)
                          acaTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                         >  1:1931202/1‑70 (MQ=255)
                            attattGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGatat                                                              <  1:1529221/47‑1 (MQ=255)
                            aTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGa                                                                 <  1:645232/44‑1 (MQ=255)
                             ttATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAtt                                                          <  1:1642185/50‑1 (MQ=255)
                             ttATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAtt                                                          <  1:1778403/50‑1 (MQ=255)
                             ttATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAtt                                                          <  1:308301/50‑1 (MQ=255)
                              tATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                               >  1:8942/1‑44 (MQ=255)
                              tATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                               >  1:3990477/1‑44 (MQ=255)
                              tATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                               >  1:1831537/1‑44 (MQ=255)
                                  ggAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCa                                                            >  1:3623444/1‑43 (MQ=255)
                                  ggAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCa                                                            >  1:3257843/1‑43 (MQ=255)
                                   gAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa                                                                 >  1:3705553/1‑37 (MQ=255)
                                   gAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa                                                                 >  1:287165/1‑37 (MQ=255)
                                     aaaCTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtttt                                                >  1:2347275/1‑52 (MQ=255)
                                      aaCTCTTTTTCACATTATTTGCATTTGAGCGTGATATc                                                             <  1:441861/38‑1 (MQ=38)
                                       aCTCTTTTTCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGaa                          <  1:3686885/71‑1 (MQ=255)
                                        ctctTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGaaa                         <  1:3308653/71‑1 (MQ=255)
                                          ctTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGa                              >  1:2990516/1‑64 (MQ=255)
                                          ctTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGa                              >  1:2252724/1‑64 (MQ=255)
                                           tttttCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGc                                             >  1:2409409/1‑49 (MQ=255)
                                           tttttCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCg                      >  1:3317074/1‑71 (MQ=255)
                                           tttttCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCg                      >  1:3508241/1‑71 (MQ=255)
                                             tttCACAGTATTTGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                     <  1:3618520/70‑1 (MQ=255)
                                             tttCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                     <  1:1191959/70‑1 (MQ=255)
                                                  caGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGc               <  1:2010960/71‑1 (MQ=255)
                                                     tATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCgg                <  1:1538575/67‑1 (MQ=255)
                                                     tATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCgg                <  1:1408860/67‑1 (MQ=255)
                                                     tATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAt             <  1:4030871/70‑1 (MQ=255)
                                                     tATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAt             <  1:1588314/70‑1 (MQ=255)
                                                     tATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAt             <  1:3820367/70‑1 (MQ=255)
                                                      aTTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCg                                              <  1:872847/37‑1 (MQ=255)
                                                      aTTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCg                                              <  1:2702091/37‑1 (MQ=255)
                                                       tttGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGc          >  1:1472507/1‑71 (MQ=255)
                                                       tttGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGc          >  1:1140248/1‑71 (MQ=255)
                                                       tttGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGc          >  1:3161124/1‑71 (MQ=255)
                                                       tttGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGc          >  1:1872261/1‑71 (MQ=255)
                                                       tttGCATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGc          >  1:3451941/1‑71 (MQ=255)
                                                       tttGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGc          >  1:1362155/1‑71 (MQ=255)
                                                       tttGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGc          >  1:2331646/1‑71 (MQ=255)
                                                        ttGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCa                                            >  1:3048840/1‑37 (MQ=255)
                                                        ttGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCa                                            >  1:3354939/1‑37 (MQ=255)
                                                         tGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCt         <  1:1212955/70‑1 (MQ=255)
                                                         tGCATTTGAGCGCGATATCATTGGCGCCTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCt         <  1:3864400/70‑1 (MQ=255)
                                                             tttGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGc                                  <  1:2263400/41‑1 (MQ=255)
                                                             tttGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTg                                <  1:3900055/43‑1 (MQ=255)
                                                             tttGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATAc    <  1:2954653/71‑1 (MQ=255)
                                                             tttGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATAc    <  1:1523787/71‑1 (MQ=255)
                                                              ttGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATa     >  1:492105/1‑69 (MQ=255)
                                                              ttGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATa     >  1:3705256/1‑69 (MQ=255)
                                                              ttGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATa     >  1:71025/1‑69 (MQ=255)
                                                              ttGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATa     >  1:3491964/1‑69 (MQ=255)
                                                              ttGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAAGCTGGATGAAACCGGTTCGGCATg            >  1:1706521/1‑63 (MQ=255)
                                                              ttGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATa     >  1:2189622/1‑69 (MQ=255)
                                                               tGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                     >  1:801604/1‑52 (MQ=255)
                                                               tGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                     >  1:2096147/1‑52 (MQ=255)
                                                               tGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                     >  1:3580398/1‑52 (MQ=255)
                                                               tGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                     >  1:3886516/1‑52 (MQ=255)
                                                               tGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAg  >  1:1338453/1‑71 (MQ=255)
                                                                     |                                                                 
CCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAG  >  minE/2912360‑2912493

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: