Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,405 A→G 7.8% L97L (CTT→CTC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4050AG 7.8% 250.6 / 9.7 89L97L (CTT→CTCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (39/43);  new base G (5/2);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTA  >  minE/2919340‑2919474
                                                                 |                                                                     
tgcGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAAc                                                                  >  1:3176586/3‑71 (MQ=255)
 gcGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAg                                                                     >  1:2482269/2‑67 (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:3655131/1‑71 (MQ=255)
                     tcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                                             >  1:1711836/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: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: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: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:3385590/71‑1 (MQ=255)
                                                             gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa     <  1:1545066/71‑1 (MQ=255)
                                                               ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGcc                     <  1:1916688/53‑1 (MQ=255)
                                                               ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg    >  1:2606645/1‑70 (MQ=255)
                                                               ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg    >  1:1435759/1‑70 (MQ=255)
                                                               ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg    >  1:3744153/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)
                                                                 |                                                                     
TCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTA  >  minE/2919340‑2919474

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: