Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,918,751 A→G 17.1% D315D (GAT→GAC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,918,7510AG17.1% 158.9 / 24.9 70D315D (GAT→GACargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base G (7/5);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATT  >  minE/2918686‑2918821
                                                                 |                                                                      
aGCAATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaa                                                                    >  1:2296540/1‑70 (MQ=255)
aGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaa                                                                    >  1:368211/1‑70 (MQ=255)
aGCAATTATTTACTGAGCGTCGAGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaa                                                                    >  1:1151857/1‑70 (MQ=255)
 gCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaaaa                                                                  <  1:1594911/71‑1 (MQ=255)
    attattTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACa                                                               <  1:2201761/71‑1 (MQ=255)
    attattTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACa                                                               <  1:1452984/71‑1 (MQ=255)
       attTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAAt                                                             >  1:1688317/1‑70 (MQ=255)
          tACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCt                                                          <  1:3499375/70‑1 (MQ=255)
          tACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                                                         >  1:1084860/1‑71 (MQ=255)
          tACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                                                         >  1:3032873/1‑71 (MQ=255)
          tACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                                                         >  1:795772/1‑71 (MQ=255)
           aCTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                                                         <  1:609998/70‑1 (MQ=255)
             tGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaa                                                                    <  1:2005096/57‑1 (MQ=255)
             tGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcg                                                      >  1:1903999/1‑71 (MQ=255)
             caaGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgcc                                                      >  1:527391/3‑70 (MQ=255)
              aaGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgccg                                                     >  1:1215440/2‑69 (MQ=255)
               aGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTcg                                                                     <  1:1579608/54‑2 (MQ=255)
               aGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcgg                                                     >  1:1836579/1‑70 (MQ=255)
               aGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcgg                                                     <  1:3623585/70‑1 (MQ=255)
                   tCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAt                                                <  1:2444571/71‑1 (MQ=255)
                     gcgACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTc                                              >  1:1460201/1‑71 (MQ=255)
                      cgACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTc                                              >  1:1284918/1‑70 (MQ=255)
                        aCCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCg                                             >  1:3350507/1‑69 (MQ=255)
                         tcATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGAt                                                             >  1:2305279/2‑52 (MQ=255)
                          cATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgccgccg                                                  >  1:511659/1‑60 (MQ=255)
                           aTCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCaaaca                                         <  1:2276540/70‑3 (MQ=255)
                            tCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGa                                            <  1:2097934/66‑1 (MQ=255)
                               ccGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtca                                    <  1:538164/71‑1 (MQ=255)
                                cGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGa                                         >  1:1843292/1‑65 (MQ=255)
                                cGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTcg                                                                     <  1:3110759/37‑2 (MQ=16)
                                 gCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcat                                   >  1:2259504/1‑70 (MQ=255)
                                 gCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcat                                   >  1:2800896/1‑70 (MQ=255)
                                    ttGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAg                                >  1:2546262/1‑70 (MQ=255)
                                     tGATAGTATGCATACGGTTTTCCGCCTGATCAAAAAc                                                                <  1:1141652/37‑1 (MQ=255)
                                     tGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCACATCAGt                               >  1:2475180/1‑70 (MQ=255)
                                     tGATAGTATGCATACGGTTTTCCACCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGt                               <  1:1433903/70‑1 (MQ=255)
                                          gTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTc                         <  1:2289895/71‑1 (MQ=255)
                                           tATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTcc                        <  1:3665075/71‑1 (MQ=255)
                                           tATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTc                         >  1:2370965/1‑70 (MQ=255)
                                             tGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcg                                                      <  1:2991960/39‑1 (MQ=255)
                                             tGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCAt                      <  1:2537598/71‑1 (MQ=255)
                                              gCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGaa                                           >  1:3323714/1‑49 (MQ=255)
                                              gCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATa                     >  1:2586964/1‑71 (MQ=255)
                                               cATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCAtcc                    >  1:3815328/1‑69 (MQ=255)
                                               cATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGa                                            <  1:231261/47‑1 (MQ=255)
                                                 tACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAtt                                               <  1:519165/42‑1 (MQ=255)
                                                 tACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATAc                    <  1:3482601/69‑1 (MQ=255)
                                                  ccGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCAcct                          >  1:456762/2‑60 (MQ=255)
                                                  ccGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCATCCCGc                 >  1:1610999/2‑71 (MQ=255)
                                                  aCGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGa                                         >  1:627506/1‑47 (MQ=255)
                                                  aCGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGc                 >  1:3433599/1‑71 (MQ=255)
                                                   cGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCATCCCGcc                <  1:2703955/71‑1 (MQ=255)
                                                   cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGAcc                                       >  1:1417919/1‑48 (MQ=255)
                                                   cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcatca                                 >  1:880587/1‑54 (MQ=255)
                                                   cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATa                     >  1:483246/1‑66 (MQ=255)
                                                   cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGcc                <  1:985325/71‑1 (MQ=255)
                                                    ggTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCa               >  1:1931823/1‑71 (MQ=255)
                                                     gTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCAt              <  1:2494147/71‑1 (MQ=255)
                                                      ttttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATg             <  1:3403002/71‑1 (MQ=255)
                                                        ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGaa                                           >  1:1991256/1‑39 (MQ=255)
                                                        ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATg             >  1:883096/1‑69 (MQ=255)
                                                        ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt            <  1:3720328/70‑1 (MQ=255)
                                                        ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt            <  1:3396146/70‑1 (MQ=255)
                                                         tCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTa           <  1:3683894/70‑1 (MQ=255)
                                                           cGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTc                         <  1:2049051/54‑1 (MQ=255)
                                                           cGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTc                         <  1:860264/54‑1 (MQ=255)
                                                             ccTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGc        <  1:3160204/69‑1 (MQ=255)
                                                             ccTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa      <  1:3791017/71‑1 (MQ=255)
                                                             ccTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa      <  1:3643457/71‑1 (MQ=255)
                                                               tGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACtt                          <  1:3791737/49‑1 (MQ=255)
                                                                 aTCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCaaa    <  1:3074562/69‑1 (MQ=255)
                                                                 aTCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAtt  >  1:3867180/1‑71 (MQ=255)
                                                                 |                                                                      
AGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATT  >  minE/2918686‑2918821

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: