Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,581 A→G 13.6% intergenic (‑309/+142) rhaD ← / ← rhaA rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,5810AG13.6% 60.6 / 4.5 33intergenic (‑309/+142)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base G (3/2);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

TCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGT  >  NC_000913/4094444‑4094722
                                                                                                                                         |                                                                                                                                             
tcagtGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGcc                                                                                                                                          <  1:310020/143‑1 (MQ=21)
tcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGcc                                                                                                                                          <  1:257180/143‑1 (MQ=35)
  aatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGa                                                                                                                                        >  1:94379/1‑143 (MQ=35)
     gCCGGATGCGACGCTGTCGCGGCTTATCAGGCCGACAACTGTGCCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCt                                                                                                                                     <  2:485361/143‑1 (MQ=11)
         ggggCGACGCTGGCGCGGCGTATCGGGCCTACAACGGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAg                                                                                                                                 <  1:26760/140‑1 (MQ=11)
                cGCTGTCGCGTCGGATCAGGCCTACAACGGTGGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGAt                                                                                                                          <  2:63053/143‑1 (MQ=255)
                                              ttGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatcaat                                                                                            <  1:371898/143‑1 (MQ=34)
                                                 ccGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGcc                                                                                         >  1:458325/1‑143 (MQ=35)
                                                                      gCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCCCTTGTGCCGCATCCGGGAATCAATATCTTATGCGCCGCTGACGCGTcc                                                                    >  1:222776/1‑142 (MQ=11)
                                                                           attgCCGGATGCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCa                                                               >  1:288860/3‑143 (MQ=12)
                                                                                  cacccGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                        <  1:83949/140‑1 (MQ=17)
                                                                                          aatcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                                                                                            <  2:333637/83‑1 (MQ=33)
                                                                                                  ccTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCCTCCTATCAGGCCCACAACTGTTGCCGCCCCg                                        >  2:443095/1‑143 (MQ=11)
                                                                                                   cTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGt                                       <  2:369583/143‑1 (MQ=34)
                                                                                                        gCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCCACAACTGTTGCCGACCCGTCGGcc                                  >  2:131845/1‑143 (MQ=25)
                                                                                                         cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTCGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCg                                 <  2:168581/143‑1 (MQ=25)
                                                                                                          gACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCgg                                >  1:104685/1‑143 (MQ=35)
                                                                                                           ccGCTGTCCCGTCTTATCAGGCCTCCAACCATAGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGa                               <  2:109069/142‑1 (MQ=255)
                                                                                                              cTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATaa                            <  1:475095/143‑1 (MQ=35)
                                                                                                                gTCGCGTCGTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAgg                          <  1:353019/143‑1 (MQ=21)
                                                                                                                    cgTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCGCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGCCTTCTCAGGCCCACCACTGTTCCCCACCCGTCGGCCGGCTGAGgcccc                      >  2:449203/1‑140 (MQ=11)
                                                                                                                    agTCTTCCCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCt                      <  1:70280/142‑1 (MQ=255)
                                                                                                                      tCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCCGGCGGACAACTGTTGCCGACCCGTCGGCCGGATAAGgcgcctg                    >  2:111599/1‑140 (MQ=11)
                                                                                                                         tATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATACGGCGCTCgcgc                 >  1:112976/1‑143 (MQ=255)
                                                                                                                          aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTCCAACTGTTGCCGACCCGTGGGCCGGATAAGGCTGTCgcgcg                >  2:405319/1‑142 (MQ=11)
                                                                                                                          aTCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgcc                >  1:359073/1‑143 (MQ=255)
                                                                                                                           tCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgccg               <  2:419209/143‑1 (MQ=255)
                                                                                                                           tCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgccg               <  2:419177/143‑1 (MQ=255)
                                                                                                                                 cTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATCAGGCGCTCGCGCCGCCTCCg         >  2:123053/1‑143 (MQ=255)
                                                                                                                                 cTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCACTCGTCAATCAATGCATGATGGGGAGCGGTTGCGTCGTATCAGGCCTACAACTGTAGCCCAGCCGCAAACCGGATAAGGGGCTCGCGCCGCCTCCg         >  2:63951/1‑143 (MQ=14)
                                                                                                                                    cAACTGTTGCCGCGCTGTAGGCCGGATAAGGCACTCGTGCCGCATCCGGCAAGCAATGCCGGATGGGACGCTGTCGCGTCTTATCAGGCCTACAACGGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCa      <  2:275873/143‑1 (MQ=14)
                                                                                                                                     aaCTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAg     >  1:166232/1‑143 (MQ=18)
                                                                                                                                        tATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGAGCCGCATCCCGCAgtgc  >  2:57742/1‑142 (MQ=255)
                                                                                                                                         |                                                                                                                                             
TCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGT  >  NC_000913/4094444‑4094722

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