Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,408 C→A 6.3% V96V (GTG→GTT argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4080CA 6.3% 297.3 / 2.3 90V96V (GTG→GTTargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (40/44);  new base A (4/2);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGC  >  minE/2919340‑2919478
                                                                    |                                                                      
tgcGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAAc                                                                      >  1:3176586/3‑71 (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)
               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:492581/71‑1 (MQ=255)
                 aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                     >  1:3046045/1‑71 (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:1105342/1‑71 (MQ=255)
                          aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                            >  1:3646189/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:1231583/1‑70 (MQ=255)
                              tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                         <  1:3021845/70‑1 (MQ=255)
                              tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                                        <  1:1195531/71‑1 (MQ=255)
                              tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                                        <  1:3527365/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:3576235/70‑1 (MQ=255)
                                            aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                           <  1:3400262/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:1014431/71‑1 (MQ=255)
                                              tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg                        >  1:2875305/1‑71 (MQ=255)
                                                ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCt                      <  1:530444/71‑1 (MQ=255)
                                                 cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt                     <  1:3216134/71‑1 (MQ=255)
                                                 cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt                     >  1:390611/1‑71 (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:3325441/1‑68 (MQ=255)
                                                         aTGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc                >  1:2753508/1‑68 (MQ=255)
                                                           gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTa                                          >  1:64676/1‑40 (MQ=255)
                                                           gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga            >  1:2474233/1‑70 (MQ=255)
                                                           gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga            >  1:1671216/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)
                                                               ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGcc                         <  1:1916688/53‑1 (MQ=255)
                                                               ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg        >  1:3744153/1‑70 (MQ=255)
                                                               ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg        >  1:1435759/1‑70 (MQ=255)
                                                               ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg        >  1:2606645/1‑70 (MQ=255)
                                                               ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg        >  1:3253438/1‑70 (MQ=255)
                                                                cAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTa      <  1:1338284/71‑1 (MQ=255)
                                                                cAAGCACGCGGGCAGTGTCTTTAATAGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGt       >  1:3808790/1‑70 (MQ=255)
                                                                  aGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAc    >  1:2019871/1‑71 (MQ=255)
                                                                  aGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAc    <  1:800302/71‑1 (MQ=255)
                                                                    cACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc  <  1:198903/71‑1 (MQ=255)
                                                                    cACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcg   >  1:1024561/1‑70 (MQ=255)
                                                                    |                                                                      
TCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGC  >  minE/2919340‑2919478

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: