Predicted mutation
evidence seq id position mutation annotation gene description
RA CP028101 514,508 G→A L2237L (TTG→TTA)  C7Y58_02630 → autotransporter domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP028101514,5080GA92.3% 19.1 / ‑3.5 13coding (6711/7290 nt)C7Y58_02630autotransporter domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/3);  minor base . (1/0);  total (10/3)
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

ATTTGGAACTAATGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTG‑GACAATATCAGGAGATATCTTTGTAGGACGTAATAGAATGCATAGAAAATACTTAGTAGTGGATGAAATATTTAATGCAAAAGCAAGATACTATACTTATGGAATAGGAGTAAAGAATGAAATAGGTAAAGAATTCAGATTAAGTGAATCATTCACATTAAGACCATATGCAGCATTAAAGTTAGAATATGGAAGAGTATCAAAGATAAGAGAAAAATCAGGAGAAATCAAGTTAG  >  CP028101/514265‑514749
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
aTTTGGAACTAATGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCAATTTGATGATAACAATAGCTTAAATTG‑Ga                                                                                                                                                                                                                                            >  2:463695/1‑251 (MQ=18)
     gAACTAATGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGTGGTATTGGTCGTTTCACAGGTATTTTTCATAATAACTTTTCGTTCAAAGCTCTCTTAAAATCACCCTCCCATATGTTATAATTTAAGGTAGTTTTCTTAAAATCATTGCCATTTGATGATAACAATAGCTT‑AATTGTGACAata                                                                                                                                                                                                                                       >  1:402895/1‑251 (MQ=2)
     gAACTAATGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTAAATTG‑GACAata                                                                                                                                                                                                                                       >  1:171546/1‑251 (MQ=18)
     gAACTAATGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTAAATTG‑GACAata                                                                                                                                                                                                                                       >  1:298711/1‑251 (MQ=18)
     gAACTAATGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTAAATTG‑GACAata                                                                                                                                                                                                                                       >  2:385827/1‑251 (MQ=18)
        cTAATGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTAAATTG‑GACAATATCa                                                                                                                                                                                                                                    >  1:59356/1‑251 (MQ=14)
        cTAATGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTAAATTG‑GACAATATCa                                                                                                                                                                                                                                    >  2:229637/1‑251 (MQ=14)
           aTGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTAAATTG‑GACAATATCAGga                                                                                                                                                                                                                                 <  1:21235/251‑1 (MQ=11)
           aTGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTAAATTG‑GACAATATCAGga                                                                                                                                                                                                                                 <  1:533890/251‑1 (MQ=11)
                                                                                                                                                                                                                                       tAACAATAGCTTAAATTG‑GACAATATCAGGAGATATCTTTGTAGGACGTAATAGAATGCATAGAAAATACTTAGTAGTGGATGAAATATTTAATGCAAAAGCAAGATACTATACTTATGGAATAGGAGTAAAGAATGAAATAGGTAAAGAATTCAGATTAAGTGAATCATTCACATTAAGACCATATGCAGCATTAAAGTTAGAATATGGAAGAGTATCAAAGATAAGAGAAAAATCAGGAGAAATCAAGt     >  1:21794/1‑251 (MQ=14)
                                                                                                                                                                                                                                       tAACAATAGCTTAAATTG‑GACAATATCAGGAGATATCTTTGTAGGACGTAATAGAATGCATAGAAAATACTTAGTAGTGGATGAAATATTTAATGCAAAAGCAAGATACTATACTTATGGAATAGGAGTAAAGAATGAAATAGGTAAAGAATTCAGATTAAGTGAATCATTCACATTAAGACCATATGCAGCATTAAAGTTAGAATATGGAAGAGTATCAAAGATAAGAGAAAAATCAGGAGAAATCAAGt     >  2:599005/1‑251 (MQ=14)
                                                                                                                                                                                                                                         aCAATAGCTTAAATTG‑GACAATATCAGGAGATATCTTTGTAGGACGTAATAGAATGCATAGAAAATACTTAGTAGTGGATGAAATATTTAATGCAAAAGCAAGATACTATACTTATGGAATAGGAGTAAAGAATGAAATAGGTAAAGAATTCAGATTAAGTGAATCATTCACATTAAGACCATATGCAGCATTAAAGTTAGAATATGGAAGAGTATCAAAGATAAGAGAAAAATCAGGAGAAATCAAGTTa   <  1:91969/251‑1 (MQ=14)
                                                                                                                                                                                                                                          cAATAGCTTAAATTG‑GACAATATCAGGAGATATCTTTGTAGGACGTAATAGAATGCATAGAAAATACTTAGTAGTGGATGAAATATTTAATGCAAAAGCAAGATACTATACTTATGGAATAGGAGTAAAGAATGAAATAGGTAAAGAATTCAGATTAAGTGAATCATTCACATTAAGACCATATGCAGCATTAAAGTTAGAATATGGAAGAGTATCAAAGATAAGAGAAAAATCAGGAGAAATCAAGTTAg  >  1:501873/1‑251 (MQ=14)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
ATTTGGAACTAATGGAGAATACAAAACAGATACAGCAGGAGTTATAGACTACAAGAATCATGCTTATGGAGTAGCTTATGTACATGAAAATGAAGATATTAAGCTAGGAAGAGGTATTGGATGGTACACAGGTATAGTTCATAATACATTTAAGTTCAAAGATATAGGAAAATCAAAAGAACATATGTTACAAGCAAAGGTAGGATTATTAAAATCAGTGCCATTTGATGATAACAATAGCTTGAATTG‑GACAATATCAGGAGATATCTTTGTAGGACGTAATAGAATGCATAGAAAATACTTAGTAGTGGATGAAATATTTAATGCAAAAGCAAGATACTATACTTATGGAATAGGAGTAAAGAATGAAATAGGTAAAGAATTCAGATTAAGTGAATCATTCACATTAAGACCATATGCAGCATTAAAGTTAGAATATGGAAGAGTATCAAAGATAAGAGAAAAATCAGGAGAAATCAAGTTAG  >  CP028101/514265‑514749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑