Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,303 A→G 13.1% N131N (AAT→AAC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,3030AG13.1% 204.1 / 19.4 84N131N (AAT→AACargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (39/34);  new base G (7/4);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGC  >  minE/2919233‑2919373
                                                                      |                                                                      
tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                                                                        >  1:94922/1‑71 (MQ=255)
 gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                                       >  1:3541285/1‑71 (MQ=255)
 gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                                       >  1:641879/1‑71 (MQ=255)
    cGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTgg                                                                    >  1:1824122/1‑71 (MQ=255)
     gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                                   <  1:3591007/71‑1 (MQ=255)
          gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCTAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg                                                               <  1:3382133/70‑1 (MQ=255)
          gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg                                                               <  1:5174129/70‑1 (MQ=255)
           ccGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                                              <  1:3721670/70‑1 (MQ=255)
           ccGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                                              <  1:4659488/70‑1 (MQ=255)
               gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc                                                         <  1:1168233/71‑1 (MQ=255)
               gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc                                                         <  1:71432/71‑1 (MQ=255)
                  aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                                       >  1:724246/1‑54 (MQ=255)
                  aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                                       >  1:682287/1‑54 (MQ=255)
                  aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                                       >  1:5669422/1‑54 (MQ=255)
                  aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                                                      <  1:2454188/71‑1 (MQ=255)
                   gtGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGcccg                                                            <  1:1420698/63‑2 (MQ=255)
                    tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                                                                        >  1:3068557/1‑51 (MQ=255)
                    tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                                                                        >  1:4413032/1‑51 (MQ=255)
                    tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                                                                        >  1:2230691/1‑51 (MQ=255)
                    tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAc                                                    >  1:3789684/1‑71 (MQ=255)
                    tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAc                                                    <  1:5521335/71‑1 (MQ=255)
                    tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAc                                                    >  1:4737010/1‑71 (MQ=255)
                    tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAc                                                    >  1:3730550/1‑71 (MQ=255)
                    tGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCacac                                                    >  1:2477149/1‑71 (MQ=255)
                     gCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCa                                                       >  1:3222573/1‑67 (MQ=255)
                      cTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAtt                                                                      >  1:5339647/1‑51 (MQ=255)
                      cTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAtt                                                                      >  1:2391817/1‑51 (MQ=255)
                      cTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAtt                                                                      >  1:1986079/1‑51 (MQ=255)
                      cTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCa                                                                 <  1:3334773/56‑1 (MQ=255)
                       tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCgg                                                 >  1:1028938/1‑71 (MQ=255)
                        ccTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCACACCgg                                                 >  1:4908013/1‑70 (MQ=255)
                         cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                                                      <  1:3892490/64‑1 (MQ=255)
                         cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCa                                               >  1:5638374/1‑71 (MQ=255)
                         cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCa                                               >  1:4673446/1‑71 (MQ=255)
                         cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCa                                               >  1:4666714/1‑71 (MQ=255)
                         cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCa                                               >  1:1223372/1‑71 (MQ=255)
                         cTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCACGTTCCACACCGGCa                                               >  1:5084020/1‑71 (MQ=255)
                          tGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCac                                                      <  1:5158794/63‑2 (MQ=255)
                            cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                                            >  1:1811018/1‑71 (MQ=255)
                            cATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCaca                                                     >  1:4321845/1‑60 (MQ=255)
                            cATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCACACCGGCACGc                                            >  1:3624200/1‑71 (MQ=255)
                            cATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCACACCGGCACGc                                            >  1:666396/1‑71 (MQ=255)
                             aTGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                                            >  1:4712972/1‑70 (MQ=255)
                             aTGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCt                                           <  1:1884845/71‑1 (MQ=255)
                               ggTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                                   <  1:4166004/45‑1 (MQ=255)
                               ggTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAg                                         >  1:4734366/1‑71 (MQ=255)
                                gTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                                              <  1:2185865/49‑1 (MQ=255)
                                gTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCt                                           <  1:5235145/68‑1 (MQ=255)
                                 tGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCg                                                  <  1:1008725/60‑1 (MQ=255)
                                   agaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTa                                     >  1:3144281/1‑71 (MQ=255)
                                       gaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTcc                                 <  1:492190/71‑1 (MQ=255)
                                       gaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTcc                                 <  1:1094994/71‑1 (MQ=255)
                                        aTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa                                                       >  1:4084968/1‑48 (MQ=255)
                                        aTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa                                                       >  1:1684915/1‑48 (MQ=255)
                                         tCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                                              >  1:1759355/1‑40 (MQ=255)
                                         tCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                                              >  1:248943/1‑40 (MQ=255)
                                          ccGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCa                                               >  1:1314175/1‑54 (MQ=255)
                                          ccGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGc                               >  1:1182776/1‑70 (MQ=255)
                                            gCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAc                                                    >  1:1855641/1‑47 (MQ=255)
                                            gCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAg                            >  1:256341/1‑71 (MQ=255)
                                             ccAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAgt                           <  1:4649019/71‑1 (MQ=255)
                                                gcagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAgtg                          >  1:4436788/1‑69 (MQ=255)
                                                 cagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                                                      <  1:2364871/40‑1 (MQ=255)
                                                   gcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCg                     <  1:540105/71‑1 (MQ=255)
                                                    cTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc                                                        >  1:1935097/1‑35 (MQ=255)
                                                      gCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCg                                                  >  1:4106739/1‑39 (MQ=255)
                                                      gCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGAc                   >  1:4554183/1‑70 (MQ=255)
                                                       cGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACaa                 <  1:4943226/71‑1 (MQ=255)
                                                       cGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACa                  <  1:1013121/70‑1 (MQ=255)
                                                       cGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACa                  <  1:5410958/70‑1 (MQ=255)
                                                       cGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACa                  <  1:2351598/70‑1 (MQ=255)
                                                        gTGGGATGGAACTCATTGGTCAGGCCATTCCATACGGGCACGCTAGCGTATTCCGCCAGTGTTTCGACa                  <  1:980047/69‑1 (MQ=255)
                                                         tGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAtc               <  1:4858629/71‑1 (MQ=255)
                                                             aTGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCt           <  1:773674/71‑1 (MQ=255)
                                                              tGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGc                                        >  1:466979/1‑41 (MQ=255)
                                                              tGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCt           >  1:2501416/1‑70 (MQ=255)
                                                              tGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCt           >  1:4717303/1‑70 (MQ=255)
                                                                gAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGAc        >  1:1907166/1‑71 (MQ=255)
                                                                  aCTCATTGGTCAGGCCATTCCATACCTGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCa      <  1:58522/71‑1 (MQ=255)
                                                                   cTCGTTGGTCAGCCCGTTCCACACCGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGccgt     <  1:2124962/71‑3 (MQ=255)
                                                                   cTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCAt     <  1:3328650/71‑1 (MQ=255)
                                                                   cTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCAt     <  1:2992503/71‑1 (MQ=255)
                                                                   cTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCAt     <  1:1589782/71‑1 (MQ=255)
                                                                    tCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCg                                <  1:4003715/43‑1 (MQ=255)
                                                                    tCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCg                                <  1:2920689/43‑1 (MQ=255)
                                                                    tCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCg                                <  1:1892150/43‑1 (MQ=255)
                                                                    tCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGAc                   >  1:212285/1‑56 (MQ=255)
                                                                     cATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAg   >  1:1217205/1‑71 (MQ=255)
                                                                      aTTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAtctc             >  1:4668015/1‑60 (MQ=255)
                                                                      aTTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGc  >  1:5674693/1‑71 (MQ=255)
                                                                      aTTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGc  >  1:5324878/1‑71 (MQ=255)
                                                                      |                                                                      
TGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGC  >  minE/2919233‑2919373

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: