Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,048 A→C 7.5% T216T (ACT→ACG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,0480AC 7.5% 266.5 / 16.2 106T216T (ACT→ACGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (52/46);  new base C (4/4);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 128 total aligned reads displayed.

CCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGC  >  minE/2918981‑2919117
                                                                   |                                                                     
cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                                     >  1:2863978/1‑70 (MQ=255)
cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                                     >  1:95989/1‑70 (MQ=255)
cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                                     >  1:323326/1‑70 (MQ=255)
cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                                     >  1:1219262/1‑70 (MQ=255)
 tcccATCGACACCCACACGTTGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTtccgt                                                                     <  1:3441947/68‑1 (MQ=255)
   ccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGc                                                                 >  1:3161924/1‑71 (MQ=255)
     aTCGACACCCACACGTCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGt                                                               <  1:2063317/71‑1 (MQ=255)
     aTCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGt                                                               <  1:3177021/71‑1 (MQ=255)
      tCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagtc                                                                    >  1:1165605/1‑65 (MQ=255)
        gACACCCACACGTCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAat                                                            >  1:1921459/1‑71 (MQ=255)
          cacCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAatat                                                          <  1:3711650/71‑1 (MQ=255)
           acCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAatat                                                          >  1:3902274/1‑70 (MQ=255)
           acCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAata                                                           >  1:1168117/1‑69 (MQ=255)
           acCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                                                         <  1:1168636/71‑1 (MQ=255)
           acCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                                                         >  1:3543402/1‑71 (MQ=255)
             ccACACGTCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAATTTTcc                                                       >  1:3005535/1‑71 (MQ=255)
             ccACACGTCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAATTTTcc                                                       <  1:2060057/71‑1 (MQ=255)
             ccACACGTCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAATTTTcc                                                       >  1:3334330/1‑71 (MQ=255)
               acacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCcca                                                     >  1:2699161/1‑71 (MQ=255)
               acacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCcca                                                     >  1:1935562/1‑71 (MQ=255)
                cacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccac                                                    >  1:2681229/1‑71 (MQ=255)
                 acaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                                     <  1:1967672/53‑1 (MQ=255)
                 acaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccac                                                    >  1:639797/1‑70 (MQ=255)
                  cgtCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAatt                                                           <  1:1630367/60‑2 (MQ=255)
                  caTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccacca                                                  <  1:3099390/71‑1 (MQ=255)
                  caTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccacca                                                  <  1:3089366/71‑1 (MQ=255)
                   aTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCAt                                                 <  1:2820243/71‑1 (MQ=255)
                    tCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                                     <  1:1141565/50‑1 (MQ=255)
                     cGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                                                         <  1:3456646/61‑1 (MQ=255)
                     cGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                                                         <  1:1350132/61‑1 (MQ=255)
                     cGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                                                         <  1:3553276/61‑1 (MQ=255)
                      ggTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtca                                                                       >  1:567349/1‑46 (MQ=255)
                          taGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGc                                                                 >  1:3874679/1‑48 (MQ=255)
                          taGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccacca                                                  <  1:2585125/63‑1 (MQ=255)
                             aTAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTTAGCGTAATATTCCCACCATTTTGCTGTGc                                       >  1:2889754/1‑71 (MQ=255)
                             aTAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccacc                                                   <  1:2545858/59‑1 (MQ=255)
                              tAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGcc                                      <  1:9899/71‑1 (MQ=255)
                              tAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGcc                                      <  1:2433784/71‑1 (MQ=255)
                              tAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGcc                                      <  1:3176763/71‑1 (MQ=255)
                              tAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGc                                       >  1:3497461/1‑70 (MQ=255)
                               aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                                     <  1:3239698/39‑1 (MQ=255)
                               aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                                     <  1:3244391/39‑1 (MQ=255)
                               aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                                     <  1:2915911/39‑1 (MQ=255)
                               aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGc                                                                 >  1:37007/1‑43 (MQ=255)
                                aaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCcca                                                     >  1:1276570/1‑54 (MQ=255)
                                aaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAg                                    >  1:3615354/1‑71 (MQ=255)
                                 aGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTcc                                                       >  1:3476320/1‑51 (MQ=255)
                                 aGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTc                                                        >  1:3414816/1‑50 (MQ=255)
                                  gTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAggg                                  >  1:1434391/1‑71 (MQ=255)
                                   tCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGgc                                 <  1:829308/71‑1 (MQ=255)
                                   tCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGgc                                 <  1:582572/71‑1 (MQ=255)
                                     aGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAata                                                           >  1:3084600/1‑43 (MQ=255)
                                     aGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAata                                                           >  1:1192768/1‑43 (MQ=255)
                                     aGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCAtttt                                              >  1:2586857/1‑56 (MQ=255)
                                     aGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGgcgc                               >  1:479238/1‑71 (MQ=255)
                                      gCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGgcg                                <  1:408460/69‑1 (MQ=255)
                                       cACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAat                                                            >  1:2086877/1‑40 (MQ=255)
                                       cACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCgg                             >  1:423687/1‑71 (MQ=255)
                                       cACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCgg                             >  1:2082138/1‑71 (MQ=255)
                                        aCCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCATGGCGCGGc                            <  1:2470564/71‑1 (MQ=255)
                                         ccTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCAt                                                 <  1:1363785/49‑1 (MQ=255)
                                         ccTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCgg                             <  1:157461/69‑1 (MQ=255)
                                         ccTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCa                           >  1:227766/1‑71 (MQ=255)
                                          cTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGgcgc                               >  1:1313973/1‑66 (MQ=255)
                                           ttCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGc                                            <  1:3341676/52‑1 (MQ=255)
                                           ttCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCAtt                         >  1:1657205/1‑71 (MQ=255)
                                            tCAACTCCCTTCGCGACATCTTCAGTCAGCGTAatat                                                          <  1:1694219/37‑1 (MQ=255)
                                            tCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTccc                                                      >  1:2213138/1‑41 (MQ=255)
                                            tCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTc                        >  1:1101509/1‑71 (MQ=255)
                                            acAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccacca                                                  >  1:1070518/2‑45 (MQ=255)
                                             cAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTcc                       >  1:698474/1‑71 (MQ=255)
                                              aaCTCCCTTCGCGACATCTTCAGTCAGCGTAATATTc                                                        >  1:2521250/1‑37 (MQ=255)
                                              aaCTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTcc                       >  1:1260552/1‑70 (MQ=255)
                                               aCTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGt                     <  1:3343767/71‑1 (MQ=255)
                                                cTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGcgcg                              <  1:3640766/61‑1 (MQ=255)
                                                cTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCAtt                         <  1:333341/66‑1 (MQ=255)
                                                cTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTa                    <  1:828876/71‑1 (MQ=255)
                                                cTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTa                    <  1:820045/71‑1 (MQ=255)
                                                   ccTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTaa                   <  1:2064240/69‑1 (MQ=255)
                                                   ccTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAAc                  <  1:391061/70‑1 (MQ=255)
                                                    cTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCa                <  1:2150435/71‑1 (MQ=255)
                                                    cTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCa                <  1:3624069/71‑1 (MQ=255)
                                                      tCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAgc              >  1:2740607/1‑71 (MQ=255)
                                                       cgcgACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAgc              <  1:3410719/70‑1 (MQ=255)
                                                           aCATCTTCAGTCAGCGTAATATTCCCACCATTTTGCt                                           <  1:2900852/37‑1 (MQ=255)
                                                           aCATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAAc                  >  1:2089322/1‑62 (MQ=255)
                                                            cATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCt        <  1:1289431/71‑1 (MQ=255)
                                                            cATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCt        >  1:1172774/1‑71 (MQ=255)
                                                               cTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGt                     <  1:3849767/55‑1 (MQ=255)
                                                               cTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAgcg             <  1:1109829/63‑1 (MQ=255)
                                                                ttCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAg                                    >  1:3532052/1‑39 (MQ=255)
                                                                ttCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAgcg             <  1:265656/62‑1 (MQ=255)
                                                                ttCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCg    >  1:3463171/1‑71 (MQ=255)
                                                                ttCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCg    >  1:878003/1‑71 (MQ=255)
                                                                  cagtcagCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTc      <  1:1768344/67‑1 (MQ=255)
                                                                  cagtcagCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCgg   <  1:3591527/70‑1 (MQ=255)
                                                                  cagtcagCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGc  >  1:850770/1‑71 (MQ=255)
                                                                   agtcagCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTc      >  1:3935135/1‑66 (MQ=255)
                                                                   agtcagCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTc      >  1:120434/1‑66 (MQ=255)
                                                                   agtcagCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGc  <  1:2087536/70‑1 (MQ=255)
                                                                   |                                                                     
CCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGC  >  minE/2918981‑2919117

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: