Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,918,982 C→G 15.3% G238G (GGG→GGC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,918,9820CG15.3% 205.2 / 25.8 85G238G (GGG→GGCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (43/29);  new base G (9/4);  total (52/33)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-01

CATCATCTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGT  >  minE/2918913‑2919050
                                                                     |                                                                    
catcatCTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTccc                                                                     >  1:1025794/1‑71 (MQ=255)
 atcatcTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTcccc                                                                    <  1:432329/71‑1 (MQ=255)
   catcTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCa                                                                  >  1:2989203/1‑71 (MQ=255)
    atcTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCa                                                                  <  1:82158/70‑1 (MQ=255)
    atcTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCAt                                                                 <  1:1746360/71‑1 (MQ=255)
    atcTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCAt                                                                 <  1:170022/71‑1 (MQ=255)
    atcTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCAt                                                                 >  1:2185739/1‑71 (MQ=255)
       ttGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGa                                                              >  1:193804/1‑71 (MQ=255)
        tGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGa                                                              >  1:2493869/1‑70 (MQ=255)
         cgcgTTCACCTGATACCCGCGCAGCAGCGCAATCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGaca                                                            <  1:1200300/68‑1 (MQ=255)
             ttCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCca                                                        >  1:2331851/1‑71 (MQ=255)
             ttCACCTGATATTCACGCAGCAATGCAATCAGTTCCGCCCATTTCTCTTTTGCTTACCCCATCGACAcc                                                          >  1:3538357/1‑69 (MQ=255)
              tCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCca                                                        >  1:481678/1‑70 (MQ=255)
              tCACCTGATACCCGCGCAGCAGCGCAATCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCca                                                        >  1:2568887/1‑70 (MQ=255)
               cACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGac                                                             >  1:1427615/1‑64 (MQ=255)
               cACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCcaca                                                      >  1:3413846/1‑71 (MQ=255)
               cACCTGATACCCGCGCAGCAGCGCAATCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCcaca                                                      >  1:1024526/1‑71 (MQ=255)
                aCCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCcaca                                                      <  1:3949957/70‑1 (MQ=255)
                 ccTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCcacaca                                                    <  1:969493/71‑1 (MQ=255)
                  cTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACAt                                                   >  1:1749662/1‑71 (MQ=255)
                  cTGATACCCGCGCAGCAGCGCAATCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCcacacgt                                                   >  1:2889099/1‑69 (MQ=255)
                    gATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCg                                                 >  1:1249276/1‑71 (MQ=255)
                     atatTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATc                                                  <  1:355764/69‑1 (MQ=255)
                     atatTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATc                                                  <  1:1851072/69‑1 (MQ=255)
                      tatTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGaca                                                            >  1:3509406/1‑58 (MQ=255)
                      tatTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCgg                                                <  1:50180/70‑1 (MQ=255)
                        ttCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGta                                              <  1:1807509/70‑1 (MQ=255)
                        ttCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGta                                              <  1:1984513/70‑1 (MQ=255)
                         tCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGa                                                              >  1:3447658/1‑53 (MQ=255)
                         ccgcGCAGCAGCGCAATCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCCACACGTCGGtata                                            <  1:1419790/68‑1 (MQ=255)
                          cACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATc                                                                >  1:595954/1‑50 (MQ=255)
                          cACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGta                                              >  1:2973078/1‑68 (MQ=255)
                          cACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGta                                              >  1:3654838/1‑68 (MQ=255)
                           aCGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTcc                                                                      >  1:2217115/1‑43 (MQ=255)
                            cGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGCCACCCACACATCGGTATAGa                                          >  1:2953382/1‑70 (MQ=255)
                             gcagcaGCGCAATCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGAt                                         >  1:651492/1‑70 (MQ=255)
                                gcaatgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAg                                     <  1:819556/71‑1 (MQ=255)
                                 caatgcaatCCGTTCCGCCAATTTCTCTTTTGCTTCCCCCATCGACACCcacaca                                                    >  1:1198134/1‑55 (MQ=255)
                                  aatgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGt                                    <  1:3871636/70‑1 (MQ=255)
                                  aatgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGt                                    >  1:302257/1‑70 (MQ=255)
                                    tgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGt                                    <  1:1206775/68‑1 (MQ=255)
                                    tgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGt                                    <  1:3139908/68‑1 (MQ=255)
                                    tgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTc                                   >  1:1119927/1‑69 (MQ=255)
                                     gcaatCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCcaca                                                      >  1:3708127/1‑49 (MQ=19)
                                     gcaatCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCcaca                                                      >  1:2594350/1‑49 (MQ=19)
                                     gcaatCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCcaca                                                      >  1:314329/1‑49 (MQ=19)
                                      caatCCGTTCCGCCCATTTCTCTTTTGCTTccccc                                                                   <  1:1752156/35‑1 (MQ=255)
                                      caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCgg                                                >  1:1299522/1‑54 (MQ=255)
                                      caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCgg                                                >  1:3543184/1‑54 (MQ=255)
                                      caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCa                               >  1:893447/1‑71 (MQ=255)
                                        atCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCAc                              >  1:2022788/1‑70 (MQ=255)
                                         tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCg                                                 >  1:71625/1‑50 (MQ=255)
                                         tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCAcc                             >  1:1399315/1‑70 (MQ=255)
                                         tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCAcc                             >  1:3408915/1‑70 (MQ=255)
                                         tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCt                            <  1:1084522/71‑1 (MQ=255)
                                         tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCt                            <  1:3320291/71‑1 (MQ=255)
                                         tCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGccc                            >  1:2377632/1‑70 (MQ=255)
                                         tCCGCTCTGCCCACTTCTCTCTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGccc                            >  1:2449673/1‑70 (MQ=255)
                                           cGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCca                                                        >  1:954254/1‑41 (MQ=255)
                                           cGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTc                          <  1:657793/71‑1 (MQ=255)
                                           cGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTc                          <  1:1222338/71‑1 (MQ=255)
                                               ccGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACt                      >  1:3360828/1‑71 (MQ=255)
                                               ccGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACt                      >  1:3549542/1‑71 (MQ=255)
                                                tgCCCACTTCTCTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGCCCTTAAc                       <  1:1736907/68‑1 (MQ=255)
                                                tgCCCACTTCTCTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGCCCTTAAc                       <  1:2711405/68‑1 (MQ=255)
                                                cGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAg                                           <  1:1730652/49‑1 (MQ=255)
                                                cGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAAc                       <  1:791799/69‑1 (MQ=255)
                                                cGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACt                      <  1:82081/70‑1 (MQ=255)
                                                 gCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTcc                    >  1:2737717/1‑71 (MQ=255)
                                                 gCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTcc                    >  1:2096706/1‑71 (MQ=255)
                                                 gCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTcc                    >  1:303929/1‑71 (MQ=255)
                                                    cATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCtt                 <  1:2408289/71‑1 (MQ=255)
                                                    cATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCtt                 >  1:1258520/1‑71 (MQ=255)
                                                     aTTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTc                >  1:865566/1‑71 (MQ=255)
                                                      tttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTc                          <  1:3778258/60‑1 (MQ=255)
                                                       ttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGt                                    <  1:2784557/49‑1 (MQ=255)
                                                        tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAg                                 >  1:3644837/1‑51 (MQ=255)
                                                        tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCaa                        >  1:271003/1‑60 (MQ=255)
                                                        tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc              <  1:3294383/70‑1 (MQ=255)
                                                            ttttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGc                                >  1:235068/1‑48 (MQ=255)
                                                            ttttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGc                                >  1:1968574/1‑48 (MQ=255)
                                                            ttttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTc                >  1:3258529/1‑64 (MQ=255)
                                                             tttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcg               >  1:155172/1‑64 (MQ=255)
                                                             tttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATc        >  1:2388753/1‑71 (MQ=255)
                                                             tttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATc        >  1:2393502/1‑71 (MQ=255)
                                                              ttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATc        <  1:731308/70‑1 (MQ=255)
                                                                gCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCtt      <  1:1767907/70‑1 (MQ=255)
                                                                 cTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtc     <  1:3931114/70‑1 (MQ=255)
                                                                  ttCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTc                          <  1:791210/48‑1 (MQ=255)
                                                                    cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc              <  1:1820010/58‑1 (MQ=255)
                                                                    cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt  >  1:1219262/1‑70 (MQ=255)
                                                                    cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt  >  1:2863978/1‑70 (MQ=255)
                                                                    cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt  >  1:95989/1‑70 (MQ=255)
                                                                    cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt  >  1:323326/1‑70 (MQ=255)
                                                                    cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcag   <  1:2110952/69‑1 (MQ=255)
                                                                     ccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCt                  >  1:326320/1‑53 (MQ=255)
                                                                     |                                                                    
CATCATCTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGT  >  minE/2918913‑2919050

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: