Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,269 G→T 14.2% L143I (CTC→ATC) ‡ argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2690GT14.2% 254.3 / 17.2 84L143I (CTC→ATC) ‡argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base T (4/8);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAT  >  minE/2919200‑2919339
                                                                     |                                                                      
gcgcgTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                                                        >  1:3389424/1‑70 (MQ=255)
     tCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATc                                                                  <  1:984147/71‑1 (MQ=255)
         cTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaaga                                                                    <  1:1332913/65‑1 (MQ=255)
           gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:3375487/1‑71 (MQ=255)
           gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:2729199/1‑71 (MQ=255)
           gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:1443836/1‑71 (MQ=255)
            cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:215778/1‑70 (MQ=255)
            cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            <  1:2524156/70‑1 (MQ=255)
              tAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                                         >  1:3820175/1‑71 (MQ=255)
              tAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                                         >  1:1566546/1‑71 (MQ=255)
               aGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCagcagc                                                        >  1:345191/1‑71 (MQ=255)
                gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                                       <  1:488704/71‑1 (MQ=255)
                gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                                       <  1:1845784/71‑1 (MQ=255)
                gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                                       <  1:678766/71‑1 (MQ=255)
                gACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                                       >  1:2324041/1‑71 (MQ=255)
                 aCCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCcagca                                                          <  1:1536491/67‑1 (MQ=255)
                  ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGc                                                     <  1:417966/71‑1 (MQ=255)
                   cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                                                      >  1:3417191/1‑53 (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:1516533/71‑1 (MQ=255)
                           tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                            <  1:291035/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:1989517/70‑3 (MQ=255)
                                   aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                                     <  1:3478967/70‑3 (MQ=255)
                                    aCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTgg                                                     >  1:705940/1‑52 (MQ=255)
                                     cGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGAt                                                                   >  1:1315962/1‑38 (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:356765/1‑71 (MQ=255)
                                           gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                            >  1:1081914/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:2438395/71‑1 (MQ=255)
                                                   aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                    <  1:1589273/71‑1 (MQ=255)
                                                    gtGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                                       <  1:2163477/34‑1 (MQ=16)
                                                    gtGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                                       <  1:2551466/34‑1 (MQ=16)
                                                    aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATa                   <  1:738670/71‑1 (MQ=255)
                                                     tGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                          >  1:278211/1‑47 (MQ=16)
                                                      gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                                             >  1:3463126/1‑43 (MQ=255)
                                                      gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAcc                 >  1:1154139/1‑71 (MQ=255)
                                                      gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAcc                 >  1:2131374/1‑71 (MQ=255)
                                                        tCCTGCATGGTGAGAAGATCCGCCAGCAGGTGCGTGGGATGGAACTCATTGGTCAgg                             >  1:1842506/1‑57 (MQ=255)
                                                        tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                              <  1:4085/56‑1 (MQ=255)
                                                        tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                              <  1:3008816/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:254901/1‑71 (MQ=255)
                                                            gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg           >  1:1567680/1‑71 (MQ=255)
                                                             cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc          <  1:1445464/71‑1 (MQ=255)
                                                             cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc          <  1:985564/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:2488977/1‑71 (MQ=255)
                                                                     gaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAt  <  1:1442393/71‑1 (MQ=255)
                                                                     |                                                                      
GCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAT  >  minE/2919200‑2919339

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: