Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,299,487 G→A 62.0% intergenic (‑412/+78) cobU ← / ← yeeX bifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,299,4870GA62.0% 33.4 / 47.6 50intergenic (‑412/+78)cobU/yeeXbifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/6);  new base A (15/16);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAA  >  Exported/1299419‑1299556
                                                                    |                                                                     
ttGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgaggat                                                                     >  1:2886363/1‑71 (MQ=255)
       acaacGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACg                                                              >  1:2418501/1‑71 (MQ=255)
                  gAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAg                                                   <  1:4801208/71‑1 (MQ=255)
                       aTCCGGTATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGAt                                                 >  1:736262/1‑68 (MQ=255)
                       aTCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgaggat                                                                     <  1:4214647/48‑1 (MQ=255)
                       aTCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACCGTg                                                           >  1:4438504/1‑58 (MQ=255)
                           ggAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATCATGACGGTGATACACAGATAAGAATc                                          <  1:3470771/71‑1 (MQ=255)
                                     cGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatgagcaat                                <  1:2815251/71‑1 (MQ=255)
                                       aCTATAATTTCGATTGGTTCTCGGATGAGGCtgatga                                                                <  1:734875/37‑1 (MQ=255)
                                         gttaATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGa                                                          >  1:4676791/3‑41 (MQ=255)
                                         gttaATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAg                                              <  1:1265391/51‑1 (MQ=255)
                                         gttaATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAg                                              <  1:16855/51‑1 (MQ=255)
                                         gttaATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATcaatgagca                                  >  1:4760695/3‑65 (MQ=255)
                                         gttaATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATcaatgagca                                  >  1:1578557/3‑65 (MQ=255)
                                         gttaATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAgg                            >  1:3712144/3‑71 (MQ=255)
                                          ttaATTTCGATTGGTTCTCGGATGAGAATGATGACgg                                                             <  1:1878994/36‑1 (MQ=255)
                                           taATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGt                          >  1:3682244/1‑71 (MQ=255)
                                            aaTTTCGATTGGTTCTCGGATGAGAATGATGCCGGTg                                                           >  1:5558383/1‑37 (MQ=255)
                                              tttCGATTGGTTCTCGGATGAGAATGATGACGGTGATa                                                        <  1:4485575/38‑1 (MQ=255)
                                              tttCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAAt                                           <  1:1576840/51‑1 (MQ=255)
                                              tttCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGc                       <  1:4329047/71‑1 (MQ=255)
                                               ttCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGcc                      >  1:1891705/1‑71 (MQ=255)
                                               ttCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGcc                      >  1:1879182/1‑71 (MQ=255)
                                               ttCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGcc                      >  1:4659974/1‑71 (MQ=255)
                                               ttCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGcc                      >  1:1138704/1‑71 (MQ=255)
                                                tCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGa                                                  >  1:5229794/1‑42 (MQ=255)
                                                tCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGa                                                  >  1:4813716/1‑42 (MQ=255)
                                                tCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGa                                                  >  1:1811707/1‑42 (MQ=255)
                                                 cGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATcaatgagcaatga                              >  1:3828495/1‑61 (MQ=255)
                                                   aTTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGc                  <  1:3413670/71‑1 (MQ=255)
                                                   aTTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGc                  <  1:2834904/71‑1 (MQ=255)
                                                    ttGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCa                 <  1:1119972/71‑1 (MQ=255)
                                                     tGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGc                       >  1:1126321/1‑64 (MQ=255)
                                                     tGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCaa                >  1:668830/1‑71 (MQ=255)
                                                     tGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCa                 <  1:4331363/70‑1 (MQ=255)
                                                      ggTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCaa                <  1:3194197/70‑1 (MQ=255)
                                                      ggTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCaa                <  1:776394/70‑1 (MQ=255)
                                                      ggTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATcaatgagcaatgag                             >  1:2492974/1‑57 (MQ=255)
                                                      ggTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATcaatgagcaatgag                             >  1:3507679/1‑57 (MQ=255)
                                                      ggTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATcaatgagcaatgag                             >  1:3660162/1‑57 (MQ=255)
                                                      ggTTCTCGGATGAGAATGATGACGGTGATACAAAGATAAGAATcaatgagcaatgag                             >  1:4400970/1‑57 (MQ=255)
                                                        ttCTCGGATGAGAATGATGACGGTGATACACAGATAAGAAt                                           <  1:5176984/41‑1 (MQ=255)
                                                        ttCTCGGATGAGAATGATGACGGTGATACACAGATAAGAAt                                           <  1:3817267/41‑1 (MQ=255)
                                                         tctcGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCt            >  1:4965426/1‑71 (MQ=255)
                                                         tctcGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCt            >  1:1631714/1‑71 (MQ=255)
                                                          ctcGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCt            >  1:4603600/1‑70 (MQ=255)
                                                          ctcGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCt            <  1:325468/70‑1 (MQ=255)
                                                            cGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatgagcaa                                 >  1:275216/1‑47 (MQ=255)
                                                            cGGATGAGAATGATGACGGTGATACACAGATAAGAATc                                          <  1:2150772/38‑1 (MQ=255)
                                                            cGGATGAGAATGATGACGGTGATACACAGATAAGAATc                                          <  1:2144172/38‑1 (MQ=255)
                                                            cGGATGAGAATGATGACGGTGATACACAGATAAGAATc                                          <  1:193263/38‑1 (MQ=255)
                                                            cGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCg                     >  1:3113920/1‑59 (MQ=255)
                                                              gatgaggatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCAtt        <  1:1463320/70‑1 (MQ=255)
                                                              gatgagaatgaTGACGGTGATACACAGATAAGAATcaatgagcaa                                 <  1:3938001/45‑1 (MQ=255)
                                                              gatgagaatgaTGACGGTGATACACAGATAAGAATcaatgagcaa                                 <  1:778066/45‑1 (MQ=255)
                                                              gatgagaatgaTGACGGTGATACACAGATAAGAATcaatgagcaa                                 <  1:2539479/45‑1 (MQ=255)
                                                              gatgagaatgaTGACGGTGATACACAGATAAGAATcaatgagcaa                                 <  1:4382338/45‑1 (MQ=255)
                                                              catgagaatgaTGACGGTGATACACAGATAAGAATcaatgagcaa                                 <  1:5486510/44‑1 (MQ=255)
                                                               atgagaatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGc                       <  1:3144439/54‑1 (MQ=255)
                                                               atgagaatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGa      <  1:2588942/71‑1 (MQ=255)
                                                                tgagaatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGaa     <  1:3050738/71‑1 (MQ=255)
                                                                tgagaatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGaa     <  1:1902230/71‑1 (MQ=255)
                                                                tgagaatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGa      <  1:2489490/70‑1 (MQ=255)
                                                                 gagaatgaTGACGGTGATACACAGATAAGAATcaa                                        >  1:1941911/1‑35 (MQ=255)
                                                                 gagaatgaTGACGGTGATACACAGATAAGAATcaa                                        >  1:2783634/1‑35 (MQ=255)
                                                                 gagaatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTg                        <  1:2689362/51‑1 (MQ=255)
                                                                 gagaatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGaa     >  1:867047/1‑70 (MQ=255)
                                                                   ggatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATaa  <  1:2367832/71‑1 (MQ=255)
                                                                    |                                                                     
TTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAA  >  Exported/1299419‑1299556

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: