Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,144,345 G→A 31.0% D206N (GAT→AAT)  argG → argininosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,144,3450GA31.0% 122.3 / 74.5 87D206N (GAT→AAT) argGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (33/27);  new base A (11/16);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

CAAAATGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGA  >  minE/2144278‑2144412
                                                                   |                                                                   
cAAAATGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATc                                                                  >  1:750543/1‑71 (MQ=255)
cAAAATGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATc                                                                  <  1:2658746/71‑1 (MQ=255)
  aaaTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTg                                                                <  1:1192767/71‑1 (MQ=255)
   aaTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                               <  1:1580226/71‑1 (MQ=255)
   aaTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                               <  1:1105003/71‑1 (MQ=255)
   aaTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                               <  1:2175003/71‑1 (MQ=255)
   aaTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                               <  1:2492069/71‑1 (MQ=255)
   aaTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                               <  1:1243837/71‑1 (MQ=255)
    aTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                               >  1:3011755/1‑70 (MQ=255)
    aTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                               >  1:3580168/1‑70 (MQ=255)
    aTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTg                                                                >  1:4010377/1‑69 (MQ=255)
             gAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCaa                                                     >  1:1668073/1‑71 (MQ=255)
             gAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCaa                                                     >  1:1111247/1‑71 (MQ=255)
             gAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCaa                                                     >  1:468186/1‑71 (MQ=255)
             gAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCaa                                                     >  1:3067744/1‑71 (MQ=255)
                 aaGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTcc                                                 <  1:2530824/71‑1 (MQ=255)
                      tACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAgg                                                                     <  1:2175890/46‑1 (MQ=255)
                      tACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAgg                                                                     <  1:1596347/46‑1 (MQ=255)
                         tCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                               <  1:1067677/49‑1 (MQ=255)
                         tCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                               <  1:582113/49‑1 (MQ=255)
                         tCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                               <  1:990924/49‑1 (MQ=255)
                         tCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCa                                          >  1:2868906/1‑70 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGa                                                                    >  1:564447/1‑42 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGAt                                                                   >  1:460427/1‑43 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGAt                                                                   >  1:3586204/1‑43 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGAt                                                                   >  1:1813361/1‑43 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATc                                                                  >  1:4092000/1‑44 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTg                                                                >  1:1185145/1‑46 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaaa                                       >  1:3014940/1‑71 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaaa                                       <  1:3718338/71‑1 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaaa                                       >  1:2221489/1‑71 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaaa                                       >  1:1309890/1‑71 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaaa                                       >  1:2478378/1‑71 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaa                                        >  1:1652193/1‑70 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaa                                        <  1:1386910/70‑1 (MQ=255)
                           cacaGACTCCAACATGCTTGGTGCAACGCATGAAGAGAAGGAt                                                                   >  1:963293/1‑43 (MQ=255)
                             caGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATc                                     >  1:2741175/1‑71 (MQ=255)
                                aCTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATGCCTCAACTCCAGCGTCAAAATCGTc                                  <  1:337590/71‑1 (MQ=255)
                                aCTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTc                                  <  1:1493088/71‑1 (MQ=255)
                                aCTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTc                                  <  1:3316588/71‑1 (MQ=255)
                                 cTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTc                                  <  1:2598872/70‑1 (MQ=255)
                                       cATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                               <  1:3447966/35‑1 (MQ=255)
                                       cATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                               <  1:993147/35‑1 (MQ=255)
                                       cATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                               <  1:563971/35‑1 (MQ=255)
                                       cATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                               <  1:355079/35‑1 (MQ=255)
                                       cATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                               <  1:3297660/35‑1 (MQ=255)
                                       cATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                               <  1:3050354/35‑1 (MQ=255)
                                       cATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCaa                                >  1:2793215/1‑66 (MQ=255)
                                         tGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAtt                         <  1:564233/71‑1 (MQ=255)
                                           cTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTa                        <  1:3579563/70‑1 (MQ=255)
                                               gTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAAt                                      >  1:3169697/1‑52 (MQ=255)
                                               gTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAAt                                      >  1:1585351/1‑52 (MQ=255)
                                               gTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAAt                                      >  1:380636/1‑52 (MQ=255)
                                                   aaCGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGa               <  1:1446405/71‑1 (MQ=255)
                                                   aaCGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGGCAAAATCGTCAACCCGATTATGGGCGTGa               <  1:3515623/71‑1 (MQ=255)
                                                      gCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTc                                  >  1:247580/1‑49 (MQ=255)
                                                        aTGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAAtt           >  1:236937/1‑70 (MQ=255)
                                                        aTGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTc          >  1:1768207/1‑71 (MQ=255)
                                                        aTGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTc          >  1:712348/1‑71 (MQ=255)
                                                        aTGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTc          >  1:1631562/1‑71 (MQ=255)
                                                        aTGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTc          >  1:1651561/1‑71 (MQ=255)
                                                        aTGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCaa                                >  1:1408142/1‑49 (MQ=255)
                                                        aTGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCaa                                >  1:491018/1‑49 (MQ=255)
                                                        aTGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCaa                                >  1:1531974/1‑49 (MQ=255)
                                                        aTGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCaa                                >  1:441695/1‑49 (MQ=255)
                                                         tGAAGCGAAGGATCTGGAATACGTCAACTCCAGCGTCAAAATCGTCAACCCGATTa                        >  1:3782180/1‑56 (MQ=255)
                                                         tGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTc                                  <  1:291275/46‑1 (MQ=255)
                                                         tGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTa                        >  1:2270747/1‑56 (MQ=255)
                                                         tGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaa                                        <  1:47821/40‑1 (MQ=255)
                                                         tGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaa                                        <  1:2392283/40‑1 (MQ=255)
                                                         tGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaa                                        <  1:576838/40‑1 (MQ=255)
                                                         tGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaa                                        <  1:2017053/40‑1 (MQ=255)
                                                         tGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaa                                        <  1:4095308/40‑1 (MQ=255)
                                                         tGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaa                                        <  1:162653/40‑1 (MQ=255)
                                                         tGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaa                                        <  1:3309164/40‑1 (MQ=255)
                                                           aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGa                           >  1:2336028/1‑51 (MQ=255)
                                                           aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTgg       >  1:334027/1‑71 (MQ=255)
                                                           aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTgg       <  1:1107149/71‑1 (MQ=255)
                                                           aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTgg       <  1:2027907/71‑1 (MQ=255)
                                                           aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTgg       <  1:3309709/71‑1 (MQ=255)
                                                           aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTgg       >  1:1498892/1‑71 (MQ=255)
                                                           aaGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATg                      >  1:1788434/1‑56 (MQ=255)
                                                              cGAAGGATCTTGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGaaa             <  1:3155840/62‑1 (MQ=255)
                                                               gaagaaTCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAAtt           >  1:22370/1‑63 (MQ=255)
                                                               gAATGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATg   <  1:3223358/71‑1 (MQ=255)
                                                               gAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATg   <  1:1026094/71‑1 (MQ=255)
                                                               gAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGAGTGAAATTCTGGGATg   <  1:2694630/71‑1 (MQ=255)
                                                                aaGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATga  <  1:3331910/71‑1 (MQ=255)
                                                                 aGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGAt    >  1:231467/1‑68 (MQ=255)
                                                                 aGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCAGATTATGGGCGTGAAATTCTGGGa     >  1:2468287/1‑67 (MQ=255)
                                                                   |                                                                   
CAAAATGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGA  >  minE/2144278‑2144412

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