Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,174 G→C 10.1% L174L (CTC→CTG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1740GC10.1% 250.2 / 15.2 89L174L (CTC→CTGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (39/41);  new base C (6/3);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

CGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTG  >  minE/2919105‑2919243
                                                                     |                                                                     
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: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)
                                                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:740140/1‑71 (MQ=255)
                                                       aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc               >  1:3935842/1‑71 (MQ=255)
                                                       aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc               <  1:910884/71‑1 (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:1578336/71‑1 (MQ=255)
                                                           ccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg           <  1:478987/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)
                                                                gCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGc      >  1:2102341/1‑71 (MQ=255)
                                                                 cTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGc                                 <  1:3952579/43‑1 (MQ=255)
                                                                 cTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCt     <  1:1992691/71‑1 (MQ=255)
                                                                 cTTCGAGCAGCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGt             <  1:365812/63‑1 (MQ=255)
                                                                 cTTCGAGCAGCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGt             <  1:361106/63‑1 (MQ=255)
                                                                  ttCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCtt    >  1:1581253/1‑71 (MQ=255)
                                                                  ttCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCtt    >  1:1117797/1‑71 (MQ=255)
                                                                   tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                                      >  1:195632/1‑36 (MQ=255)
                                                                   tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                >  1:2467971/1‑42 (MQ=255)
                                                                   tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg           >  1:2810286/1‑63 (MQ=255)
                                                                   tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg           >  1:1166324/1‑63 (MQ=255)
                                                                    cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                                      <  1:1986271/35‑1 (MQ=255)
                                                                    cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt   <  1:839445/70‑1 (MQ=255)
                                                                    cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt   <  1:299123/70‑1 (MQ=255)
                                                                    cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCt     >  1:3850122/1‑68 (MQ=255)
                                                                    cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCt     >  1:3835298/1‑68 (MQ=255)
                                                                     gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg  <  1:3146009/70‑1 (MQ=255)
                                                                     gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg  <  1:1436368/70‑1 (MQ=255)
                                                                     |                                                                     
CGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTG  >  minE/2919105‑2919243

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: