Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,153 A→C 6.3% G181G (GGT→GGG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1530AC 6.3% 152.9 / 4.0 62G181G (GGT→GGGargIornithine 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.73e-01

GCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACC  >  minE/2919085‑2919219
                                                                    |                                                                  
gcgcGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGAcc                                                                  >  1:2818034/1‑71 (MQ=255)
    ggCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTa                                                              <  1:1871865/71‑1 (MQ=255)
     gCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTaa                                                             <  1:1956572/71‑1 (MQ=255)
          ccGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCGAATCCAGACCGGTAAGcgcc                                                        >  1:3204820/1‑71 (MQ=255)
          ccGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgcc                                                        >  1:2933541/1‑71 (MQ=255)
          ccGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgcc                                                        >  1:3102944/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:2202388/1‑71 (MQ=255)
                 cAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCtt                                                 >  1:2574056/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:3397154/44‑1 (MQ=255)
                                                     cAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                        <  1:283685/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)
                                                                    |                                                                  
GCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACC  >  minE/2919085‑2919219

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: