Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,288 C→G 8.9% T136T (ACG→ACC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2880CG 8.9% 224.0 / 8.8 78T136T (ACG→ACCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (32/39);  new base G (2/5);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAA  >  minE/2919218‑2919358
                                                                      |                                                                      
ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGc                                                                        <  1:417966/71‑1 (MQ=255)
 cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg                                                                       >  1:898756/1‑71 (MQ=255)
      tcctttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa                                                                  <  1:1672665/68‑1 (MQ=255)
      tCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa                                                                  <  1:1495109/71‑1 (MQ=255)
       cATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt                                                                 >  1:3709877/1‑71 (MQ=255)
        atctCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTg                                                                <  1:2544862/68‑1 (MQ=255)
        aTTTCGTTGAACGCTTTGCCTGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                                                                <  1:1517941/71‑1 (MQ=255)
         tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGTGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                               <  1:2800/71‑1 (MQ=255)
         tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                               <  1:291035/71‑1 (MQ=255)
         tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                               <  1:1516533/71‑1 (MQ=255)
         tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGAATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                               <  1:2940543/71‑1 (MQ=255)
            cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                                            >  1:464517/1‑71 (MQ=255)
            cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                                            >  1:3500525/1‑71 (MQ=255)
             gTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                                             <  1:3462692/69‑1 (MQ=255)
              ttGAACGCTTTGCCGGGCAAATGGTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                               >  1:16463/1‑66 (MQ=255)
              ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                                          <  1:1671315/71‑1 (MQ=255)
              ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                                          <  1:1014133/71‑1 (MQ=255)
               tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGTATGGAACTCa                                                         >  1:242989/1‑71 (MQ=255)
               tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                                          <  1:856089/70‑1 (MQ=255)
               tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                                                         >  1:632586/1‑71 (MQ=255)
                gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATTGAACTCa                                                         <  1:2928175/70‑1 (MQ=255)
                 aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                               <  1:3392923/63‑1 (MQ=255)
                 aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                        <  1:935183/70‑1 (MQ=255)
                 aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                                                        <  1:3478967/70‑3 (MQ=255)
                 aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                                                        <  1:1989517/70‑3 (MQ=255)
                    gCTTTGCCGGGCAATTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                    <  1:2727425/71‑1 (MQ=255)
                    gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                    <  1:2405156/71‑1 (MQ=255)
                    gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                    <  1:56103/71‑1 (MQ=255)
                      tttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                    >  1:25113/1‑69 (MQ=255)
                      cttGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTCGAACTCGTTGGt                                                    >  1:260038/2‑69 (MQ=255)
                       ttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                                                 >  1:1152492/1‑71 (MQ=255)
                        tGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg                                                >  1:983456/1‑71 (MQ=255)
                         gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                               >  1:1081914/1‑71 (MQ=255)
                         gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                               >  1:356765/1‑71 (MQ=255)
                           cGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGcc                                              <  1:1194635/70‑1 (MQ=255)
                            gggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAt                                            <  1:2248350/71‑1 (MQ=255)
                             ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt                                           <  1:3517193/71‑1 (MQ=255)
                             ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt                                           <  1:3875150/71‑1 (MQ=255)
                              gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc                                          >  1:2506764/1‑71 (MQ=255)
                               cAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc                                          >  1:3647491/1‑70 (MQ=255)
                               cAAATGCTCCTGCATGGTGAGAAGATCCGCAAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc                                         <  1:2330912/71‑1 (MQ=255)
                                aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc                                         >  1:1394035/1‑70 (MQ=255)
                                aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa                                        >  1:3419366/1‑71 (MQ=255)
                                 aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                                                              >  1:2164908/1‑48 (MQ=255)
                                 aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                                                              >  1:301951/1‑48 (MQ=255)
                                 aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAtt                                                       >  1:2429364/1‑55 (MQ=255)
                                 aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                                       <  1:1589273/71‑1 (MQ=255)
                                 aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                                       <  1:2438395/71‑1 (MQ=255)
                                  aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATa                                      <  1:738670/71‑1 (MQ=255)
                                   tGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                                             >  1:278211/1‑47 (MQ=16)
                                    gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                                                                >  1:3463126/1‑43 (MQ=255)
                                    gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAcc                                    >  1:2131374/1‑71 (MQ=255)
                                    gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAcc                                    >  1:1154139/1‑71 (MQ=255)
                                      tCCTGCATGGTGAGAAGATCCGCCAGCAGGTGCGTGGGATGGAACTCATTGGTCAgg                                                >  1:1842506/1‑57 (MQ=255)
                                      tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                                                 <  1:3008816/56‑1 (MQ=255)
                                      tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                                                 <  1:4085/56‑1 (MQ=255)
                                       ccTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATTGAACTCATTGGTCAg                                                 <  1:3099690/55‑1 (MQ=255)
                                        cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCa                                >  1:2599004/1‑71 (MQ=255)
                                        cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGGCAGGCCATTCCATACCgg                                  <  1:2504979/69‑1 (MQ=255)
                                        cTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCACACCGGc                                 <  1:2359793/70‑1 (MQ=255)
                                          gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg                              >  1:1567680/1‑71 (MQ=255)
                                          gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg                              >  1:254901/1‑71 (MQ=255)
                                           cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                             <  1:985564/71‑1 (MQ=255)
                                           cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                             <  1:1445464/71‑1 (MQ=255)
                                            aTGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                             >  1:1852288/1‑70 (MQ=255)
                                                tGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCg                        <  1:236448/71‑1 (MQ=255)
                                                  agaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                        <  1:2142807/37‑1 (MQ=255)
                                                  agaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTa                      <  1:3870150/71‑1 (MQ=255)
                                                   gaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAt                     <  1:1442393/71‑1 (MQ=255)
                                                   gaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAt                     >  1:2488977/1‑71 (MQ=255)
                                                     agaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTc                   >  1:2220404/1‑71 (MQ=255)
                                                       aTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTc                                                   <  1:627250/37‑1 (MQ=255)
                                                        tCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGc                >  1:2926612/1‑71 (MQ=255)
                                                          cGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCa              <  1:2358888/71‑1 (MQ=255)
                                                             cagcagCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGcc               <  1:1370370/67‑1 (MQ=255)
                                                                cagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGttt        >  1:2596451/1‑71 (MQ=255)
                                                                cagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGttt        >  1:2651682/1‑71 (MQ=255)
                                                                 agcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTc       <  1:568277/71‑1 (MQ=255)
                                                                   cTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGttt        <  1:2716573/68‑1 (MQ=255)
                                                                      cGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACaa  >  1:3275705/1‑71 (MQ=255)
                                                                      |                                                                      
CCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAA  >  minE/2919218‑2919358

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: