Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,651,047 A→G V223A (GTC→GCC)  ynfC ← YnfC family lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,651,0470AG95.0% 61.8 / ‑4.1 20V223A (GTC→GCC) ynfCYnfC family lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (11/8);  total (12/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 = 9.95e-01

TCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTA  >  NZ_CP009273/1650979‑1651131
                                                                    |                                                                                    
tcCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGGCGGCGGGTTTTACTCCTtcat                                                                 >  1:116751/1‑90 (MQ=255)
 ccTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtcatc                                                                >  2:397828/1‑90 (MQ=255)
         cGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGa                                                        >  2:480292/1‑90 (MQ=255)
         cGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGGCGGCGAGTTTTACTCCTTCATCAACAATGa                                                        >  2:180140/1‑90 (MQ=255)
                         ttAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGGCGGCGGGtttt                                                                           >  1:37705/1‑55 (MQ=255)
                         ttAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGGCGGCGGGtttt                                                                           <  2:37705/55‑1 (MQ=255)
                         ttAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGGCGGCGAGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCAc                                        <  1:180140/90‑1 (MQ=255)
                                     gatcgTATTTTTGATGGTGTAAACCCGTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTg                            >  2:404776/1‑90 (MQ=255)
                                       tcgTATTTTTGATGGTGTAAACCCGTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGAc                          >  1:207414/1‑90 (MQ=255)
                                                   tGGTGTAAACCCGTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTcac              <  1:480292/90‑1 (MQ=255)
                                                     gtgtAAACCCGTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGAc                          <  2:112927/76‑1 (MQ=255)
                                                     gtgtAAACCCGTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGAc                          >  1:112927/1‑76 (MQ=255)
                                                        tAAACCCGTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGctc         <  2:116751/90‑1 (MQ=255)
                                                        tAAACCCGTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGctc         <  1:442664/90‑1 (MQ=255)
                                                               gTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTcac              >  1:8014/1‑78 (MQ=255)
                                                               gTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTcac              <  1:75394/78‑1 (MQ=255)
                                                               gTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTcac              >  2:75394/1‑78 (MQ=255)
                                                               gTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTcac              <  2:8014/78‑1 (MQ=255)
                                                               gTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa  >  2:182064/1‑90 (MQ=255)
                                                               gTTCGGCGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa  >  1:68479/1‑90 (MQ=255)
                                                                    |                                                                                    
TCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTA  >  NZ_CP009273/1650979‑1651131

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