Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,399 G→C 7.8% R99R (CGC→CGG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,3990GC 7.8% 262.5 / 9.8 89R99R (CGC→CGGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (38/44);  new base C (5/2);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGA  >  minE/2919329‑2919468
                                                                      |                                                                     
cGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATg                                                                       <  1:2342932/71‑1 (MQ=255)
 gCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGc                                                                      >  1:1302964/1‑71 (MQ=255)
   tAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGc                                                                      <  1:515452/69‑1 (MQ=255)
   tAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                                    >  1:2230889/1‑71 (MQ=255)
     gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                                    >  1:1440037/1‑69 (MQ=255)
     gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                                  <  1:1156042/71‑1 (MQ=255)
     gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAc                                                                   >  1:54930/1‑70 (MQ=255)
     gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAc                                                                   >  1:2425850/1‑70 (MQ=255)
     gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAc                                                                   >  1:3137877/1‑70 (MQ=255)
     gCGTATTCCGCCAGCGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                                    >  1:3579612/1‑69 (MQ=255)
           tgcGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAAc                                                            >  1:3176586/3‑71 (MQ=255)
            gcGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAg                                                               >  1:2482269/2‑67 (MQ=255)
             cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                                  <  1:3164253/63‑1 (MQ=255)
             cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                                           <  1:1019544/70‑1 (MQ=255)
             cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                                           <  1:1303261/70‑1 (MQ=255)
             cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                                           <  1:1708370/70‑1 (MQ=255)
              gCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcgcg                                                         <  1:2143831/71‑1 (MQ=255)
                cAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCggg                                                       <  1:1226227/71‑1 (MQ=255)
                  gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCa                                                     <  1:3019389/71‑1 (MQ=255)
                  gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCa                                                     <  1:2081044/71‑1 (MQ=255)
                  gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCATCATACATGCGACCAAGCACGCGGGCa                                                     <  1:1911017/71‑1 (MQ=255)
                  gcgtTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAACCCgcgcg                                                     <  1:778990/69‑2 (MQ=255)
                    gtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgt                                                   >  1:576050/1‑71 (MQ=255)
                     tttCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtg                                                  <  1:1409972/71‑1 (MQ=255)
                      ttCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtg                                                  >  1:157666/1‑70 (MQ=255)
                       tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgt                                                   >  1:2757035/1‑68 (MQ=255)
                        cGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGc                                                                      <  1:701256/48‑1 (MQ=255)
                        cGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                                    >  1:410678/1‑50 (MQ=255)
                         gACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                                  <  1:3801225/51‑1 (MQ=255)
                          aCAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAc                                                                   >  1:1237543/1‑49 (MQ=255)
                          aCAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCtt                                              >  1:3666389/1‑70 (MQ=255)
                            aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAc                                                            >  1:257113/1‑54 (MQ=255)
                            aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCttt                                             >  1:937260/1‑69 (MQ=255)
                            aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                           <  1:492581/71‑1 (MQ=255)
                            aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                           >  1:3046045/1‑71 (MQ=255)
                            aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                           <  1:312579/71‑1 (MQ=255)
                              tctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                         >  1:2754572/1‑71 (MQ=255)
                                tcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGc                                                      >  1:3645628/1‑56 (MQ=255)
                                tcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                                       >  1:3655131/1‑71 (MQ=255)
                                tcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                                       >  1:1711836/1‑71 (MQ=255)
                                 ccTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAc                                      >  1:1750873/1‑71 (MQ=255)
                                  cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                         <  1:2676158/67‑1 (MQ=255)
                                  cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCg                                        <  1:3462241/68‑1 (MQ=255)
                                  cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                                     >  1:2337930/1‑71 (MQ=255)
                                  cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAc                                      >  1:1923432/1‑70 (MQ=255)
                                   tGACCATAGCCGCGATACTGAATACCGTCATACATGc                                                                      <  1:1708387/37‑1 (MQ=255)
                                   tGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCt                                               <  1:1139500/60‑1 (MQ=255)
                                     aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCt                                               >  1:1550077/1‑58 (MQ=255)
                                     aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActc                                    <  1:342752/69‑1 (MQ=255)
                                     aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                  >  1:1105342/1‑71 (MQ=255)
                                     aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                  >  1:3646189/1‑71 (MQ=255)
                                      ccATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCttt                                 <  1:1678526/71‑1 (MQ=255)
                                       cATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCtt                                              <  1:282461/57‑1 (MQ=255)
                                        aTAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                               <  1:1988056/71‑1 (MQ=255)
                                         tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                               <  1:3021845/70‑1 (MQ=255)
                                         tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                               >  1:1231583/1‑70 (MQ=255)
                                         tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                              <  1:3527365/71‑1 (MQ=255)
                                         tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                              <  1:1195531/71‑1 (MQ=255)
                                         tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTCTCTTTAATCGACTCTTTAt                               >  1:2264886/1‑70 (MQ=255)
                                          aGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTa                                            <  1:2481641/56‑1 (MQ=255)
                                            ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGa                             >  1:2516378/1‑69 (MQ=255)
                                            ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc                            >  1:287660/1‑70 (MQ=255)
                                            ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc                            >  1:3877993/1‑70 (MQ=255)
                                             cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAcc                           >  1:3140073/1‑70 (MQ=255)
                                               cgATACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAAt                        >  1:3343308/1‑71 (MQ=255)
                                               cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt                        <  1:1501430/71‑1 (MQ=255)
                                               cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt                        <  1:1072865/71‑1 (MQ=255)
                                                 aTACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt                      <  1:3140774/71‑1 (MQ=255)
                                                  tACTTAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt                      <  1:3182861/70‑1 (MQ=255)
                                                  tACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTg                     >  1:1151001/1‑71 (MQ=255)
                                                      gAAAACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAc                                      <  1:44526/50‑1 (MQ=255)
                                                       aaaaCCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                 <  1:2421608/67‑1 (MQ=255)
                                                       aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTgg                    >  1:1196747/1‑67 (MQ=255)
                                                       aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                 <  1:3400262/70‑1 (MQ=255)
                                                       aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                 <  1:1588042/70‑1 (MQ=255)
                                                       aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                 <  1:3576235/70‑1 (MQ=255)
                                                        aTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                               >  1:3646598/1‑55 (MQ=255)
                                                         tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                              <  1:1969164/55‑1 (MQ=255)
                                                         tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg              >  1:2875305/1‑71 (MQ=255)
                                                         tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg              <  1:1014431/71‑1 (MQ=255)
                                                           ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCt            <  1:530444/71‑1 (MQ=255)
                                                            cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt           >  1:390611/1‑71 (MQ=255)
                                                            cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt           <  1:3216134/71‑1 (MQ=255)
                                                              tCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTgg         <  1:3916074/71‑1 (MQ=255)
                                                              tCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCgg         >  1:3047201/1‑71 (MQ=255)
                                                                 tACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc      <  1:3711769/71‑1 (MQ=255)
                                                                   cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                  <  1:1304227/41‑1 (MQ=255)
                                                                   cATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGcccta    <  1:1101366/71‑3 (MQ=255)
                                                                    aTGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc      >  1:3325441/1‑68 (MQ=255)
                                                                    aTGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc      >  1:2753508/1‑68 (MQ=255)
                                                                      gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTa                                >  1:64676/1‑40 (MQ=255)
                                                                      gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga  >  1:1671216/1‑70 (MQ=255)
                                                                      gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga  >  1:286668/1‑70 (MQ=255)
                                                                      gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga  >  1:2474233/1‑70 (MQ=255)
                                                                      |                                                                     
CGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGA  >  minE/2919329‑2919468

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: