Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,168 A→C 5.9% A176A (GCT→GCG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1680AC 5.9% 177.8 / 5.0 68A176A (GCT→GCGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (26/38);  new base C (2/2);  total (28/40)
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.55e-01

CCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGA  >  minE/2919101‑2919235
                                                                   |                                                                   
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:672418/69‑1 (MQ=255)
                           gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                                         <  1:1649156/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:3903685/71‑1 (MQ=255)
                                        aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg                          <  1:2070343/71‑1 (MQ=255)
                                          gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGa                                                       <  1:278062/40‑1 (MQ=255)
                                          gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGa                                                       <  1:2652221/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:1135382/66‑1 (MQ=255)
                                               tCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                        <  1:2544179/66‑1 (MQ=255)
                                               tCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                        <  1:1816826/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:3430471/71‑1 (MQ=255)
                                                        gTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAt          >  1:2721947/1‑71 (MQ=255)
                                                         tcaGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCAtc         >  1:3373137/3‑70 (MQ=255)
                                                           aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgcgc                                  >  1:2090193/1‑44 (MQ=255)
                                                           aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc       >  1:2368877/1‑71 (MQ=255)
                                                           aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc       >  1:73158/1‑71 (MQ=255)
                                                            gcgcCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCa                 <  1:1950686/60‑1 (MQ=255)
                                                            gcgcCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCa                 <  1:1818742/60‑1 (MQ=255)
                                                             cgccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt             >  1:3649841/1‑63 (MQ=255)
                                                               ccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGtt    >  1:929025/1‑70 (MQ=255)
                                                                cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCAcc                            <  1:3324485/45‑1 (MQ=255)
                                                                cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGa  <  1:2055691/71‑1 (MQ=255)
                                                                cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGa  <  1:1296743/71‑1 (MQ=255)
                                                                  cAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                        >  1:632986/1‑47 (MQ=255)
                                                                  cAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                        >  1:1252383/1‑47 (MQ=255)
                                                                   |                                                                   
CCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGA  >  minE/2919101‑2919235

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: