Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,953,437 Δ1 bp 29.4% pseudogene (178/243 nt) yifN ← PemK toxin family pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,953,4370C.29.4% 26.9 / 13.9 17pseudogene (178/243 nt)yifNPemK toxin family pseudogene
Reads supporting (aligned to +/- strand):  ref base C (4/7);  new base . (2/3);  total (6/11)
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

AGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAAGAATTGCAGGTCATTGATAACATT  >  CP009273/3953309‑3953572
                                                                                                                                |                                                                                                                                       
aGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTaa                                                                                                                     <  2:294009/149‑1 (MQ=255)
            ttAGGGCCTCTATGGGCTACAAAGTCATTGTCGGGGATGGCTGTGTTATTGTCAATTTATGAGCGATATGGCATAGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAGCATTGATTAGAGTACCCTCTGt                                                                                                         >  1:61967/1‑149 (MQ=255)
                         gggCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATAAGCAATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTg                                                                                           >  1:173298/1‑149 (MQ=255)
                          ggCTACAAAGTCATTTTCGGGTATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGt                                                                                          <  1:114227/149‑1 (MQ=255)
                                 aaGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATCTGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATc                                                                                    <  2:102950/149‑1 (MQ=255)
                                      attattGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGGATTGATTGACTTGATGGCCGCCCTCTGTATGTCTGCAATCATCTTTGTGTCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTc                                                                              <  1:140921/149‑1 (MQ=255)
                                      aTTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGttt                                                                               >  2:13930/1‑149 (MQ=255)
                                                   gTTATTTTAGTGTCAATTTATGAGCCATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACttgtt                                                                  <  2:169848/149‑1 (MQ=255)
                                                             tGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTc                                                        <  1:13930/149‑1 (MQ=255)
                                                                             atatGGCATTGATTGACTTGATGACCGCCCTCTGTATTTATGCAACCAACAATATAACAAAAAATAATTAAGATACCCTCTGTAAGTCCATGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCtttagaaaaacaaa                                       <  2:124521/149‑12 (MQ=255)
                                                                                              ttGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATCAACTATGAAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAAGaa                       >  1:59171/1‑149 (MQ=255)
                                                                                                gatgaCCGCCCTCTGTATGTCTGCAACCAACTATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAGCGAATCTGCTGATCAAAGAATTg                    >  2:47247/1‑149 (MQ=255)
                                                                                                   gaCCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAAGAATTGCa                  <  2:210488/149‑1 (MQ=255)
                                                                                                    aCCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAAGAATTGCAg                 <  1:269812/149‑1 (MQ=255)
                                                                                                                aTGTCTGCAACCAACTATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAAGAATTGCAGGTCATTGATAACa    <  2:173298/149‑1 (MQ=255)
                                                                                                                 tGTCTGCAACCTACTCATGTGGCAAACATTGTTTAAGAGACCCTTTGTAAGTTCTTGCTGTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAAGAATTGCAGGTCATTGATAACa    <  2:42568/149‑1 (MQ=255)
                                                                                                                   tCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAAGAATTGCAGGTCATTGATAACAtt  >  1:208711/1‑149 (MQ=255)
                                                                                                                                |                                                                                                                                       
AGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAAGAATTGCAGGTCATTGATAACATT  >  CP009273/3953309‑3953572

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