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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4110GC 5.8% 266.1 / 4.6 86R95R (CGC→CGGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (40/41);  new base C (3/2);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAG  >  minE/2919342‑2919481
                                                                     |                                                                      
cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                                                        <  1:1019544/70‑1 (MQ=255)
cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                                                        <  1:1708370/70‑1 (MQ=255)
cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                                                        <  1:1303261/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)
             aCAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCtt                                                           >  1:3666389/1‑70 (MQ=255)
               aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCttt                                                          >  1:937260/1‑69 (MQ=255)
               aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                        <  1:492581/71‑1 (MQ=255)
               aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                        <  1:312579/71‑1 (MQ=255)
               aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                        >  1:3046045/1‑71 (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)
                      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:1195531/71‑1 (MQ=255)
                            tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                                           <  1:3527365/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:3877993/1‑70 (MQ=255)
                               ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc                                         >  1:287660/1‑70 (MQ=255)
                                cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAcc                                        >  1:3140073/1‑70 (MQ=255)
                                  cgATACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAAt                                     >  1:3343308/1‑71 (MQ=255)
                                  cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt                                     <  1:1072865/71‑1 (MQ=255)
                                  cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt                                     <  1:1501430/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:3576235/70‑1 (MQ=255)
                                          aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                              <  1:3400262/70‑1 (MQ=255)
                                          aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                              <  1:1588042/70‑1 (MQ=255)
                                           aTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                            >  1:3646598/1‑55 (MQ=255)
                                            tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                                           <  1:1969164/55‑1 (MQ=255)
                                            tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg                           <  1:1014431/71‑1 (MQ=255)
                                            tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg                           >  1:2875305/1‑71 (MQ=255)
                                              ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCt                         <  1:530444/71‑1 (MQ=255)
                                               cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt                        <  1:3216134/71‑1 (MQ=255)
                                               cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt                        >  1:390611/1‑71 (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:2753508/1‑68 (MQ=255)
                                                       aTGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc                   >  1:3325441/1‑68 (MQ=255)
                                                         gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTa                                             >  1:64676/1‑40 (MQ=255)
                                                         gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga               >  1:286668/1‑70 (MQ=255)
                                                         gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga               >  1:1671216/1‑70 (MQ=255)
                                                         gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga               >  1:2474233/1‑70 (MQ=255)
                                                          cGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGAt              >  1:628717/1‑70 (MQ=255)
                                                           gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                            <  1:3747826/39‑1 (MQ=255)
                                                           gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa            <  1:3385590/71‑1 (MQ=255)
                                                           gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa            <  1:1545066/71‑1 (MQ=255)
                                                             ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGcc                            <  1:1916688/53‑1 (MQ=255)
                                                             ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg           >  1:3744153/1‑70 (MQ=255)
                                                             ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg           >  1:1435759/1‑70 (MQ=255)
                                                             ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg           >  1:2606645/1‑70 (MQ=255)
                                                             ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg           >  1:3253438/1‑70 (MQ=255)
                                                              cAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTa         <  1:1338284/71‑1 (MQ=255)
                                                              cAAGCACGCGGGCAGTGTCTTTAATAGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGt          >  1:3808790/1‑70 (MQ=255)
                                                                aGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAc       >  1:2019871/1‑71 (MQ=255)
                                                                aGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAc       <  1:800302/71‑1 (MQ=255)
                                                                  cACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc     <  1:198903/71‑1 (MQ=255)
                                                                  cACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcg      >  1:1024561/1‑70 (MQ=255)
                                                                    cgcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATc                                    <  1:3360299/38‑1 (MQ=255)
                                                                    cgcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa   <  1:1212140/71‑1 (MQ=255)
                                                                     gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg  >  1:3674753/1‑71 (MQ=255)
                                                                     gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg  >  1:2676850/1‑71 (MQ=255)
                                                                     |                                                                      
CGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAG  >  minE/2919342‑2919481

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: