Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,018 A→C 11.5% A226A (GCT→GCG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,0180AC11.5% 249.7 / 31.7 113A226A (GCT→GCGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (56/44);  new base C (6/7);  total (62/51)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 117 total aligned reads displayed.

CAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGC  >  minE/2918951‑2919088
                                                                   |                                                                      
caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCa                                                                     >  1:893447/1‑71 (MQ=255)
   tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCAcc                                                                   >  1:3408915/1‑70 (MQ=255)
   tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCAcc                                                                   >  1:1399315/1‑70 (MQ=255)
   tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCt                                                                  <  1:1084522/71‑1 (MQ=255)
   tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCt                                                                  <  1:3320291/71‑1 (MQ=255)
   tCCGCTCTGCCCACTTCTCTCTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGccc                                                                  >  1:2449673/1‑70 (MQ=255)
     cGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTc                                                                <  1:657793/71‑1 (MQ=255)
     cGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTc                                                                <  1:1222338/71‑1 (MQ=255)
         ccGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACt                                                            >  1:3549542/1‑71 (MQ=255)
         ccGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACt                                                            >  1:3360828/1‑71 (MQ=255)
          tgCCCACTTCTCTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGCCCTTAAc                                                             <  1:2711405/68‑1 (MQ=255)
          tgCCCACTTCTCTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGCCCTTAAc                                                             <  1:1736907/68‑1 (MQ=255)
          cGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAAc                                                             <  1:791799/69‑1 (MQ=255)
          cGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACt                                                            <  1:82081/70‑1 (MQ=255)
           gCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTcc                                                          >  1:303929/1‑71 (MQ=255)
           gCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTcc                                                          >  1:2737717/1‑71 (MQ=255)
           gCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTcc                                                          >  1:2096706/1‑71 (MQ=255)
              cATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCtt                                                       >  1:1258520/1‑71 (MQ=255)
              cATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCtt                                                       <  1:2408289/71‑1 (MQ=255)
               aTTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTc                                                      >  1:865566/1‑71 (MQ=255)
                  tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAg                                                                       >  1:3644837/1‑51 (MQ=255)
                  tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc                                                    <  1:3294383/70‑1 (MQ=255)
                      ttttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGc                                                                      >  1:1968574/1‑48 (MQ=255)
                      ttttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGc                                                                      >  1:235068/1‑48 (MQ=255)
                      ttttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTc                                                      >  1:3258529/1‑64 (MQ=255)
                       tttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcg                                                     >  1:155172/1‑64 (MQ=255)
                       tttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATc                                              >  1:2393502/1‑71 (MQ=255)
                       tttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATc                                              >  1:2388753/1‑71 (MQ=255)
                        ttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATc                                              <  1:731308/70‑1 (MQ=255)
                           cTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtc                                           <  1:3931114/70‑1 (MQ=255)
                            ttCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTc                                                                <  1:791210/48‑1 (MQ=255)
                              cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc                                                    <  1:1820010/58‑1 (MQ=255)
                              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)
                              cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcag                                         <  1:2110952/69‑1 (MQ=255)
                               tcccATCGACACCCACACGTTGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTtccgt                                        <  1:3441947/68‑1 (MQ=255)
                               ccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCt                                                        >  1:326320/1‑53 (MQ=255)
                                cccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTc                                                                <  1:490392/44‑1 (MQ=255)
                                cccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCt                                                        <  1:1616654/52‑1 (MQ=255)
                                 ccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCt                                                                  <  1:3658055/41‑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)
                                    tCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGt                                  >  1:235799/1‑70 (MQ=255)
                                     cGACAGCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGAc                                                 >  1:3020854/1‑54 (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:3334330/1‑71 (MQ=255)
                                           ccACACGTCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAATTTTcc                          <  1:2060057/71‑1 (MQ=255)
                                            cacacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCtt                                            <  1:376757/52‑1 (MQ=255)
                                             acacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCcca                        >  1:1314309/1‑71 (MQ=255)
                                             acacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCcca                        >  1:3515815/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)
                                              cacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccac                       >  1:355753/1‑71 (MQ=255)
                                               acgtCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAATTTTCccgcc                      >  1:3201090/4‑71 (MQ=255)
                                               acaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGa                                                  >  1:2655714/1‑43 (MQ=255)
                                               acaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                        <  1:1967672/53‑1 (MQ=255)
                                               acaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccac                       >  1:639797/1‑70 (MQ=255)
                                                cgtCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAatt                              <  1:1630367/60‑2 (MQ=255)
                                                cgtCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAATTTTCccgccg                     <  1:860477/69‑2 (MQ=255)
                                                caTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccacca                     <  1:3099390/71‑1 (MQ=255)
                                                 aTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtca                                          <  1:1145214/49‑1 (MQ=255)
                                                 aTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCAt                    <  1:2820243/71‑1 (MQ=255)
                                                 aTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCAt                    <  1:3451319/71‑1 (MQ=255)
                                                  tCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                        <  1:1141565/50‑1 (MQ=255)
                                                   cGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                            <  1:3553276/61‑1 (MQ=255)
                                                   cGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                            <  1:3456646/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)
                                                         aGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCtt                                            >  1:2022260/1‑39 (MQ=255)
                                                         aGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCtt                                            >  1:2641522/1‑39 (MQ=255)
                                                         aGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCtt                                            >  1:2905724/1‑39 (MQ=255)
                                                          gATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCtgt            <  1:3665501/70‑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:2915911/39‑1 (MQ=255)
                                                             aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                        <  1:3244391/39‑1 (MQ=255)
                                                             aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGc                                    >  1:37007/1‑43 (MQ=255)
                                                              aaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCcca                        >  1:1276570/1‑54 (MQ=255)
                                                               aGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTc                           >  1:3414816/1‑50 (MQ=255)
                                                                gTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAggg     >  1:1434391/1‑71 (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)
                                                                   |                                                                      
CAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGC  >  minE/2918951‑2919088

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: