Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,765,154 A→G N95D (AAT→GAT)  ydiL → YdiL family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,765,1540AG96.8% 99.3 / ‑3.3 31N95D (AAT→GAT) ydiLYdiL family protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (16/14);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

TGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTACAACAATAATGAACGGATGCTGACACAA  >  NZ_CP009273/1765032‑1765286
                                                                                                                          |                                                                                                                                    
tGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTgg                                                                                                                        >  2:417117/1‑137 (MQ=255)
       gCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTAt                                                                                                               <  2:839109/139‑1 (MQ=255)
               aaaaaCATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGt                                                                                                       <  1:737618/139‑1 (MQ=255)
                 aaaCATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAg                                                                                                     <  1:516543/139‑1 (MQ=255)
                    cATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAgctg                                                                                                  <  1:642682/139‑1 (MQ=255)
                                  catAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGca                                                                                    >  1:94430/1‑139 (MQ=255)
                                    tAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGcaca                                                                                  <  2:121085/139‑1 (MQ=255)
                                           aataaataACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTg                                                                           >  1:205116/1‑139 (MQ=255)
                                           aataaataACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTg                                                                           >  1:779478/1‑139 (MQ=255)
                                                                       aTGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGa                                                >  1:66039/1‑138 (MQ=255)
                                                                            ttttttCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGaa                                          <  1:461676/139‑1 (MQ=255)
                                                                             tttttCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAAt                                         >  1:316744/1‑139 (MQ=255)
                                                                                   cccGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCgt                                   <  2:316744/139‑1 (MQ=255)
                                                                                        cTTAACCGCATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTa                              <  2:705778/139‑1 (MQ=255)
                                                                                                cATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAg                                                                        <  2:316413/89‑1 (MQ=255)
                                                                                                cATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAg                                                                        >  1:599573/1‑89 (MQ=255)
                                                                                                cATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAg                                                                        <  2:599573/89‑1 (MQ=255)
                                                                                                cATTTCAACGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAg                                                                        >  1:316413/1‑89 (MQ=255)
                                                                                                      aaCGAGTCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTACAACAATAATGaa                 <  1:417117/138‑1 (MQ=255)
                                                                                                            tCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGt                                                                                                       <  2:853133/46‑1 (MQ=255)
                                                                                                            tCTACCCTGATGGCGATTTTATTGACTGGAAGATCTATCAATCGGt                                                                                                       >  1:853133/1‑46 (MQ=255)
                                                                                                                     aTGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGCATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTACAACAATAATGAACGGATGCTGACacaa  >  1:441209/1‑138 (MQ=255)
                                                                                                                     aTGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCca                                                    >  2:271197/1‑88 (MQ=255)
                                                                                                                     aTGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCca                                                    <  1:271197/88‑1 (MQ=255)
                                                                                                                     aTGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCATAGAAAGGAATTGTCGTTGTTACAACAATAATGAACGGATGCTGACacaa  >  2:427143/1‑138 (MQ=255)
                                                                                                                     aTGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTacaac                          <  2:246365/114‑1 (MQ=255)
                                                                                                                     aTGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTacaac                          >  1:246365/1‑114 (MQ=255)
                                                                                                                     aTGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTACAACAATAATGAACGGATGCTGACacaa  >  1:640251/1‑138 (MQ=255)
                                                                                                                     aTGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTACAACAATAATGAACGGATGCTGACacaa  >  2:530029/1‑138 (MQ=255)
                                                                                                                     aTGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTACAACAATAATGAACGGATGCTGACacaa  >  1:369654/1‑138 (MQ=255)
                                                                                                                     aTGGCGATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTACAACAATAATGAACGGATGCTGACacaa  <  2:794998/138‑1 (MQ=255)
                                                                                                                          |                                                                                                                                    
TGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAAGGAATTGTCGTTGTTACAACAATAATGAACGGATGCTGACACAA  >  NZ_CP009273/1765032‑1765286

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