Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,673,542 A→G K287E (AAG→GAG)  ydgH → DUF1471 family protein YdgH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,673,5420AG98.7% 241.8 / ‑5.8 78K287E (AAG→GAG) ydgHDUF1471 family protein YdgH
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (25/52);  total (25/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

ATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATA  >  NZ_CP009273/1673455‑1673627
                                                                                       |                                                                                     
atgatgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCggag                                                                                     <  2:316638/90‑1 (MQ=255)
  gatgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCggagaa                                                                                   <  2:2981346/90‑1 (MQ=255)
   atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                                  <  2:882429/90‑1 (MQ=255)
   atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                                  <  1:1483545/90‑1 (MQ=255)
   atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                                  <  1:2986481/90‑1 (MQ=255)
     ggTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAgg                                                                                <  1:379690/90‑1 (MQ=255)
     ggTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAgg                                                                                >  2:1042178/1‑90 (MQ=255)
     ggTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAgg                                                                                <  1:2344894/90‑1 (MQ=255)
       tCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTg                                                                              >  1:2464787/1‑90 (MQ=255)
        cccGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGc                                                                             <  2:1817196/90‑1 (MQ=255)
          cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                           <  1:2453640/90‑1 (MQ=255)
          cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                           <  2:1005177/90‑1 (MQ=255)
          cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                           <  2:1609313/90‑1 (MQ=255)
          cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                           >  2:1905886/1‑90 (MQ=255)
          cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                           <  2:2256885/90‑1 (MQ=255)
          cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                           <  2:489980/90‑1 (MQ=255)
           gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          <  2:184754/90‑1 (MQ=255)
           gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          <  2:2321365/90‑1 (MQ=255)
           gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          <  2:1319752/90‑1 (MQ=255)
           gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          <  1:2745164/90‑1 (MQ=255)
           gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          <  2:36580/90‑1 (MQ=255)
           gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          <  1:598835/90‑1 (MQ=255)
           gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          >  2:958270/1‑90 (MQ=255)
             tCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGt                                                                        <  1:1042178/90‑1 (MQ=255)
                aCAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTAtt                                                                     <  2:1169554/90‑1 (MQ=255)
                aCAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTAtt                                                                     <  2:2571693/90‑1 (MQ=255)
                aCAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTAtt                                                                     <  1:378154/90‑1 (MQ=255)
                aCAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTAtt                                                                     >  1:138926/1‑90 (MQ=255)
                aCAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTAtt                                                                     >  1:2994699/1‑90 (MQ=255)
                   gCATCAAATTCTCTGGCAGCTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTAcc                                                                  <  1:1663772/90‑1 (MQ=255)
                   gCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTAcc                                                                  >  1:697689/1‑90 (MQ=255)
                   gCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTAcc                                                                  >  1:2616184/1‑90 (MQ=255)
                   gCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTAcc                                                                  >  1:2802019/1‑90 (MQ=255)
                             ctctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          <  1:578405/72‑1 (MQ=255)
                             ctctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          >  2:578405/1‑72 (MQ=255)
                             ctctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCg                                                        >  2:1809690/1‑90 (MQ=255)
                              tctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGc                                                       >  1:596481/1‑90 (MQ=255)
                              tctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGc                                                       >  2:361774/1‑90 (MQ=255)
                                      cTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCa                                               <  2:1500341/90‑1 (MQ=255)
                                         cGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                            <  1:492271/58‑1 (MQ=255)
                                         cGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                            >  2:492271/1‑58 (MQ=255)
                                             aaCTTGACCAAGGTCCCCTATAAGGTTGCAAACCGCGCCGGGGAAAAAGGTGCTAAGTTTTACCACACCCCCCGCCAGTGGCAGGAACGt                                        <  2:2450726/90‑1 (MQ=255)
                                             aaCATGCCCGAAGTCTCCTATCGGGTTGCAAACCGGGCGGGGGAGAAAGGTGAAAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGt                                        <  1:1640051/90‑1 (MQ=255)
                                             aaCATGCCCGAAGTCCCCTATCAGGTTGCAAAACGTGCCGCGGAGAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGt                                        <  2:2221306/90‑1 (MQ=255)
                                             aaCATGACCGAGGTCTCCTATCGGGTTGCGAAACGGGCCGGGGAAAAAGGTGCTAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGt                                        <  1:916799/90‑1 (MQ=255)
                                             aaCATGACCGAAGTCTCCTATCGGGTTGCGAAACGGGCCGGGGAAAAAGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGt                                        <  1:766143/90‑1 (MQ=255)
                                             aaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGGGGAAAAAGGTGCTAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGt                                        <  1:984346/90‑1 (MQ=255)
                                             aaCATGACCGAAGTCTCCTATCAGGTTGCAAAACGTGCCGGGGAAAAAGGTGCAAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGt                                        <  1:284306/90‑1 (MQ=255)
                                                aTGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                            >  1:2587301/1‑51 (MQ=255)
                                                aTGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                            <  2:2587301/51‑1 (MQ=255)
                                                 tGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGGGGAGAAGGGTGCTAATTATTACCACATCCCCCGCCAGTGGCAGGAACGTGGta                                    <  2:2578620/90‑1 (MQ=255)
                                                   aCCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGgaat                                  >  2:1428684/1‑90 (MQ=255)
                                                     cGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGtaataa                                <  2:68588/90‑1 (MQ=255)
                                                       aaGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAAcc                              <  2:596481/90‑1 (MQ=255)
                                                        aGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          >  2:335736/1‑45 (MQ=255)
                                                        aGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          <  1:335736/45‑1 (MQ=255)
                                                           ctcCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGAc                          >  1:105626/1‑90 (MQ=255)
                                                            tcCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGAcc                         >  1:1900318/1‑90 (MQ=255)
                                                             ccTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGGGCTAAGTATTACCACATCACCCGCCAGTGGCAGGCACGTGGTAATAAACTTACCg                        >  2:2112691/1‑90 (MQ=255)
                                                               tATCGGGTTGCGAACCGTGCCGGGGAAAAAGGTGCTAAGTATTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc                      <  2:2023515/90‑1 (MQ=255)
                                                               tATCGGGTTGCAAAACGTGCGGAAAAAAAAGGTGCTAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc                      <  1:2441746/90‑1 (MQ=255)
                                                               tATCAGGTTGCGAAAGGTGCGGGGGAGAAAGGTGCAAAGTTTTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc                      <  1:3006276/90‑1 (MQ=255)
                                                               tATCAGGTTGCGAAACGTGCCGGGGAGAAAGGTGCTAAGTTTTACCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc                      <  1:1631151/90‑1 (MQ=255)
                                                               tATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc                      <  2:105626/90‑1 (MQ=255)
                                                                    ggTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAAACTGACCGTCAgcgc                 >  1:2475798/1‑90 (MQ=255)
                                                                    ggTTGCGAAACGTGCCGCGGAAAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                 <  2:380991/90‑1 (MQ=255)
                                                                    ggTTGCAAACCGTGCGGGGGAAAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                 <  2:1136479/90‑1 (MQ=255)
                                                                    ggTTGCAAACCGGGCCGCGGAAAAGGGTGTTAAGTTTTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                 <  1:1905886/90‑1 (MQ=255)
                                                                    ggTTGCAAAACGGGCCGCGGAGAAAGGTGCTAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                 <  1:272425/90‑1 (MQ=255)
                                                                     gTTGCGAAACGTGCTGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCa                >  1:1136479/1‑90 (MQ=255)
                                                                     gTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGGGGTAATAACCTGACCGTCAGCGCa                >  2:1965197/1‑90 (MQ=255)
                                                                       tGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGa              >  2:2835553/1‑90 (MQ=255)
                                                                       tGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAGCGCAGa              >  1:2119443/1‑90 (MQ=255)
                                                                         cGAAAGGTGCCGGGGAAAAAGGTGCAAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc            <  2:2119443/90‑1 (MQ=255)
                                                                         cGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc            <  2:804515/90‑1 (MQ=255)
                                                                         cGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc            <  2:298349/90‑1 (MQ=255)
                                                                         cGAAACGTGCCGCGGAGAAAGGTGCAAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc            <  2:2066872/90‑1 (MQ=255)
                                                                          gAAACGTGCCGCGGAAAAAGGTGCTAATTTTTCCCCCATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCt           <  2:884497/90‑1 (MQ=255)
                                                                                tGCCGGGGAAAAAGGTGCAAATTATTACCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataa     <  2:2994699/90‑1 (MQ=255)
                                                                                   cgcgGAGAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaata  <  2:379470/90‑1 (MQ=255)
                                                                                       |                                                                                     
ATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATA  >  NZ_CP009273/1673455‑1673627

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