Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,589 C→T 14.4% S36N (AGC→AAC) ‡ argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,5890CT14.4% 207.9 / 22.7 83S36N (AGC→AAC) ‡argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (25/46);  new base T (3/9);  total (28/55)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAG  >  minE/2919521‑2919659
                                                                    |                                                                      
tACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                                                      <  1:2555359/71‑1 (MQ=255)
tACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                                                      <  1:3063819/71‑1 (MQ=255)
  cgagtcgagtCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                                                     <  1:1964499/70‑1 (MQ=255)
  cgagtcgagtCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                                                     >  1:1144197/1‑70 (MQ=255)
     gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                                                 <  1:156901/71‑1 (MQ=255)
     gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                                                 <  1:458861/71‑1 (MQ=255)
     gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                                                 <  1:3491272/71‑1 (MQ=255)
     gtcgagtcTTTTTCGAAGATGAGCCCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                                                  >  1:1430430/1‑70 (MQ=255)
      tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTa                                                                <  1:2460278/71‑1 (MQ=255)
       cgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                                                      <  1:932425/64‑1 (MQ=255)
        gagtcTTTTTGGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATc                                                              >  1:2406720/1‑71 (MQ=255)
        gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGATTCTGTACTTCCTTGCCATTTTTTTTATc                                                              >  1:1382560/1‑71 (MQ=255)
        gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTAt                                                               <  1:1185952/70‑1 (MQ=255)
        gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                                              <  1:917204/71‑1 (MQ=255)
          gtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCg                                                             >  1:737502/1‑70 (MQ=255)
          gtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCgg                                                            <  1:1258608/71‑1 (MQ=255)
           tcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGc                                                           <  1:3193547/71‑1 (MQ=255)
             tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTc                                                                   >  1:1812235/1‑61 (MQ=255)
               tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtt                                                         <  1:2106005/69‑1 (MQ=255)
               tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtt                                                         <  1:3306119/69‑1 (MQ=255)
               tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCt                                                          <  1:260631/68‑1 (MQ=255)
               tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCt                                                          <  1:2574720/68‑1 (MQ=255)
               tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCt                                                          <  1:3732892/68‑1 (MQ=255)
               tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                                                        <  1:1762441/70‑1 (MQ=255)
               tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                                                        <  1:1167246/70‑1 (MQ=255)
                ttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtttg                                                       <  1:875383/70‑2 (MQ=255)
                ttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTc                                                       <  1:918403/70‑1 (MQ=255)
                   gAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                                                     <  1:2989852/53‑1 (MQ=255)
                   gAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                                              >  1:2414015/1‑60 (MQ=255)
                      gATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTc                                                                   >  1:3413156/1‑52 (MQ=255)
                      gATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcg                                                >  1:1783812/1‑71 (MQ=255)
                      gATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcg                                                <  1:1160/71‑1 (MQ=255)
                       aTGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                                         <  1:3525389/61‑1 (MQ=255)
                       aTGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcgc                                               <  1:1868592/71‑1 (MQ=255)
                       aTGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcgc                                               >  1:104555/1‑71 (MQ=255)
                        tGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGc                                                           <  1:1810735/58‑1 (MQ=255)
                        tGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtttt                                                                    >  1:3398013/1‑49 (MQ=255)
                         gAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                                                     >  1:79207/1‑47 (MQ=255)
                         gAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                                                     >  1:2509847/1‑47 (MQ=255)
                          aGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGc                                            >  1:89350/1‑71 (MQ=255)
                           gcgcGATGTTTTTACCAGTGTGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                                           <  1:2880690/71‑1 (MQ=255)
                           gcgcGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                                           <  1:2816680/71‑1 (MQ=255)
                           gcgcGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                                           <  1:2324083/71‑1 (MQ=255)
                            cgcgATGTTTTTAGCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTa                                          <  1:1296647/71‑1 (MQ=255)
                            cgcgATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGcaa                                          <  1:3578651/71‑3 (MQ=255)
                              cgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTc                                                       >  1:726269/1‑56 (MQ=255)
                                aTGTTTTTACCAGTGAGGTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTa                                          <  1:728554/67‑1 (MQ=255)
                                 tGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCg                                                             >  1:3055440/1‑47 (MQ=255)
                                 tGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCg                                                             >  1:194072/1‑47 (MQ=255)
                                  gTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                                                        >  1:1636867/1‑51 (MQ=255)
                                   tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                                                  <  1:2845223/40‑1 (MQ=255)
                                   tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                                                  <  1:792723/40‑1 (MQ=255)
                                   tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                                                  <  1:451192/40‑1 (MQ=255)
                                   tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCa                                                      <  1:1537903/52‑1 (MQ=255)
                                   tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag                                   <  1:2830672/71‑1 (MQ=255)
                                   tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgca                                    <  1:2776959/70‑1 (MQ=255)
                                     tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag                                   >  1:455139/1‑69 (MQ=255)
                                          cAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTg                            <  1:1053806/71‑1 (MQ=255)
                                            gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt                          <  1:2841737/71‑1 (MQ=255)
                                            gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt                          <  1:857016/71‑1 (MQ=255)
                                             tGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt                          >  1:1997174/1‑70 (MQ=255)
                                             tGAGTTTGGCTTCTTCTTTACCGCTTTTATTATCGGCTTTCAGCTTCGCGGCt                                           >  1:2001694/1‑53 (MQ=255)
                                              gAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt                          >  1:3628590/1‑69 (MQ=255)
                                               aGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg                         <  1:683609/69‑1 (MQ=255)
                                                gTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTct                                                <  1:111335/45‑2 (MQ=255)
                                                      cttcttcttTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTc                    <  1:1127622/67‑1 (MQ=255)
                                                       ttcttcttTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                                           <  1:429318/43‑1 (MQ=255)
                                                       ttcttcttTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagc                                  <  1:2887964/52‑1 (MQ=255)
                                                       ttcttcttTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGt                           <  1:1969946/59‑1 (MQ=255)
                                                        tcttcttTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt                          >  1:327109/1‑59 (MQ=255)
                                                         cttcttTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTg               <  1:1344883/69‑1 (MQ=255)
                                                         cttcttTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGc             >  1:3092359/1‑71 (MQ=255)
                                                         cttcttTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTCTTGAGTTCAGCTgg              <  1:3197157/70‑1 (MQ=255)
                                                         cttcttTACCGCTTTTCTTATCGGCTTTCAACTTCGCGGCTaa                                         >  1:3609696/1‑43 (MQ=255)
                                                             tttACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGcagcag                                <  1:1546572/48‑1 (MQ=255)
                                                               tACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGaaa       >  1:2900650/1‑71 (MQ=255)
                                                               tACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGaaa       >  1:3251835/1‑71 (MQ=255)
                                                               tACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGaa        <  1:2746118/70‑1 (MQ=255)
                                                                aCCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGaaa       <  1:391322/70‑1 (MQ=255)
                                                                 ccGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATc     <  1:840109/71‑1 (MQ=255)
                                                                 ccGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATc     >  1:120022/1‑71 (MQ=255)
                                                                 ccGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATc     <  1:211394/71‑1 (MQ=255)
                                                                   gCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAg                       <  1:2384591/51‑1 (MQ=255)
                                                                    cTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAg  <  1:714676/71‑1 (MQ=255)
                                                                    cTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAg  >  1:243207/1‑71 (MQ=255)
                                                                    |                                                                      
TACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAG  >  minE/2919521‑2919659

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