Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,153 A→C 9.7% G181G (GGT→GGG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1530AC 9.7% 156.9 / 14.3 72G181G (GGT→GGGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (32/33);  new base C (4/3);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

GGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGC  >  minE/2919083‑2919222
                                                                      |                                                                     
gggCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGa                                                                       >  1:639488/1‑71 (MQ=255)
 ggCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGAc                                                                      <  1:1041255/71‑1 (MQ=255)
 ggCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGAc                                                                      <  1:3343534/71‑1 (MQ=255)
 ggCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGAc                                                                      >  1:3722810/1‑71 (MQ=255)
    gcgGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCgg                                                                   <  1:3039677/71‑1 (MQ=255)
     cgGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCgg                                                                   <  1:3044176/70‑1 (MQ=255)
     cgGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGt                                                                  <  1:3384901/71‑1 (MQ=255)
      ggCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGt                                                                  >  1:1910813/1‑70 (MQ=255)
       gCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGt                                                                  >  1:1177411/1‑69 (MQ=255)
            ccGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgcc                                                           >  1:322476/1‑71 (MQ=255)
             cGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgcc                                                           <  1:2202825/70‑1 (MQ=255)
              gTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgccg                                                          >  1:2350324/1‑70 (MQ=255)
                aaCCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgccg                                                          >  1:941793/1‑68 (MQ=255)
                aaCCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgccg                                                          >  1:641342/1‑68 (MQ=255)
                aaCCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAg                                                       >  1:2556612/1‑71 (MQ=255)
                aaCCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGACGCAg                                                       >  1:3089947/1‑71 (MQ=255)
                 aCCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAg                                                       <  1:543678/70‑1 (MQ=255)
                 aCCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGc                                                      >  1:1150364/1‑71 (MQ=255)
                 aCCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGc                                                      <  1:3053122/71‑1 (MQ=255)
                   cAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCtt                                                    >  1:1277108/1‑71 (MQ=255)
                      cgcAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGa                                                 <  1:1326754/71‑1 (MQ=255)
                       gcAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAg                                                <  1:2420820/71‑1 (MQ=255)
                        cAGCTTCCGGCCAGCACGCTTGTGGCGCGCCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAg                                                >  1:3120473/1‑70 (MQ=255)
                            ttCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCg                                                  >  1:923315/1‑64 (MQ=255)
                             tCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATc                                            >  1:1532503/1‑69 (MQ=255)
                             tCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                           >  1:1002485/1‑70 (MQ=255)
                              ccGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGcgcc                                                           >  1:2759945/1‑53 (MQ=255)
                              ccGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                           <  1:1766788/69‑1 (MQ=255)
                              ccGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGaa                                         >  1:3638961/1‑71 (MQ=255)
                               cGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAAt                                        >  1:1093903/1‑71 (MQ=255)
                                   cAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                    <  1:258419/71‑1 (MQ=255)
                                    aGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                    <  1:3182796/70‑1 (MQ=255)
                                       aCGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtg                                >  1:3831688/1‑71 (MQ=255)
                                         gCTTGTGGCGCGACCAGACGCAAATCCGGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgtt                              <  1:750023/71‑1 (MQ=255)
                                         gCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgtt                              <  1:2695455/71‑1 (MQ=255)
                                          cTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgtt                              <  1:225363/70‑1 (MQ=255)
                                           ttGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGc                                               <  1:410862/52‑1 (MQ=255)
                                             gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                            >  1:1885405/1‑69 (MQ=255)
                                             gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                            <  1:2005970/69‑1 (MQ=255)
                                             gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATgttgtt                           <  1:2893317/70‑1 (MQ=255)
                                              tgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAc                         <  1:553622/71‑1 (MQ=255)
                                                gcgcGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGGATCGAATTGCCCATGTTGTTAcgc                       >  1:3063692/1‑71 (MQ=255)
                                                gcgcGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCg                                           <  1:2546399/51‑1 (MQ=255)
                                                 cgcgTCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttg                             <  1:1486235/64‑1 (MQ=255)
                                                  gcgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAt                                 <  1:2133395/59‑1 (MQ=255)
                                                  gcgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgcgc                     <  1:2179825/71‑1 (MQ=255)
                                                  gcgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgcgc                     <  1:2175231/71‑1 (MQ=255)
                                                   cgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATgttgtt                           >  1:783572/1‑64 (MQ=255)
                                                     aacAGGCGCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATc                  >  1:242975/3‑71 (MQ=255)
                                                     aCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGt                   >  1:1056823/1‑70 (MQ=255)
                                                     aCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTc                  <  1:1261596/71‑1 (MQ=255)
                                                      ccAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGa                                                 <  1:1829194/39‑1 (MQ=255)
                                                         gACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCt              >  1:282821/1‑71 (MQ=255)
                                                         gACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCt              >  1:3557828/1‑71 (MQ=255)
                                                          gcGCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGccc              >  1:359579/2‑69 (MQ=255)
                                                          ccGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg             >  1:652588/2‑71 (MQ=255)
                                                          aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg             >  1:281883/1‑71 (MQ=255)
                                                           cGCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGccc              <  1:1487388/69‑2 (MQ=255)
                                                           cGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGc            <  1:1517398/71‑1 (MQ=255)
                                                            gCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGcg           >  1:2437121/1‑70 (MQ=255)
                                                             cAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCAt          <  1:2048631/71‑1 (MQ=255)
                                                              agaTCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTa         <  1:1072696/69‑1 (MQ=255)
                                                              aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa           >  1:2616502/1‑69 (MQ=255)
                                                              aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa         <  1:1711077/71‑1 (MQ=255)
                                                              aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa         <  1:3293749/71‑1 (MQ=255)
                                                               aaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAg        <  1:3413603/71‑1 (MQ=255)
                                                                aTCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTgcgcgca                    <  1:3216249/58‑2 (MQ=25)
                                                                aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCGATGTTGTTACGCGCGTCACCTGCATAg        <  1:1724833/70‑1 (MQ=255)
                                                                aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGa       >  1:2904147/1‑71 (MQ=255)
                                                                aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGa       >  1:992121/1‑71 (MQ=255)
                                                                   cAGCCCGGTCAGCGCGGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGAc      >  1:2948535/1‑69 (MQ=255)
                                                                   cAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                 >  1:1661717/1‑58 (MQ=255)
                                                                     gACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGc  >  1:1582406/1‑71 (MQ=255)
                                                                      aCCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCa    <  1:383810/68‑1 (MQ=255)
                                                                      |                                                                     
GGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGC  >  minE/2919083‑2919222

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: