Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,168,050 A→G Q151R (CAG→CGG)  ycfJ → PF05433 family protein YcfJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,168,0500AG88.9% 45.3 / 0.0 18Q151R (CAG→CGG) ycfJPF05433 family protein YcfJ
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (9/7);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

CGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCC  >  NC_000913/1167806‑1168294
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
cGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCGGCAggg                                                                                                                                                                                                                                                <  1:166258/251‑1 (MQ=255)
                                                 aTCAGTTTGGTGGTGGGCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGGGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCAGGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGc                                                                                                                                                                                               <  1:149675/251‑1 (MQ=255)
                                                       ttGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATcggcag                                                                                                                                                                                                                                                  >  1:15171/1‑194 (MQ=255)
                                                       ttGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATcggcag                                                                                                                                                                                                                                                  <  2:15171/194‑1 (MQ=255)
                                                                  cgcgGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGAtttt                                                                                                                                                                              >  1:49786/1‑251 (MQ=255)
                                                                                  ccACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAAACAGATCCAGGGCTCTCTCCAGGAACGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGtata                                                                                                                                                               >  1:306043/1‑250 (MQ=255)
                                                                                                 cgcTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTg                                                                                                                                               <  1:77208/251‑1 (MQ=255)
                                                                                                       gtggtgGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACACAGAATAACAAGGCTGTACTCTGCAAttt                                                                                                                                         >  1:122072/1‑251 (MQ=255)
                                                                                                                                ccAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCGGCAATTTGGCCCCTCATTCGCTCATGCTGAgg                                                                                                                <  1:104721/251‑1 (MQ=255)
                                                                                                                                            ccAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCtttttt                                                                                                       <  1:15174/248‑1 (MQ=255)
                                                                                                                                                     cGATACTTACACGACTACGCAACAGCGTTGTAAAACGGCGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATACGGGTACGCCGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGCATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTAttt                                                                                           >  1:33330/1‑251 (MQ=255)
                                                                                                                                                                         aaCAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGGAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGa                                                                                                                                                                                                   <  1:140088/127‑1 (MQ=255)
                                                                                                                                                                         aaCAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGGAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGa                                                                                                                                                                                                   >  2:140088/1‑127 (MQ=255)
                                                                                                                                                                                                                    cGGTTATGATGTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCt                                                                                                            <  2:84657/171‑1 (MQ=255)
                                                                                                                                                                                                                    cGGTTACGATGTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCt                                                                                                            >  1:84657/1‑171 (MQ=255)
                                                                                                                                                                                                                     ggTTATGATGTGAACTATAAGATTGGCGATCGGCCGGGCAACATCCGCATGGACCGCGATCAGGGTCAGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACCAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGcc                           >  2:131583/1‑250 (MQ=255)
                                                                                                                                                                                                                           gatgTGACCTATAAGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTcg                     >  1:119007/1‑251 (MQ=255)
                                                                                                                                                                                                                            atgTGACCTATAAGATTGGCGATAGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTcgc                    >  1:198866/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aGATTGGCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCg        >  1:243767/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gCGATCGGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGAGCCCGCTAGAGAGCAACGGGCAACTGATTTGGAATAACAAAGTATAACAAGGCTGTACGCTGCAATTTGGCCCCTCATTCTCTCAGTCTGAGGGGCTTTTTTTGCGGCTTATTTCACCAGTTCGGGCCATAAACGCAAAGGCGTGCCGGCAATCTGCATCAGTTTTGCCAGCGATGCACCGTCGCGTGCGCTGAGCGACATcc  <  2:33330/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
CGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTGCAATTTGGCCCCTCATTCGCTCAGGCTGAGGGGCTTTTTTTGCGACTTATTTCACCAGTTCGGGCCATAAACGCAAAGTCGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCC  >  NC_000913/1167806‑1168294

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