Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,273 A→G 14.3% D141D (GAT→GAC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2730AG14.3% 210.3 / 22.8 84D141D (GAT→GACargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base G (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 = 1.00e+00

TCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTC  >  minE/2919205‑2919341
                                                                    |                                                                    
tCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATc                                                                    <  1:984147/71‑1 (MQ=255)
    cTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaaga                                                                      <  1:1332913/65‑1 (MQ=255)
      gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                              >  1:2729199/1‑71 (MQ=255)
      gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                              >  1:3375487/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:1566546/1‑71 (MQ=255)
         tAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                                           >  1:3820175/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)
              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:3500525/1‑71 (MQ=255)
                         cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                           >  1:464517/1‑71 (MQ=255)
                          gTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                            <  1:3462692/69‑1 (MQ=255)
                           ttGAACGCTTTGCCGGGCAAATGGTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                              >  1:16463/1‑66 (MQ=255)
                           ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                         <  1:1014133/71‑1 (MQ=255)
                           ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                         <  1:1671315/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:3875150/71‑1 (MQ=255)
                                          ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt                          <  1:3517193/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)
                                               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: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: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: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)
                                                                    |                                                                    
TCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTC  >  minE/2919205‑2919341

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: