Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,291 G→C 8.7% P135P (CCC→CCG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2910GC 8.7% 226.9 / 9.7 80P135P (CCC→CCGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base C (2/5);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAT  >  minE/2919224‑2919359
                                                                   |                                                                    
tcctttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa                                                                   <  1:1672665/68‑1 (MQ=255)
tCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa                                                                   <  1:1495109/71‑1 (MQ=255)
 cATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt                                                                  >  1:3709877/1‑71 (MQ=255)
  atctCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTg                                                                 <  1:2544862/68‑1 (MQ=255)
  aTTTCGTTGAACGCTTTGCCTGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                                                                 <  1:1517941/71‑1 (MQ=255)
   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGTGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                <  1:2800/71‑1 (MQ=255)
   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                <  1:291035/71‑1 (MQ=255)
   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                <  1:1516533/71‑1 (MQ=255)
   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGAATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                <  1:2940543/71‑1 (MQ=255)
      cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                                             >  1:464517/1‑71 (MQ=255)
      cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                                             >  1:3500525/1‑71 (MQ=255)
       gTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                                              <  1:3462692/69‑1 (MQ=255)
        ttGAACGCTTTGCCGGGCAAATGGTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                >  1:16463/1‑66 (MQ=255)
        ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                                           <  1:1014133/71‑1 (MQ=255)
        ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                                           <  1:1671315/71‑1 (MQ=255)
         tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGTATGGAACTCa                                                          >  1:242989/1‑71 (MQ=255)
         tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                                           <  1:856089/70‑1 (MQ=255)
         tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                                                          >  1:632586/1‑71 (MQ=255)
          gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATTGAACTCa                                                          <  1:2928175/70‑1 (MQ=255)
           aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                <  1:3392923/63‑1 (MQ=255)
           aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                         <  1:935183/70‑1 (MQ=255)
           aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                                                         <  1:1989517/70‑3 (MQ=255)
           aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                                                         <  1:3478967/70‑3 (MQ=255)
              gCTTTGCCGGGCAATTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                     <  1:2727425/71‑1 (MQ=255)
              gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                     <  1:2405156/71‑1 (MQ=255)
              gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                     <  1:56103/71‑1 (MQ=255)
                tttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                     >  1:25113/1‑69 (MQ=255)
                cttGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTCGAACTCGTTGGt                                                     >  1:260038/2‑69 (MQ=255)
                 ttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                                                  >  1:1152492/1‑71 (MQ=255)
                  tGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg                                                 >  1:983456/1‑71 (MQ=255)
                   gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                                >  1:1081914/1‑71 (MQ=255)
                   gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                                >  1:356765/1‑71 (MQ=255)
                     cGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGcc                                               <  1:1194635/70‑1 (MQ=255)
                      gggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAt                                             <  1:2248350/71‑1 (MQ=255)
                       ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt                                            <  1:3517193/71‑1 (MQ=255)
                       ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt                                            <  1:3875150/71‑1 (MQ=255)
                        gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc                                           >  1:2506764/1‑71 (MQ=255)
                         cAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc                                           >  1:3647491/1‑70 (MQ=255)
                         cAAATGCTCCTGCATGGTGAGAAGATCCGCAAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc                                          <  1:2330912/71‑1 (MQ=255)
                          aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc                                          >  1:1394035/1‑70 (MQ=255)
                          aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa                                         >  1:3419366/1‑71 (MQ=255)
                           aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                                                               >  1:301951/1‑48 (MQ=255)
                           aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                                                               >  1:2164908/1‑48 (MQ=255)
                           aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAtt                                                        >  1:2429364/1‑55 (MQ=255)
                           aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                                        <  1:1589273/71‑1 (MQ=255)
                           aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                                        <  1:2438395/71‑1 (MQ=255)
                            aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATa                                       <  1:738670/71‑1 (MQ=255)
                             tGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                                              >  1:278211/1‑47 (MQ=16)
                              gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                                                                 >  1:3463126/1‑43 (MQ=255)
                              gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAcc                                     >  1:2131374/1‑71 (MQ=255)
                              gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAcc                                     >  1:1154139/1‑71 (MQ=255)
                                tCCTGCATGGTGAGAAGATCCGCCAGCAGGTGCGTGGGATGGAACTCATTGGTCAgg                                                 >  1:1842506/1‑57 (MQ=255)
                                tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                                                  <  1:3008816/56‑1 (MQ=255)
                                tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                                                  <  1:4085/56‑1 (MQ=255)
                                 ccTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATTGAACTCATTGGTCAg                                                  <  1:3099690/55‑1 (MQ=255)
                                  cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCa                                 >  1:2599004/1‑71 (MQ=255)
                                  cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGGCAGGCCATTCCATACCgg                                   <  1:2504979/69‑1 (MQ=255)
                                  cTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCACACCGGc                                  <  1:2359793/70‑1 (MQ=255)
                                    gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg                               >  1:1567680/1‑71 (MQ=255)
                                    gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg                               >  1:254901/1‑71 (MQ=255)
                                     cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                              <  1:985564/71‑1 (MQ=255)
                                     cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                              <  1:1445464/71‑1 (MQ=255)
                                      aTGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                              >  1:1852288/1‑70 (MQ=255)
                                          tGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCg                         <  1:236448/71‑1 (MQ=255)
                                            agaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                         <  1:2142807/37‑1 (MQ=255)
                                            agaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTa                       <  1:3870150/71‑1 (MQ=255)
                                             gaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAt                      >  1:2488977/1‑71 (MQ=255)
                                             gaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAt                      <  1:1442393/71‑1 (MQ=255)
                                               agaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTc                    >  1:2220404/1‑71 (MQ=255)
                                                 aTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTc                                                    <  1:627250/37‑1 (MQ=255)
                                                  tCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGc                 >  1:2926612/1‑71 (MQ=255)
                                                    cGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCa               <  1:2358888/71‑1 (MQ=255)
                                                       cagcagCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGcc                <  1:1370370/67‑1 (MQ=255)
                                                          cagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGttt         >  1:2596451/1‑71 (MQ=255)
                                                          cagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGttt         >  1:2651682/1‑71 (MQ=255)
                                                           agcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTc        <  1:568277/71‑1 (MQ=255)
                                                             cTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGttt         <  1:2716573/68‑1 (MQ=255)
                                                                cGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACaa   >  1:3275705/1‑71 (MQ=255)
                                                                agTGGGATGGAACTCATTGGTCAGGCCATTCCATACAGGCACGCTAGCGt                        >  1:3782737/2‑50 (MQ=255)
                                                                 gTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAt                      >  1:3475825/1‑51 (MQ=255)
                                                                   gggATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAt  >  1:1176296/1‑69 (MQ=255)
                                                                   |                                                                    
TCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAT  >  minE/2919224‑2919359

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