Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,328,974 G→C P384A (CCA→GCA)  yfaA ← DUF2138 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,328,9740GC90.5% 61.3 / 0.4 21P384A (CCA→GCA) yfaADUF2138 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base C (12/7);  total (14/7)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

GAGTCTGGCTTACCGGCAAAACACCTTCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGGCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGAGTCTTC  >  NZ_CP009273/2328896‑2329043
                                                                              |                                                                     
gAGTCTGGCTTACCGGCAAAACACCTTCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCt                                                            >  2:840727/1‑90 (MQ=255)
gAGTCTGGCTTACCGGCAAAACACCTTCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCt                                                            >  2:751619/1‑90 (MQ=255)
       gCTTACCGGCAAAACACCTTCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGc                                                     >  1:34786/1‑90 (MQ=255)
            ccGGCAAAACACCTTCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGc                                                >  2:101727/1‑90 (MQ=255)
                      acCTTCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTg                                      <  1:632787/90‑1 (MQ=255)
                       ccTTCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGa                                     >  2:637698/1‑90 (MQ=255)
                        cTTCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGAc                                    >  1:301434/1‑90 (MQ=255)
                        cTTCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGAc                                    >  2:309961/1‑90 (MQ=255)
                         ttCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGAcc                                   >  1:533226/1‑90 (MQ=255)
                                     gctttCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCgg                       <  1:639813/90‑1 (MQ=255)
                                                    aaTATTTTGCGTAAACAGTTTCCCTGGCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGa        >  1:519942/1‑90 (MQ=255)
                                                     atatTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGAg       <  1:210291/90‑1 (MQ=255)
                                                     atatTTTGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGAg       <  2:34786/90‑1 (MQ=255)
                                                          ttGCGTAAACCGTTTCCCTGCCAACTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGAGTCTTc  <  2:573535/90‑1 (MQ=255)
                                                          ttGCGTAAACAGTTTCCCTGGCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGAGTCTTc  >  1:733010/1‑90 (MQ=255)
                                                          ttGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGAcc                                   >  2:11535/1‑57 (MQ=255)
                                                          ttGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGAcc                                   <  1:11535/57‑1 (MQ=255)
                                                          ttGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGAGTCTTc  <  1:98565/90‑1 (MQ=255)
                                                          ttGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGAGTCTTc  >  1:557725/1‑90 (MQ=255)
                                                          ttGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGAGTCTTc  >  2:826206/1‑90 (MQ=255)
                                                          ttGCGTAAACAGTTTCCCTGCCAATTGCGCCTGTTCGGCAGGGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGAGGCTTc  >  1:338185/1‑90 (MQ=255)
                                                                              |                                                                     
GAGTCTGGCTTACCGGCAAAACACCTTCTGGCGCTTTGCTTTCGTGCGCACCAATATTTTGCGTAAACAGTTTCCCTGGCAATTGCGCCTGTTCGGCAGTGCCATCAAACTGACCGACAAACAGCGGGGTTTGCAATTTTGAGTCTTC  >  NZ_CP009273/2328896‑2329043

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