Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,474,881 G→A 30.8% intergenic (+25/‑103) yfdQ → / → yfdR CPS‑53 (KpLE1) prophage; uncharacterized protein/CPS‑53 (KpLE1) prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,474,8810GA30.8% 36.6 / 19.1 26intergenic (+25/‑103)yfdQ/yfdRCPS‑53 (KpLE1) prophage; uncharacterized protein/CPS‑53 (KpLE1) prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base G (9/9);  new base A (4/4);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAA  >  NC_000913/2474735‑2474983
                                                                                                                                                  |                                                                                                      
ttCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCgg                                                                                                       <  2:582166/148‑1 (MQ=255)
                                       tGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCa                                                                                                  >  4:582593/1‑114 (MQ=255)
                                       tGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCa                                                                                                  <  3:582593/114‑1 (MQ=255)
                                            aaCGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCg                                                                                                        <  4:955844/103‑1 (MQ=255)
                                            aaCGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCg                                                                                                        >  3:955844/1‑103 (MQ=255)
                                                         aTCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCttt                                              >  1:1094714/1‑148 (MQ=255)
                                                           cTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTa                                            <  2:1094714/148‑1 (MQ=255)
                                                              cTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAt                                                                                                 >  4:38510/1‑92 (MQ=255)
                                                              cTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAt                                                                                                 <  3:38510/92‑1 (MQ=255)
                                                                          ttCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCggggg                                                                                                    <  3:109518/77‑1 (MQ=255)
                                                                          ttCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCggggg                                                                                                    >  4:109518/1‑77 (MQ=255)
                                                                               aGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCa                                                                    >  4:750011/1‑104 (MQ=255)
                                                                               aGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCa                                                                    <  3:750011/104‑1 (MQ=255)
                                                                                    gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                                  >  4:555003/1‑85 (MQ=255)
                                                                                    gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                                  <  3:555003/85‑1 (MQ=255)
                                                                                                ccTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtt                                                                                                >  4:241476/1‑59 (MQ=38)
                                                                                                ccTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtt                                                                                                <  3:241476/59‑1 (MQ=38)
                                                                                                   ttATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCtggg                                                                                                     >  2:362482/1‑51 (MQ=255)
                                                                                                   ttATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCgggg                                                                                                     <  1:362482/51‑1 (MQ=255)
                                                                                                              ttCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATc                                                                     <  1:1126152/72‑1 (MQ=255)
                                                                                                              ttCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATc                                                                     >  2:1126152/1‑72 (MQ=255)
                                                                                                                  ccGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAt                                                                                                 >  3:816827/1‑40 (MQ=25)
                                                                                                                  ccGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAt                                                                                                 <  4:816827/40‑1 (MQ=25)
                                                                                                                    gCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATAc              >  3:293433/1‑121 (MQ=255)
                                                                                                                    gCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATAc              <  4:293433/121‑1 (MQ=255)
                                                                                                                          tttCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                       <  1:266142/106‑1 (MQ=255)
                                                                                                                          tttCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                       >  2:266142/1‑106 (MQ=255)
                                                                                                                               ttACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCt                                                >  2:86921/1‑76 (MQ=255)
                                                                                                                               ttACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCt                                                <  1:86921/76‑1 (MQ=255)
                                                                                                                                tACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                       >  3:311635/1‑100 (MQ=255)
                                                                                                                                tACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                       <  4:311635/100‑1 (MQ=255)
                                                                                                                                              cTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACaaa  >  2:1003057/1‑107 (MQ=255)
                                                                                                                                              cTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACaaa  <  1:1003057/107‑1 (MQ=255)
                                                                                                                                                  |                                                                                                      
TTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAA  >  NC_000913/2474735‑2474983

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: