Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,919,784 G→C 8.3% intergenic (‑89/‑73) argI ← / → yjgD ornithine carbamoyltransferase 1/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,919,7840GC 8.3% 461.5 / 22.5 157intergenic (‑89/‑73)argI/yjgDornithine carbamoyltransferase 1/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (86/58);  new base C (7/6);  total (93/64)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 158 total aligned reads displayed.

GCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAAT  >  Exported/2919716‑2919851
                                                                    |                                                                   
gCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTc                                                                   <  1:3988443/71‑1 (MQ=255)
     aaTGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTcaca                                                                >  1:2115807/1‑69 (MQ=255)
     aaTGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTcaca                                                                >  1:3737516/1‑69 (MQ=255)
     aaTGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaaa                                                              <  1:725143/71‑1 (MQ=255)
       tGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTcaca                                                                <  1:4754993/67‑1 (MQ=255)
        gAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGt                                                           >  1:2808177/1‑71 (MQ=255)
          aTTAAAATTCACTTTATATGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTg                                                          <  1:552584/70‑1 (MQ=255)
            tAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAg                                                       >  1:4909873/1‑71 (MQ=255)
             aaaaTTCACTTTATATGTGTAATTATTCATTTGCACCCCCATTTCACAATTCTTTCTTACAAAGGTGGAgg                                                      >  1:4112108/1‑71 (MQ=255)
              aaaTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaa                                                               >  1:3431875/1‑61 (MQ=255)
              aaaTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                                     <  1:4741111/71‑1 (MQ=255)
                tttCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaaa                                                              <  1:5648880/59‑1 (MQ=255)
                aTTCACTTTATATGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTGGAGGCaa                                                   <  1:2416299/71‑1 (MQ=255)
                 ttCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                   >  1:346257/1‑70 (MQ=255)
                 ttCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                   >  1:3072625/1‑70 (MQ=255)
                 ttCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAAt                                                  >  1:5313446/1‑71 (MQ=255)
                  tCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                                            >  1:4259186/1‑60 (MQ=255)
                  tCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                   >  1:5297771/1‑69 (MQ=255)
                   cactttATATGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTGGa                                                        >  1:4876526/4‑63 (MQ=255)
                   cAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                                            >  1:5261381/1‑59 (MQ=255)
                    aaTTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTc                                               <  1:2953688/71‑1 (MQ=255)
                     aTTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGt                                                                    <  1:4381462/49‑1 (MQ=255)
                     aTTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCa                                              <  1:496189/71‑1 (MQ=255)
                     aTTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCa                                              <  1:667106/71‑1 (MQ=255)
                      tttATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTg                                                          <  1:1656882/58‑1 (MQ=255)
                      tttATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                                     >  1:4966007/1‑63 (MQ=255)
                        tataTGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTGGAGGCaaa                                                  >  1:2163919/1‑63 (MQ=255)
                         atatGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtgg                                          >  1:3991974/1‑71 (MQ=255)
                          tatGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtggt                                         <  1:4119257/71‑1 (MQ=255)
                              gatgatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                                            >  1:1215417/1‑48 (MQ=255)
                               atgatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGgtgtgt                                     >  1:525556/1‑70 (MQ=255)
                                 gatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                                            >  1:1964015/1‑45 (MQ=255)
                                  attattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTCGTGTGTGAc                                  >  1:416787/1‑70 (MQ=255)
                                  attattCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTgg                                                         >  1:4554212/1‑47 (MQ=255)
                                  attattCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTgg                                                         >  1:3582467/1‑47 (MQ=255)
                                   ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaa                                                               <  1:1363484/40‑1 (MQ=255)
                                   ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                   >  1:2205611/1‑52 (MQ=255)
                                   ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                   >  1:1820770/1‑52 (MQ=255)
                                   ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                   >  1:2698251/1‑52 (MQ=255)
                                   ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                   >  1:1234272/1‑52 (MQ=255)
                                   ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                   >  1:5496535/1‑52 (MQ=255)
                                   ttattCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTg                                                          <  1:3561749/45‑1 (MQ=255)
                                    tattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaa                                                               <  1:3648732/39‑1 (MQ=255)
                                    tattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                                     <  1:2569352/49‑1 (MQ=255)
                                      ttCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAt                               >  1:2149397/1‑69 (MQ=255)
                                       tCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtg                                           >  1:3078145/1‑56 (MQ=255)
                                       tCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAt                               >  1:5002463/1‑68 (MQ=255)
                                       tCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAt                               >  1:1000203/1‑68 (MQ=255)
                                         aTTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                          >  1:3890303/1‑71 (MQ=255)
                                         aTTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                          >  1:851574/1‑71 (MQ=255)
                                           ttGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTg                                                          >  1:3220703/1‑37 (MQ=255)
                                            tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgt                                            >  1:3778628/1‑50 (MQ=255)
                                            tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTAt                       >  1:5344715/1‑71 (MQ=255)
                                            tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTAt                       >  1:2416517/1‑71 (MQ=255)
                                            tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAAAAGAGTAt                       >  1:4699726/1‑71 (MQ=255)
                                             gCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTg                                                          <  1:2487863/35‑1 (MQ=255)
                                              cAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                   <  1:2657415/57‑1 (MQ=255)
                                               aaGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGgtgtgt                                     <  1:3225505/54‑1 (MQ=255)
                                               aaGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCgg                    >  1:3692972/1‑71 (MQ=255)
                                                aGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGc                   <  1:5530268/71‑1 (MQ=255)
                                                 gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                                     <  1:1558451/36‑1 (MQ=255)
                                                 gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                                     <  1:4891607/36‑1 (MQ=255)
                                                 gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGAc                                  <  1:1005353/55‑1 (MQ=255)
                                                   cTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAg                 >  1:2026998/1‑70 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa               >  1:4540978/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa               >  1:299969/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa               >  1:5600443/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa               >  1:1160747/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa               >  1:4820641/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa               >  1:1570164/1‑71 (MQ=255)
                                                      aaGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGAc              <  1:4459273/70‑1 (MQ=255)
                                                      aaGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACa             >  1:5501872/1‑71 (MQ=255)
                                                       aGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCa                  >  1:3494638/1‑65 (MQ=255)
                                                       aGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACAt            <  1:4154781/71‑1 (MQ=255)
                                                       aGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACAt            <  1:1517455/71‑1 (MQ=255)
                                                           tAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgt                                            <  1:2837860/35‑1 (MQ=255)
                                                           tAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAt                               <  1:1136269/48‑1 (MQ=255)
                                                           tAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                          >  1:4209560/1‑53 (MQ=255)
                                                            aaaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                   <  1:4790631/43‑1 (MQ=255)
                                                            aaaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTaa         >  1:3393215/1‑69 (MQ=255)
                                                             aaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTaa         >  1:615444/1‑68 (MQ=255)
                                                             aaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAagag      <  1:4443757/71‑1 (MQ=255)
                                                             aaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAagag      <  1:1123610/71‑1 (MQ=255)
                                                             aaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAaga       <  1:1972154/70‑1 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGAc                                  >  1:5053565/1‑42 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAagag      >  1:1193032/1‑70 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAagag      >  1:3029495/1‑70 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAgg     <  1:1066473/71‑1 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAgg     <  1:1490840/71‑1 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAgg     <  1:1880655/71‑1 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                   >  1:3165966/1‑39 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                   >  1:1345603/1‑39 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                   >  1:5594556/1‑39 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACa                                 >  1:351258/1‑41 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGc                   <  1:2594342/55‑1 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGaa   >  1:1812894/1‑71 (MQ=255)
                                                                   tGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                          >  1:730605/1‑45 (MQ=255)
                                                                   tGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTa                        >  1:3804617/1‑47 (MQ=255)
                                                                   tGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTa                        >  1:3155701/1‑47 (MQ=255)
                                                                   tGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTaa         <  1:5169190/62‑1 (MQ=255)
                                                                    |                                                                   
GCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAAT  >  Exported/2919716‑2919851

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: