Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,786,991 T→C D18G (GAT→GGT)  yibB ← protein YibB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,786,9910TC89.5% 45.1 / ‑0.9 19D18G (GAT→GGT) yibBprotein YibB
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (8/9);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

AAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA  >  NZ_CP009273/3786906‑3787073
                                                                                     |                                                                                  
aaGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                                <  1:462304/90‑1 (MQ=255)
             cGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCa                                                                   <  1:265580/90‑1 (MQ=255)
                      tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  2:169554/90‑1 (MQ=255)
                       aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAg                                                         >  1:600514/1‑90 (MQ=255)
                       aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAg                                                         >  1:625053/1‑90 (MQ=255)
                         cTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       >  1:529352/1‑90 (MQ=255)
                                         cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  1:750406/90‑1 (MQ=255)
                                         cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  2:131225/90‑1 (MQ=255)
                                          tGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatga                                      >  1:605528/1‑90 (MQ=255)
                                          tGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatga                                      >  1:619214/1‑90 (MQ=255)
                                             aCGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGAttt                                   >  1:609474/1‑90 (MQ=255)
                                                aGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAt                                <  2:140408/90‑1 (MQ=255)
                                                aGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAATATGCAGTGATAATTGTGGTGGATGATTTCAtt                               >  2:472068/1‑90 (MQ=255)
                                                  cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  2:697477/90‑1 (MQ=255)
                                                  cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  2:619214/90‑1 (MQ=255)
                                                        tttCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                        <  1:680129/90‑1 (MQ=255)
                                                           cACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGa                     <  2:289926/90‑1 (MQ=255)
                                                                      ttattaGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAAtata          <  2:83863/90‑1 (MQ=255)
                                                                              aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGa  >  1:150611/1‑90 (MQ=255)
                                                                                     |                                                                                  
AAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA  >  NZ_CP009273/3786906‑3787073

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