Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,045,818 C→T 6.5% intergenic (+580/+215) yqiK → / ← rfaE conserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,045,8180CT 6.5% 210.4 / 3.9 76intergenic (+580/+215)yqiK/rfaEconserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base C (31/40);  new base T (2/3);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

AAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGT  >  Exported/2045750‑2045887
                                                                    |                                                                      
aaaGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCt                                                                      <  1:1302570/71‑1 (MQ=255)
      ttGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaa                                                                <  1:4422296/71‑1 (MQ=255)
       tGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCg                                                                  >  1:565428/1‑68 (MQ=255)
         cgcgGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaagaa                                                             <  1:3636028/71‑1 (MQ=255)
         cgcgGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag                                                               >  1:1446895/1‑69 (MQ=255)
               gTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGcc                                                      >  1:4650558/1‑71 (MQ=255)
                  cGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                              >  1:4429500/1‑61 (MQ=255)
                      gcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGta                                                 >  1:5442407/1‑69 (MQ=255)
                       cTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtata                                               <  1:1887124/70‑1 (MQ=255)
                          atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                <  1:2333700/66‑1 (MQ=255)
                          atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                            >  1:3343442/1‑70 (MQ=255)
                          atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAg                                             <  1:5466060/69‑1 (MQ=255)
                           tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTg                                                   <  1:4827154/62‑1 (MQ=255)
                           tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTg                                                   <  1:1552405/62‑1 (MQ=255)
                           tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGaa                                          <  1:2717963/71‑1 (MQ=255)
                             aacaaTCGGGGGGCAGCAAGGGGCTGAAACTGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGa                                           <  1:892892/68‑1 (MQ=255)
                             aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag                                                               >  1:822613/1‑49 (MQ=255)
                             aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAg                                        >  1:2466718/1‑71 (MQ=255)
                             aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAg                                        >  1:5640262/1‑71 (MQ=255)
                              ctaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGaaa                                         <  1:4132819/67‑1 (MQ=255)
                              acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTg                                                   <  1:5090951/59‑1 (MQ=255)
                              acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTg                                                   <  1:711585/59‑1 (MQ=255)
                              acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTg                                                   <  1:1271467/59‑1 (MQ=255)
                                     gggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag                                                               <  1:914440/41‑1 (MQ=255)
                                     gggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag                                                               <  1:3182772/41‑1 (MQ=255)
                                        ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                            <  1:4051813/56‑1 (MQ=255)
                                        ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                            <  1:276982/56‑1 (MQ=255)
                                        ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                            <  1:346786/56‑1 (MQ=255)
                                        ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                            <  1:2562837/56‑1 (MQ=255)
                                         gcagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                            <  1:165441/55‑1 (MQ=255)
                                         gcagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGa                                           <  1:733399/56‑1 (MQ=255)
                                          cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAgg                                                          >  1:2605999/1‑40 (MQ=40)
                                          cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAg                           >  1:830053/1‑71 (MQ=255)
                                           agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTa                                   <  1:4329956/62‑1 (MQ=255)
                                           agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTa                                   <  1:2531398/62‑1 (MQ=255)
                                            gcaAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                              <  1:4260276/35‑1 (MQ=40)
                                            gcaAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                              <  1:3957118/35‑1 (MQ=40)
                                            gcaAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                              <  1:4963063/35‑1 (MQ=40)
                                            gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                            <  1:4661619/52‑1 (MQ=255)
                                            gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAgg                                       <  1:3000704/57‑1 (MQ=255)
                                            gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTatg                                 >  1:4328361/1‑63 (MQ=255)
                                              aaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGt                                     >  1:781119/1‑57 (MQ=255)
                                               aGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTa                                   >  1:3736880/1‑58 (MQ=255)
                                               aGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                        <  1:657273/69‑1 (MQ=255)
                                                 gggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGTCCTTGTATAAGGAAAGGGTTa                                   <  1:1977594/56‑1 (MQ=255)
                                                 acgCTGAAACGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGt                      >  1:303292/3‑70 (MQ=0)
                                                  ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                <  1:57129/42‑1 (MQ=255)
                                                  ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                   >  1:552372/1‑71 (MQ=255)
                                                  ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                   >  1:5627124/1‑71 (MQ=255)
                                                  ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                   >  1:3187602/1‑71 (MQ=255)
                                                  cgCTGAAACGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGtca                    <  1:1786500/70‑1 (MQ=1)
                                                    cTGAAACGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGaaa                                         >  1:3653992/1‑47 (MQ=14)
                                                    cTGAAACGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCatg                                 <  1:144998/55‑1 (MQ=14)
                                                    cTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATaa                                              >  1:442617/1‑42 (MQ=255)
                                                      gAAACGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGaa                                          >  1:1492151/1‑44 (MQ=14)
                                                      gAAACGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACtcat                   >  1:1889915/1‑67 (MQ=14)
                                                      gAAACGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACtcat                   >  1:1896669/1‑67 (MQ=14)
                                                      gAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGtt                                    <  1:1190765/50‑1 (MQ=255)
                                                       aaaCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCatgat                               <  1:4403377/54‑1 (MQ=14)
                                                       aaaCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCatgat                               <  1:3787364/54‑1 (MQ=14)
                                                       aaaCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcatcg              <  1:565443/71‑3 (MQ=11)
                                                       aaaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc              >  1:4899675/1‑71 (MQ=255)
                                                        aaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAg                           >  1:3535738/1‑57 (MQ=255)
                                                         aCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCAtgatga                             >  1:3633700/1‑54 (MQ=14)
                                                         aCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCAtgatga                             >  1:3590871/1‑54 (MQ=14)
                                                         aCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAATGAAAGGGTTATGATGAAGCTCGt                      <  1:5193989/61‑1 (MQ=255)
                                                            ggAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCatgat                               <  1:2622584/49‑1 (MQ=14)
                                                            ggAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtg         >  1:1155474/1‑70 (MQ=2)
                                                            ggAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgtt         >  1:3793091/1‑71 (MQ=255)
                                                             gAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGgt          >  1:3680797/1‑70 (MQ=11)
                                                             gAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtgc        >  1:3332284/1‑69 (MQ=2)
                                                             gAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg        <  1:2081641/71‑1 (MQ=255)
                                                              aaaGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtgc        >  1:4164927/1‑68 (MQ=2)
                                                              aaaGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtgc        >  1:3132334/1‑68 (MQ=2)
                                                                aGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGt                                     <  1:3505689/40‑1 (MQ=0)
                                                                aGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGt                                     <  1:22867/40‑1 (MQ=0)
                                                                   cccTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACt             <  1:1253870/60‑1 (MQ=255)
                                                                    ccTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt  >  1:3907027/1‑70 (MQ=255)
                                                                    |                                                                      
AAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGT  >  Exported/2045750‑2045887

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: