Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,402 A→C 7.9% G98G (GGT→GGG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4020AC 7.9% 214.2 / 13.5 89G98G (GGT→GGGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (38/44);  new base C (5/2);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

TAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAA  >  minE/2919332‑2919471
                                                                      |                                                                     
tAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                                       >  1:2230889/1‑71 (MQ=255)
  gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                                       >  1:1440037/1‑69 (MQ=255)
  gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                                     <  1:1156042/71‑1 (MQ=255)
  gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAc                                                                      >  1:54930/1‑70 (MQ=255)
  gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAc                                                                      >  1:3137877/1‑70 (MQ=255)
  gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAc                                                                      >  1:2425850/1‑70 (MQ=255)
  gCGTATTCCGCCAGCGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                                       >  1:3579612/1‑69 (MQ=255)
        tgcGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAAc                                                               >  1:3176586/3‑71 (MQ=255)
         gcGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAg                                                                  >  1:2482269/2‑67 (MQ=255)
          cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                                     <  1:3164253/63‑1 (MQ=255)
          cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                                              <  1:1019544/70‑1 (MQ=255)
          cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                                              <  1:1708370/70‑1 (MQ=255)
          cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                                              <  1:1303261/70‑1 (MQ=255)
           gCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcgcg                                                            <  1:2143831/71‑1 (MQ=255)
             cAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCggg                                                          <  1:1226227/71‑1 (MQ=255)
               gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCa                                                        <  1:3019389/71‑1 (MQ=255)
               gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCa                                                        <  1:2081044/71‑1 (MQ=255)
               gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCATCATACATGCGACCAAGCACGCGGGCa                                                        <  1:1911017/71‑1 (MQ=255)
               gcgtTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAACCCgcgcg                                                        <  1:778990/69‑2 (MQ=255)
                 gtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgt                                                      >  1:576050/1‑71 (MQ=255)
                  tttCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtg                                                     <  1:1409972/71‑1 (MQ=255)
                   ttCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtg                                                     >  1:157666/1‑70 (MQ=255)
                    tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgt                                                      >  1:2757035/1‑68 (MQ=255)
                     cGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                                       >  1:410678/1‑50 (MQ=255)
                      gACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                                     <  1:3801225/51‑1 (MQ=255)
                       aCAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAc                                                                      >  1:1237543/1‑49 (MQ=255)
                       aCAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCtt                                                 >  1:3666389/1‑70 (MQ=255)
                         aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAc                                                               >  1:257113/1‑54 (MQ=255)
                         aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCttt                                                >  1:937260/1‑69 (MQ=255)
                         aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                              <  1:312579/71‑1 (MQ=255)
                         aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                              >  1:3046045/1‑71 (MQ=255)
                         aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                              <  1:492581/71‑1 (MQ=255)
                           tctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                            >  1:2754572/1‑71 (MQ=255)
                             tcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGc                                                         >  1:3645628/1‑56 (MQ=255)
                             tcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                                          >  1:1711836/1‑71 (MQ=255)
                             tcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                                          >  1:3655131/1‑71 (MQ=255)
                              ccTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAc                                         >  1:1750873/1‑71 (MQ=255)
                               cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                            <  1:2676158/67‑1 (MQ=255)
                               cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCg                                           <  1:3462241/68‑1 (MQ=255)
                               cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                                        >  1:2337930/1‑71 (MQ=255)
                               cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAc                                         >  1:1923432/1‑70 (MQ=255)
                                tGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCt                                                  <  1:1139500/60‑1 (MQ=255)
                                  aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCt                                                  >  1:1550077/1‑58 (MQ=255)
                                  aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActc                                       <  1:342752/69‑1 (MQ=255)
                                  aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                     >  1:3646189/1‑71 (MQ=255)
                                  aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                     >  1:1105342/1‑71 (MQ=255)
                                   ccATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCttt                                    <  1:1678526/71‑1 (MQ=255)
                                    cATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCtt                                                 <  1:282461/57‑1 (MQ=255)
                                     aTAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                  <  1:1988056/71‑1 (MQ=255)
                                      tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                  <  1:3021845/70‑1 (MQ=255)
                                      tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                  >  1:1231583/1‑70 (MQ=255)
                                      tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                                 <  1:3527365/71‑1 (MQ=255)
                                      tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                                 <  1:1195531/71‑1 (MQ=255)
                                      tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTCTCTTTAATCGACTCTTTAt                                  >  1:2264886/1‑70 (MQ=255)
                                       aGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTa                                               <  1:2481641/56‑1 (MQ=255)
                                         ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGa                                >  1:2516378/1‑69 (MQ=255)
                                         ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc                               >  1:287660/1‑70 (MQ=255)
                                         ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc                               >  1:3877993/1‑70 (MQ=255)
                                          cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAcc                              >  1:3140073/1‑70 (MQ=255)
                                            cgATACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAAt                           >  1:3343308/1‑71 (MQ=255)
                                            cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt                           <  1:1072865/71‑1 (MQ=255)
                                            cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt                           <  1:1501430/71‑1 (MQ=255)
                                              aTACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt                         <  1:3140774/71‑1 (MQ=255)
                                               tACTTAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt                         <  1:3182861/70‑1 (MQ=255)
                                               tACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTg                        >  1:1151001/1‑71 (MQ=255)
                                                   gAAAACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAc                                         <  1:44526/50‑1 (MQ=255)
                                                    aaaaCCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                    <  1:2421608/67‑1 (MQ=255)
                                                    aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTgg                       >  1:1196747/1‑67 (MQ=255)
                                                    aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                    <  1:3400262/70‑1 (MQ=255)
                                                    aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                    <  1:3576235/70‑1 (MQ=255)
                                                    aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                    <  1:1588042/70‑1 (MQ=255)
                                                     aTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                  >  1:3646598/1‑55 (MQ=255)
                                                      tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                                 <  1:1969164/55‑1 (MQ=255)
                                                      tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg                 >  1:2875305/1‑71 (MQ=255)
                                                      tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg                 <  1:1014431/71‑1 (MQ=255)
                                                        ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCt               <  1:530444/71‑1 (MQ=255)
                                                         cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt              >  1:390611/1‑71 (MQ=255)
                                                         cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt              <  1:3216134/71‑1 (MQ=255)
                                                           tCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTgg            <  1:3916074/71‑1 (MQ=255)
                                                           tCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCgg            >  1:3047201/1‑71 (MQ=255)
                                                              tACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc         <  1:3711769/71‑1 (MQ=255)
                                                                cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                     <  1:1304227/41‑1 (MQ=255)
                                                                cATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGcccta       <  1:1101366/71‑3 (MQ=255)
                                                                 aTGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc         >  1:2753508/1‑68 (MQ=255)
                                                                 aTGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc         >  1:3325441/1‑68 (MQ=255)
                                                                   gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTa                                   >  1:64676/1‑40 (MQ=255)
                                                                   gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga     >  1:1671216/1‑70 (MQ=255)
                                                                   gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga     >  1:2474233/1‑70 (MQ=255)
                                                                   gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga     >  1:286668/1‑70 (MQ=255)
                                                                    cGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGAt    >  1:628717/1‑70 (MQ=255)
                                                                     gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                  <  1:3747826/39‑1 (MQ=255)
                                                                     gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa  <  1:1545066/71‑1 (MQ=255)
                                                                     gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa  <  1:3385590/71‑1 (MQ=255)
                                                                      |                                                                     
TAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAA  >  minE/2919332‑2919471

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: