Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,159 A→C 10.0% L179L (CTT→CTG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1590AC10.0% 158.8 / 14.4 70L179L (CTT→CTGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base C (4/3);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

CCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTT  >  minE/2919095‑2919229
                                                                |                                                                      
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:992121/1‑71 (MQ=255)
                                                    aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGa              >  1:2904147/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)
                                                          acaggTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCg        >  1:18140/4‑71 (MQ=255)
                                                          aCCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCa           <  1:383810/68‑1 (MQ=255)
                                                           ccGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt       >  1:3004126/1‑71 (MQ=255)
                                                           ccGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt       <  1:205712/71‑1 (MQ=255)
                                                            cGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAc                                <  1:1661539/45‑1 (MQ=255)
                                                            cGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt       >  1:3549354/1‑70 (MQ=255)
                                                                aaGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGAc             >  1:2726544/1‑60 (MQ=255)
                                                                aaGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTc      >  1:3090568/1‑67 (MQ=255)
                                                                aaGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAttt  >  1:721565/1‑71 (MQ=255)
                                                                aaGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAtt   >  1:3856891/1‑70 (MQ=255)
                                                                |                                                                      
CCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTT  >  minE/2919095‑2919229

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: