Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,168 A→C 10.1% A176A (GCT→GCG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1680AC10.1% 179.4 / 17.2 79A176A (GCT→GCGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base C (5/3);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

AACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACG  >  minE/2919099‑2919238
                                                                     |                                                                      
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)
              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:2175231/71‑1 (MQ=255)
                                  gcgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgcgc                                     <  1:2179825/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:3557828/1‑71 (MQ=255)
                                         gACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCt                              >  1:282821/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:3293749/71‑1 (MQ=255)
                                              aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                         <  1:1711077/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)
                                                          gtcaGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCAt             >  1:371926/4‑71 (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)
                                                             aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc          <  1:910884/71‑1 (MQ=255)
                                                             aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc          >  1:740140/1‑71 (MQ=255)
                                                             aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc          >  1:3935842/1‑71 (MQ=255)
                                                             aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGAATAGACCAGCGTCATTTc          >  1:2800822/1‑71 (MQ=255)
                                                               cgccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGGCACCTGCATAGACCAGCGTCATTTCGt        <  1:1808402/71‑1 (MQ=255)
                                                                gccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGtt       <  1:199087/71‑1 (MQ=255)
                                                                 ccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg      <  1:478987/71‑1 (MQ=255)
                                                                 ccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg      <  1:1578336/71‑1 (MQ=255)
                                                                  cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTAGGCGCGTCACCt                              <  1:1501341/46‑1 (MQ=255)
                                                                   gcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGttaac    >  1:3249707/1‑68 (MQ=255)
                                                                     aGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACg  >  1:910126/1‑71 (MQ=255)
                                                                     |                                                                      
AACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACG  >  minE/2919099‑2919238

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: