Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,348 T→C 5.7% T116T (ACA→ACG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,3480TC 5.7% 246.3 / 4.8 87T116T (ACA→ACGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (30/52);  new base C (3/2);  total (33/54)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAG  >  minE/2919280‑2919418
                                                                    |                                                                      
agcagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAgtgt                                                                      <  1:2925766/71‑1 (MQ=255)
agcagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAgtgt                                                                      <  1:1736284/71‑1 (MQ=255)
 gcagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGtt                                                                     <  1:3705595/71‑1 (MQ=255)
 gcagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGtt                                                                     <  1:3167202/71‑1 (MQ=255)
 gcagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGtt                                                                     <  1:2877452/71‑1 (MQ=255)
 gcagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGtt                                                                     <  1:2741176/71‑1 (MQ=255)
 gcagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGtt                                                                     <  1:1732293/71‑1 (MQ=255)
   agcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTc                                                                   >  1:3347864/1‑71 (MQ=255)
   agcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTc                                                                   <  1:2061197/71‑1 (MQ=255)
        cGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACa                                                               <  1:2567537/70‑1 (MQ=255)
          tGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAtc                                                            <  1:434097/71‑1 (MQ=255)
            ggATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAtctc                                                          <  1:99996/71‑1 (MQ=255)
            ggATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAtctc                                                          >  1:3712463/1‑71 (MQ=255)
             gATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAtctc                                                          <  1:3494664/70‑1 (MQ=255)
               tGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCt                                                        >  1:839963/1‑70 (MQ=255)
                      cATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAg                                                >  1:2147751/1‑71 (MQ=255)
                       aTTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGtt                                                                     <  1:2981549/49‑1 (MQ=255)
                       aTTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGtt                                                                     <  1:3606300/49‑1 (MQ=255)
                       aTTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAtct                                                           <  1:1593429/59‑1 (MQ=255)
                       aTTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAtct                                                           <  1:2049546/59‑1 (MQ=255)
                       aTTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGc                                               >  1:3316682/1‑71 (MQ=255)
                       aTTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGc                                               >  1:2857317/1‑71 (MQ=255)
                       aTTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGc                                               >  1:540008/1‑71 (MQ=255)
                        ttGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTc                                                                   <  1:1595805/50‑1 (MQ=255)
                        ttGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTc                                                                   <  1:1162512/50‑1 (MQ=255)
                        ttGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACa                                                               <  1:706691/54‑1 (MQ=255)
                        ttGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGcc                                              <  1:2147108/71‑1 (MQ=255)
                          ggTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTc                                                                   <  1:2416373/48‑1 (MQ=255)
                          ggTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCcgc                                            >  1:1400550/1‑71 (MQ=255)
                           gTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAtct                                                           <  1:3125336/55‑1 (MQ=255)
                           gTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAtct                                                           <  1:3627548/55‑1 (MQ=255)
                              aGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAt                                                             <  1:1413042/50‑1 (MQ=255)
                              aGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAt                                                             <  1:2569396/50‑1 (MQ=255)
                              aGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAt                                                             <  1:4051714/50‑1 (MQ=255)
                              aGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCt                                                        >  1:4022679/1‑55 (MQ=255)
                              aGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGa                                          >  1:1298777/1‑69 (MQ=255)
                              aGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGAt                                         <  1:945720/70‑1 (MQ=255)
                              aGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATa                                        <  1:2571718/71‑1 (MQ=255)
                                  cATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTg                                                       >  1:1467381/1‑52 (MQ=255)
                                   aTTCCATACCTGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGaa                                   <  1:1788375/71‑1 (MQ=255)
                                   aTTCCATACCGGCACGCTAGCGTATTCCGCCAgtgt                                                                      <  1:645198/36‑1 (MQ=255)
                                   aTTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACt                                      <  1:2951198/68‑1 (MQ=255)
                                   aTTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTg                                     <  1:3825256/69‑1 (MQ=255)
                                   aTTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTg                                     <  1:3489124/69‑1 (MQ=255)
                                    ttCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGAc                                                                >  1:2827083/1‑41 (MQ=255)
                                    ttCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACaa                                                              <  1:1775977/43‑1 (MQ=255)
                                    ttCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCa                                                   >  1:766114/1‑54 (MQ=255)
                                    ttCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGa                                          <  1:3564948/63‑1 (MQ=255)
                                    ttCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGa                                          <  1:3525161/63‑1 (MQ=255)
                                    ttCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGaa                                   >  1:2744934/1‑70 (MQ=255)
                                    ttCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGaa                                   >  1:3602864/1‑70 (MQ=255)
                                      ccATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAt                                                             >  1:1310585/1‑42 (MQ=255)
                                      ccATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGa                                          >  1:587775/1‑61 (MQ=255)
                                        aTACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCg                              <  1:547756/71‑1 (MQ=255)
                                         tACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGAc                                                                <  1:1652988/36‑1 (MQ=255)
                                         tACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACa                                                               <  1:498675/37‑1 (MQ=255)
                                         tACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGAc                                                     <  1:2679199/47‑1 (MQ=255)
                                         tACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATAc                                <  1:3470373/68‑1 (MQ=255)
                                         tACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATAGCGt                             <  1:1016432/71‑1 (MQ=255)
                                         caCCGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGt                             <  1:3084572/70‑1 (MQ=255)
                                          aCCGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGt                             >  1:257623/1‑70 (MQ=255)
                                           ccGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                           >  1:3326762/1‑71 (MQ=255)
                                           ccGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                           >  1:2075023/1‑71 (MQ=255)
                                           ccGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                           >  1:2365229/1‑71 (MQ=255)
                                            cGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCa                           >  1:2466837/1‑70 (MQ=255)
                                            cGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCa                           >  1:651248/1‑70 (MQ=255)
                                              gCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATAc                        <  1:3357700/71‑1 (MQ=255)
                                              gCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATAc                        <  1:2119899/71‑1 (MQ=255)
                                                  gCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                           <  1:4011572/64‑1 (MQ=255)
                                                   cTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGaa                                   >  1:1388457/1‑55 (MQ=255)
                                                   cTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCg                   <  1:1613966/71‑1 (MQ=255)
                                                      gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAc                 >  1:3695835/1‑70 (MQ=255)
                                                       cGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                <  1:2830463/70‑1 (MQ=255)
                                                       cGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                >  1:1433522/1‑70 (MQ=255)
                                                         tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGAt                                         <  1:1795073/43‑1 (MQ=255)
                                                         tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTg                                     <  1:3521037/47‑1 (MQ=255)
                                                         tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTg                                     <  1:1050227/47‑1 (MQ=255)
                                                         tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGt                             <  1:2525052/55‑1 (MQ=255)
                                                         tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                <  1:168086/68‑1 (MQ=255)
                                                         tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAg             <  1:1179099/71‑1 (MQ=255)
                                                          aTTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAg             >  1:1365459/1‑70 (MQ=255)
                                                           ttCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCa           >  1:1606350/1‑71 (MQ=255)
                                                           ttCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCa           >  1:449873/1‑71 (MQ=255)
                                                            tCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGa                                          >  1:2001516/1‑39 (MQ=255)
                                                             ccGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg         >  1:2587752/1‑71 (MQ=255)
                                                              cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                <  1:1292680/63‑1 (MQ=255)
                                                              cGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAAcc         <  1:1021139/70‑2 (MQ=255)
                                                                ccAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCAt                          <  1:1134790/51‑1 (MQ=255)
                                                                  aGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCg                   >  1:2856499/1‑56 (MQ=255)
                                                                  aGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                  >  1:3363395/1‑57 (MQ=255)
                                                                   gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCa   <  1:253905/71‑1 (MQ=255)
                                                                    tgtTTCGATAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAg  >  1:2243959/1‑71 (MQ=255)
                                                                    tgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAg  >  1:2286494/1‑71 (MQ=255)
                                                                    |                                                                      
AGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAG  >  minE/2919280‑2919418

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