Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,174 G→C 7.2% L174L (CTC→CTG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1740GC 7.2% 209.1 / 3.9 68L174L (CTC→CTGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (25/38);  new base C (3/2);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGC  >  minE/2919106‑2919244
                                                                    |                                                                      
gcAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGACGCAGCTTCGAg                                                                      >  1:1889231/1‑71 (MQ=255)
 cAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGc                                                                     <  1:3771981/71‑1 (MQ=255)
  aGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGc                                                                     <  1:118685/70‑1 (MQ=255)
   tcTTCCGGCCAGCAGGCTTTCGGGGCCAACAGGCGCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCAt                                                                   >  1:800963/2‑71 (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:1724167/71‑1 (MQ=255)
                 cGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgt                                                     <  1:2825293/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)
                                                    tcaGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCAtc                  >  1:3373137/3‑70 (MQ=255)
                                                      aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgcgc                                           >  1:2090193/1‑44 (MQ=255)
                                                      aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                >  1:73158/1‑71 (MQ=255)
                                                      aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                >  1:2368877/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:1296743/71‑1 (MQ=255)
                                                           cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGa           <  1:2055691/71‑1 (MQ=255)
                                                             cAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                 >  1:632986/1‑47 (MQ=255)
                                                             cAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                 >  1:1252383/1‑47 (MQ=255)
                                                               gCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                 <  1:3639030/45‑1 (MQ=255)
                                                               gCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGc       >  1:866859/1‑71 (MQ=255)
                                                                 ttCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                               >  1:1838680/1‑45 (MQ=255)
                                                                  tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                                       >  1:2707584/1‑36 (MQ=255)
                                                                  tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGc                                  >  1:3163876/1‑41 (MQ=255)
                                                                  tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                 >  1:444015/1‑42 (MQ=255)
                                                                  tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                 >  1:2039151/1‑42 (MQ=255)
                                                                  tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                 >  1:3159334/1‑42 (MQ=255)
                                                                  tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt    >  1:2638871/1‑71 (MQ=255)
                                                                   cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt    <  1:3383346/70‑1 (MQ=255)
                                                                    gAGCATCGAATTGCCCATGTTGTTAGGCGTGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc  <  1:1053910/71‑1 (MQ=255)
                                                                    gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc  >  1:3029943/1‑71 (MQ=255)
                                                                    |                                                                      
GCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGC  >  minE/2919106‑2919244

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: