Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,786,139 T→C intergenic (‑228/+47) kbl ← / ← yibB glycine C‑acetyltransferase/protein YibB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,786,1390TC100.0% 47.4 / NA 14intergenic (‑228/+47)kbl/yibBglycine C‑acetyltransferase/protein YibB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/9);  total (5/9)

TGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTT  >  NZ_CP009273/3786057‑3786225
                                                                                  |                                                                                      
tGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCt                                                                                     >  1:108299/1‑86 (MQ=255)
tGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCt                                                                                     <  2:108299/86‑1 (MQ=255)
tGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTAtct                                                                                 <  2:452987/90‑1 (MQ=255)
                ggCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTAtc                                                                                  >  1:123301/1‑73 (MQ=255)
                ggCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTAtc                                                                                  <  2:123301/73‑1 (MQ=255)
                ggCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTATCTCAGTTGGTATCAAAAc                                                                 >  2:483820/1‑90 (MQ=255)
                  cATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTATCTCAGTTGGTATCAAAACGc                                                               <  1:173799/90‑1 (MQ=255)
                     aaTGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTATCTCAGTTGGTATCAAAACGCTAc                                                            <  2:107063/90‑1 (MQ=255)
                     aaTGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTATCTCAGTTGGTATCAAAACGCTAc                                                            <  2:167797/90‑1 (MQ=255)
                       tgatgaTTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTATCTCAGTTGGTATCAAAACGCTAccc                                                          <  2:61584/90‑1 (MQ=255)
                              tgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATg                                                   <  1:483820/90‑1 (MQ=255)
                                attTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGa                                                 >  1:59818/1‑90 (MQ=255)
                                     aCTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTCCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaa                                            <  1:221407/90‑1 (MQ=255)
                                                  aaaaaGAAAAAATGATGCCAAAGCATCATTTTCCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTg                               >  1:439021/1‑90 (MQ=255)
                                                                               tttcCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt  <  2:409440/90‑1 (MQ=255)
                                                                               tttcCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt  <  2:439021/90‑1 (MQ=255)
                                                                                  |                                                                                      
TGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTT  >  NZ_CP009273/3786057‑3786225

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