Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,030 T→G 7.3% G222G (GGA→GGC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,0300TG 7.3% 314.9 / 18.7 123G222G (GGA→GGCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (67/47);  new base G (4/5);  total (71/52)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

Only 100 of 126 total aligned reads displayed.

CCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTA  >  minE/2918960‑2919099
                                                                      |                                                                     
ccGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACt                                                                       >  1:3549542/1‑71 (MQ=255)
ccGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACt                                                                       >  1:3360828/1‑71 (MQ=255)
 cGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACt                                                                       <  1:82081/70‑1 (MQ=255)
  gCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTcc                                                                     >  1:303929/1‑71 (MQ=255)
  gCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTcc                                                                     >  1:2737717/1‑71 (MQ=255)
     cATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCtt                                                                  <  1:2408289/71‑1 (MQ=255)
     cATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCtt                                                                  >  1:1258520/1‑71 (MQ=255)
      aTTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTc                                                                 >  1:865566/1‑71 (MQ=255)
         tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc                                                               <  1:3294383/70‑1 (MQ=255)
              tttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcg                                                                >  1:155172/1‑64 (MQ=255)
              tttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATc                                                         >  1:2393502/1‑71 (MQ=255)
               ttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATc                                                         <  1:731308/70‑1 (MQ=255)
                 gCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCtt                                                       <  1:1767907/70‑1 (MQ=255)
                     cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc                                                               <  1:1820010/58‑1 (MQ=255)
                     cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                   >  1:95989/1‑70 (MQ=255)
                     cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                   >  1:2863978/1‑70 (MQ=255)
                     cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                   >  1:323326/1‑70 (MQ=255)
                     cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcag                                                    <  1:2110952/69‑1 (MQ=255)
                      ccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCt                                                                   >  1:326320/1‑53 (MQ=255)
                       cccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCt                                                                   <  1:1616654/52‑1 (MQ=255)
                        ccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGc                                               >  1:3161924/1‑71 (MQ=255)
                        ccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGc                                               >  1:2549709/1‑71 (MQ=255)
                          aTCGACACCCACACGTCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGt                                             <  1:2063317/71‑1 (MQ=255)
                           tCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGt                                             >  1:235799/1‑70 (MQ=255)
                            cGACAGCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGAc                                                            >  1:3020854/1‑54 (MQ=255)
                             gACACCCACACGTCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAat                                          >  1:1921459/1‑71 (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)
                                   cacacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCtt                                                       <  1:376757/52‑1 (MQ=255)
                                    acacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCcca                                   >  1:1314309/1‑71 (MQ=255)
                                    acacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCcca                                   >  1:2699161/1‑71 (MQ=255)
                                     cacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccac                                  >  1:2681229/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)
                                       caTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccacca                                <  1:3089366/71‑1 (MQ=255)
                                       caTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccacca                                <  1:3099390/71‑1 (MQ=255)
                                        aTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtca                                                     <  1:1145214/49‑1 (MQ=255)
                                        aTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCAt                               <  1:3451319/71‑1 (MQ=255)
                                         tCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                   <  1:1141565/50‑1 (MQ=255)
                                          cGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                                       <  1:3456646/61‑1 (MQ=255)
                                          cGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                                       <  1:3553276/61‑1 (MQ=255)
                                          cGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                                       <  1:1350132/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: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:3176763/71‑1 (MQ=255)
                                                   tAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGcc                    <  1:2433784/71‑1 (MQ=255)
                                                   tAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGcc                    <  1:9899/71‑1 (MQ=255)
                                                   tAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGc                     >  1:3497461/1‑70 (MQ=255)
                                                    aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                   <  1:2915911/39‑1 (MQ=255)
                                                    aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                   <  1:3239698/39‑1 (MQ=255)
                                                    aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                                   <  1:3244391/39‑1 (MQ=255)
                                                    aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGc                                               >  1:37007/1‑43 (MQ=255)
                                                     aaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAata                                         >  1:3357466/1‑48 (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:582572/71‑1 (MQ=255)
                                                        tCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGgc               <  1:829308/71‑1 (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:2082138/1‑71 (MQ=255)
                                                            cACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCgg           >  1:423687/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)
                                                                 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:387851/1‑71 (MQ=255)
                                                                  cAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTcc     >  1:698474/1‑71 (MQ=255)
                                                                   aaCTCCCTTCGCGACATCTTCAGTCAGCGTAATATTc                                      >  1:1259898/1‑37 (MQ=255)
                                                                   aaCTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTcc     >  1:1260552/1‑70 (MQ=255)
                                                                   aaCTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTcc     >  1:1183515/1‑70 (MQ=255)
                                                                    aCTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCCTTCCGt   <  1:195/71‑1 (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:2075453/71‑1 (MQ=255)
                                                                     cTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTa  <  1:820045/71‑1 (MQ=255)
                                                                     cTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTa  <  1:828876/71‑1 (MQ=255)
                                                                      |                                                                     
CCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTA  >  minE/2918960‑2919099

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: