Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,159 A→C 6.5% L179L (CTT→CTG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1590AC 6.5% 162.3 / 6.5 62L179L (CTT→CTGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (24/34);  new base C (2/2);  total (26/36)
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.46e-01

GGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAT  >  minE/2919089‑2919227
                                                                      |                                                                    
ggCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTa                                                                      <  1:1871865/71‑1 (MQ=255)
 gCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTaa                                                                     <  1:1956572/71‑1 (MQ=255)
      ccGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCGAATCCAGACCGGTAAGcgcc                                                                >  1:3204820/1‑71 (MQ=255)
      ccGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgcc                                                                >  1:3102944/1‑71 (MQ=255)
      ccGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgcc                                                                >  1:2933541/1‑71 (MQ=255)
      ccGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgcc                                                                >  1:1547302/1‑71 (MQ=255)
        gTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgccg                                                               >  1:3976433/1‑70 (MQ=255)
         tAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgccgc                                                              <  1:1639513/70‑1 (MQ=255)
            ccAGCGCAGCTTCCGGCCAGCACGCTTGTTGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCt                                                          >  1:1133181/1‑71 (MQ=255)
            ccAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCt                                                          >  1:3464485/1‑71 (MQ=255)
             cAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCtt                                                         >  1:2574056/1‑71 (MQ=255)
             cAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCtt                                                         >  1:2202388/1‑71 (MQ=255)
              aGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTc                                                        <  1:415287/71‑1 (MQ=255)
              aGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTc                                                        <  1:2189107/71‑1 (MQ=255)
                cgcAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTc                                                        >  1:871257/1‑69 (MQ=255)
                 gcAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGACGCAGCTTCGAg                                                     >  1:1889231/1‑71 (MQ=255)
                  cAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTc                                                        >  1:2031368/1‑67 (MQ=255)
                  cAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGc                                                    <  1:3771981/71‑1 (MQ=255)
                   aGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGc                                                    <  1:118685/70‑1 (MQ=255)
                    tcTTCCGGCCAGCAGGCTTTCGGGGCCAACAGGCGCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCAt                                                  >  1:800963/2‑71 (MQ=255)
                     cTTCCGGCCAGCAGGCTTTCGGGGCCAACAGGCGCAGATCCAGCCCGGTCAGcgccgccgc                                                           >  1:1686392/1‑61 (MQ=255)
                     cTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTc                                                        >  1:396156/1‑64 (MQ=255)
                       tCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCtt                                                         >  1:2828900/1‑61 (MQ=255)
                       tCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                                >  1:1909483/1‑70 (MQ=255)
                        ccGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGa                                               <  1:3970053/70‑1 (MQ=255)
                          ggCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAAtt                                            <  1:2481428/71‑1 (MQ=255)
                             cAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                         <  1:616015/71‑1 (MQ=255)
                                  cGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgt                                    <  1:2825293/71‑1 (MQ=255)
                                  cGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgt                                    <  1:1724167/71‑1 (MQ=255)
                                     ttGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                                <  1:586403/56‑1 (MQ=255)
                                       gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                                 <  1:1649156/69‑1 (MQ=255)
                                       gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                                 <  1:672418/69‑1 (MQ=255)
                                       gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATgttgtt                                <  1:3894669/70‑1 (MQ=255)
                                        tgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAc                              <  1:3365917/71‑1 (MQ=255)
                                          gcgcGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgc                            >  1:4046496/1‑71 (MQ=255)
                                              gACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGt                        <  1:1252158/71‑1 (MQ=255)
                                                ccAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                      <  1:865859/71‑1 (MQ=255)
                                                acAGGCGCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATcg                      <  1:1698637/70‑2 (MQ=255)
                                                 cAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGa                                                      >  1:740496/1‑38 (MQ=255)
                                                 cAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                                <  1:283685/44‑1 (MQ=255)
                                                 cAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                                <  1:3397154/44‑1 (MQ=255)
                                                  aGACGCAAATCCAGACCGGTAAGCGCCGCATCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCAcc                    >  1:636375/1‑71 (MQ=255)
                                                    aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg                  <  1:2070343/71‑1 (MQ=255)
                                                    aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg                  <  1:3903685/71‑1 (MQ=255)
                                                      gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGa                                               <  1:2652221/40‑1 (MQ=255)
                                                      gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGa                                               <  1:278062/40‑1 (MQ=255)
                                                      gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTg                                           >  1:2229634/1‑44 (MQ=255)
                                                      gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTc                       <  1:2806903/64‑1 (MQ=255)
                                                      gCAAATCCAGACCGGTAAGCGCCGCACCTTCGAGCATCGAATTGCCCATGTTGTTAc                              >  1:3790474/1‑57 (MQ=255)
                                                       cAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                                <  1:2841360/38‑1 (MQ=255)
                                                       cAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGTGCGTCACCTGCa                >  1:3907997/1‑70 (MQ=255)
                                                       cAAATCCAGACCGGTAAGAGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                >  1:1794130/1‑70 (MQ=255)
                                                        aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGaa                                              >  1:3457089/1‑39 (MQ=255)
                                                        aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa              <  1:3379729/71‑1 (MQ=255)
                                                         aaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcg                             <  1:3282971/55‑1 (MQ=255)
                                                          aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                         <  1:1036567/42‑1 (MQ=255)
                                                           tCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                <  1:1816826/66‑1 (MQ=255)
                                                           tCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                <  1:1135382/66‑1 (MQ=255)
                                                           tCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                <  1:2544179/66‑1 (MQ=255)
                                                            ccAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAt                                      <  1:968806/43‑1 (MQ=255)
                                                             cAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGAcc          <  1:2526445/70‑1 (MQ=255)
                                                                  cGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTc    <  1:714916/71‑1 (MQ=255)
                                                                    gTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAt  >  1:2721947/1‑71 (MQ=255)
                                                                    gTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAt  <  1:3430471/71‑1 (MQ=255)
                                                                      |                                                                    
GGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAT  >  minE/2919089‑2919227

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: