Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,183 A→G 10.8% N171N (AAT→AAC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1830AG10.8% 202.4 / 14.6 74N171N (AAT→AACargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (27/39);  new base G (5/3);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAA  >  minE/2919115‑2919251
                                                                    |                                                                    
ggCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAAtt                                                                    <  1:2481428/71‑1 (MQ=255)
   cAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                                                 <  1:616015/71‑1 (MQ=255)
        cGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgt                                                            <  1:2825293/71‑1 (MQ=255)
        cGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgt                                                            <  1:1724167/71‑1 (MQ=255)
             gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                                                         <  1:672418/69‑1 (MQ=255)
             gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                                                         <  1:1649156/69‑1 (MQ=255)
             gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATgttgtt                                                        <  1:3894669/70‑1 (MQ=255)
              tgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAc                                                      <  1:3365917/71‑1 (MQ=255)
                gcgcGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgc                                                    >  1:4046496/1‑71 (MQ=255)
                    gACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGt                                                <  1:1252158/71‑1 (MQ=255)
                      ccAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                                              <  1:865859/71‑1 (MQ=255)
                      acAGGCGCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATcg                                              <  1:1698637/70‑2 (MQ=255)
                        aGACGCAAATCCAGACCGGTAAGCGCCGCATCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCAcc                                            >  1:636375/1‑71 (MQ=255)
                          aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg                                          <  1:3903685/71‑1 (MQ=255)
                          aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg                                          <  1:2070343/71‑1 (MQ=255)
                            gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTg                                                                   >  1:2229634/1‑44 (MQ=255)
                            gCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTc                                               <  1:2806903/64‑1 (MQ=255)
                            gCAAATCCAGACCGGTAAGCGCCGCACCTTCGAGCATCGAATTGCCCATGTTGTTAc                                                      >  1:3790474/1‑57 (MQ=255)
                             cAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGTGCGTCACCTGCa                                        >  1:3907997/1‑70 (MQ=255)
                             cAAATCCAGACCGGTAAGAGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                        >  1:1794130/1‑70 (MQ=255)
                              aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGaa                                                                      >  1:3457089/1‑39 (MQ=255)
                              aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                                      <  1:3379729/71‑1 (MQ=255)
                               aaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcg                                                     <  1:3282971/55‑1 (MQ=255)
                                aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                                                 <  1:1036567/42‑1 (MQ=255)
                                 tCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                        <  1:1816826/66‑1 (MQ=255)
                                 tCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                        <  1:2544179/66‑1 (MQ=255)
                                 tCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                        <  1:1135382/66‑1 (MQ=255)
                                  ccAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAt                                                              <  1:968806/43‑1 (MQ=255)
                                   cAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGAcc                                  <  1:2526445/70‑1 (MQ=255)
                                        cGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTc                            <  1:714916/71‑1 (MQ=255)
                                          gTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAt                          <  1:3430471/71‑1 (MQ=255)
                                          gTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAt                          >  1:2721947/1‑71 (MQ=255)
                                           tcaGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCAtc                         >  1:3373137/3‑70 (MQ=255)
                                             aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgcgc                                                  >  1:2090193/1‑44 (MQ=255)
                                             aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                       >  1:2368877/1‑71 (MQ=255)
                                             aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                       >  1:73158/1‑71 (MQ=255)
                                              gcgcCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCa                                 <  1:1950686/60‑1 (MQ=255)
                                              gcgcCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCa                                 <  1:1818742/60‑1 (MQ=255)
                                               cgccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt                             >  1:3649841/1‑63 (MQ=255)
                                                 ccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGtt                    >  1:929025/1‑70 (MQ=255)
                                                  cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCAcc                                            <  1:3324485/45‑1 (MQ=255)
                                                  cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGa                  <  1:2055691/71‑1 (MQ=255)
                                                  cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGa                  <  1:1296743/71‑1 (MQ=255)
                                                    cAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                        >  1:632986/1‑47 (MQ=255)
                                                    cAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                        >  1:1252383/1‑47 (MQ=255)
                                                      gCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                        <  1:3639030/45‑1 (MQ=255)
                                                      gCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGc              >  1:866859/1‑71 (MQ=255)
                                                        ttCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                                      >  1:1838680/1‑45 (MQ=255)
                                                         tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                                              >  1:2707584/1‑36 (MQ=255)
                                                         tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGc                                         >  1:3163876/1‑41 (MQ=255)
                                                         tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                        >  1:3159334/1‑42 (MQ=255)
                                                         tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                        >  1:444015/1‑42 (MQ=255)
                                                         tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                        >  1:2039151/1‑42 (MQ=255)
                                                         tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt           >  1:2638871/1‑71 (MQ=255)
                                                          ccaGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTg          >  1:2298655/3‑71 (MQ=255)
                                                          ccaGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTg          >  1:573383/3‑71 (MQ=255)
                                                          cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt           <  1:3383346/70‑1 (MQ=255)
                                                           gAGCATCGAATTGCCCATGTTGTTAGGCGTGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc         <  1:1053910/71‑1 (MQ=255)
                                                           gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc         >  1:3029943/1‑71 (MQ=255)
                                                            aGCATCGAATTGCCCATGTTGTTAAGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc        <  1:2289205/71‑1 (MQ=255)
                                                             gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCg       <  1:902106/71‑1 (MQ=255)
                                                              cATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGct                                        <  1:959661/37‑2 (MQ=16)
                                                              cATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc         <  1:3393915/68‑1 (MQ=255)
                                                                tCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGttaa                  >  1:3008743/1‑55 (MQ=255)
                                                                tCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                   <  1:1649990/56‑1 (MQ=255)
                                                                tCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                   <  1:3627368/56‑1 (MQ=255)
                                                                  gAATTGCCCATGTTGTTACGCTCGTCACCTGCATAGACCAGGGTCATTTCGTTGAACGCTTTGCCGGGCaa  <  1:4103696/71‑1 (MQ=255)
                                                                  gAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                                      >  1:1329128/1‑35 (MQ=255)
                                                                  gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGaa                 >  1:547197/1‑56 (MQ=255)
                                                                  gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa   <  1:598495/70‑1 (MQ=255)
                                                                  gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa   >  1:481889/1‑70 (MQ=255)
                                                                  gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa   <  1:2921634/70‑1 (MQ=255)
                                                                  gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa   <  1:132593/70‑1 (MQ=255)
                                                                  gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGAGTCATTTCGTTGAACGCTTTGCCGGGCa   >  1:2397820/1‑70 (MQ=255)
                                                                  gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCa   >  1:93207/1‑70 (MQ=255)
                                                                    aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt                             >  1:393616/1‑42 (MQ=255)
                                                                    aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt                             >  1:1146533/1‑42 (MQ=255)
                                                                    aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGc    >  1:1354108/1‑67 (MQ=255)
                                                                    |                                                                    
GGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAA  >  minE/2919115‑2919251

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: