Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,237 C→T 47.0% A153A (GCG→GCA argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2370CT47.0% 13.8 / 86.3 102A153A (GCG→GCAargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (38/16);  new base T (20/28);  total (58/44)
Fisher's exact test for biased strand distribution p-value = 4.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

Only 100 of 102 total aligned reads displayed.

GCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCA  >  minE/2919176‑2919303
                                                             |                                                                  
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:1072827/68‑1 (MQ=255)
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:1939703/68‑1 (MQ=255)
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:1879251/68‑1 (MQ=255)
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:2081378/68‑1 (MQ=255)
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:1644672/68‑1 (MQ=255)
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:1644369/68‑1 (MQ=255)
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:482966/68‑1 (MQ=255)
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:525201/68‑1 (MQ=255)
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:734643/68‑1 (MQ=255)
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:874462/68‑1 (MQ=255)
     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCggg                                                         >  1:1548677/1‑68 (MQ=255)
     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCggg                                                         >  1:805939/1‑68 (MQ=255)
     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCggg                                                         >  1:1155254/1‑68 (MQ=255)
       aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                       <  1:2069804/68‑1 (MQ=255)
        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                     >  1:2386815/1‑69 (MQ=255)
        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                     >  1:1090229/1‑69 (MQ=255)
        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                     >  1:2104895/1‑69 (MQ=255)
        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                     >  1:2297010/1‑69 (MQ=255)
        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                     >  1:2299895/1‑69 (MQ=255)
        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                     >  1:1647396/1‑69 (MQ=255)
        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                     >  1:2719747/1‑69 (MQ=255)
        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                     >  1:3251603/1‑69 (MQ=255)
        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                     >  1:359157/1‑69 (MQ=255)
        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaa                                                      >  1:2282201/1‑68 (MQ=255)
            ccATGTTGTTACGCGCGTTACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGc                                                  <  1:1594738/68‑1 (MQ=255)
            ccATGTTGTTACGCGCGTTACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGc                                                  <  1:158994/68‑1 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                >  1:704842/1‑69 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                >  1:2533861/1‑69 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                >  1:2328560/1‑69 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                >  1:3184669/1‑69 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                >  1:3381440/1‑69 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                >  1:38789/1‑69 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                >  1:957977/1‑69 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                >  1:1538937/1‑69 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                >  1:1138392/1‑69 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                >  1:1736846/1‑69 (MQ=255)
             cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGCGCAAATGCTc                                                >  1:2671091/1‑69 (MQ=255)
                 ttgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGc                                            >  1:1889113/1‑69 (MQ=255)
                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCa                                                       <  1:156312/57‑1 (MQ=255)
                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCa                                                       <  1:1574607/57‑1 (MQ=255)
                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCa                                                       <  1:2651936/57‑1 (MQ=255)
                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCa                                                       <  1:2647841/57‑1 (MQ=255)
                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCa                                                       <  1:2560099/57‑1 (MQ=255)
                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCa                                                       <  1:1593198/57‑1 (MQ=255)
                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCa                                                       <  1:1621420/57‑1 (MQ=255)
                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGTAAAt                                                    <  1:855115/52‑1 (MQ=255)
                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCAAAt                                                    <  1:2256612/52‑1 (MQ=255)
                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCAAAt                                                    <  1:1830965/52‑1 (MQ=255)
                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCAAAt                                                    <  1:2018641/52‑1 (MQ=255)
                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCAAAt                                                    <  1:3252431/52‑1 (MQ=255)
                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCAAAt                                                    <  1:3267956/52‑1 (MQ=255)
                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCAAAt                                                    <  1:3324886/52‑1 (MQ=255)
                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGCCGGGCAAAt                                                    <  1:223443/52‑1 (MQ=255)
                          gcgTCACCTGCATAGACCAGCGTCATTTCGATGAATGCTTTGCCGGGCAAAt                                                    <  1:1587465/52‑1 (MQ=255)
                            gTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaag                                 <  1:1847192/69‑1 (MQ=255)
                            gTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaag                                 <  1:576188/69‑1 (MQ=255)
                            gTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaag                                 <  1:2032714/69‑1 (MQ=255)
                               aCCTGCATAGACCAGCGTCATTTCGTTGAATGCTTTGACGGGCAAAt                                                    <  1:827500/47‑1 (MQ=255)
                                 cTGCATAGACCAGCGTCATTTCGTTGAATGCTTTg                                                              <  1:3078532/35‑1 (MQ=255)
                                  tGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                             >  1:778813/1‑67 (MQ=255)
                                  tGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                             >  1:151643/1‑67 (MQ=255)
                                          ccAGCGTCATTTCGTTGAACGCTTTGCCTGGCAAATGCTCCTGCAt                                          >  1:16580/1‑46 (MQ=255)
                                          ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCAt                                          >  1:1552741/1‑46 (MQ=255)
                                          ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCAt                                          >  1:684310/1‑46 (MQ=255)
                                          ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCAt                                          >  1:2705489/1‑46 (MQ=255)
                                          ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCAt                                          >  1:1790155/1‑46 (MQ=255)
                                          ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCAt                                          >  1:770269/1‑46 (MQ=255)
                                           cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                  <  1:702588/69‑1 (MQ=255)
                                           cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                  <  1:929533/69‑1 (MQ=255)
                                           cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                  <  1:3318113/69‑1 (MQ=255)
                                            aGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                  <  1:2970336/68‑1 (MQ=255)
                                            aGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                  <  1:3025486/68‑1 (MQ=255)
                                            aGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                  <  1:1159634/68‑1 (MQ=255)
                                                  aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa           >  1:1397775/1‑69 (MQ=255)
                                                  aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa           >  1:387736/1‑69 (MQ=255)
                                                  aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa           >  1:1879786/1‑69 (MQ=255)
                                                  aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa           >  1:2547192/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTgaggag                                 >  1:1590967/1‑39 (MQ=37)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:2644717/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:1033824/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:2132764/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:1079628/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:21368/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:1743030/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:1686910/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:241786/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:773252/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:725419/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:2836573/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:671556/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:569326/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:2918987/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:3196566/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:3360658/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:1557457/1‑69 (MQ=255)
                                                          gAATGCTTTGCCGGGCAAATGCTACTGCATGGTGAGGAGATCCGCCAGCAGCTGCGTGGGGTGGAACTc   >  1:29096/1‑69 (MQ=255)
                                                            aCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa  <  1:3390329/68‑1 (MQ=255)
                                                            aCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa  <  1:2621723/68‑1 (MQ=255)
                                                            aCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa  <  1:2914540/68‑1 (MQ=255)
                                                            aCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa  <  1:3418362/68‑1 (MQ=255)
                                                             |                                                                  
GCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCA  >  minE/2919176‑2919303

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: