Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,318,468 A→C A609A (GCT→GCG ftsH ← ATP‑dependent zinc metalloprotease FtsH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,318,4680AC97.5% 265.4 / ‑1.6 81A609A (GCT→GCGftsHATP‑dependent zinc metalloprotease FtsH
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (44/35);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-02

CATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCAT  >  NZ_CP009273/3318384‑3318553
                                                                                    |                                                                                     
cATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGcc                                                                                  <  2:645628/90‑1 (MQ=255)
 aTGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCt                                                                                 <  1:510903/90‑1 (MQ=255)
   ggTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTgg                                                                               >  1:562505/1‑90 (MQ=255)
   ggTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTgg                                                                               >  1:1082265/1‑90 (MQ=255)
   ggTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTgg                                                                               >  2:815380/1‑90 (MQ=255)
   ggTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTgg                                                                               >  2:511360/1‑90 (MQ=255)
   ggTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTgg                                                                               >  2:49113/1‑90 (MQ=255)
   ggTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTgg                                                                               >  1:1346585/1‑90 (MQ=255)
      gtTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttc                                                                            <  1:464945/90‑1 (MQ=255)
        tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttctt                                                                          >  2:356313/1‑90 (MQ=255)
        tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttctt                                                                          >  1:950150/1‑90 (MQ=255)
        tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttctt                                                                          >  2:1158899/1‑90 (MQ=255)
        tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttctt                                                                          >  2:1246090/1‑90 (MQ=255)
        tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttctt                                                                          >  2:910797/1‑90 (MQ=255)
        tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttctt                                                                          >  1:1210366/1‑90 (MQ=255)
        tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttctt                                                                          >  1:1284159/1‑90 (MQ=255)
           ccGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTccc                                                                       <  2:1210366/90‑1 (MQ=255)
             gggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAg                                                                     >  1:590950/1‑90 (MQ=255)
              ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  1:1379421/1‑90 (MQ=255)
              ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  2:570467/1‑90 (MQ=255)
              ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  2:1009943/1‑90 (MQ=255)
              ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  2:1075689/1‑90 (MQ=255)
              ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  2:555424/1‑90 (MQ=255)
              ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  2:1512752/1‑90 (MQ=255)
              ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  1:1549158/1‑90 (MQ=255)
              ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  2:28461/1‑90 (MQ=255)
              ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  1:1306653/1‑90 (MQ=255)
                    gCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTg                                                              >  1:952453/1‑90 (MQ=255)
                     cGTCCGCGTTCCACCACCCGGAGGAGGAGCTTTTGGCCTCCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGc                                                            <  2:495874/90‑1 (MQ=255)
                      gTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGc                                                            <  2:1368193/90‑1 (MQ=255)
                      gTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGc                                                            <  1:1075689/90‑1 (MQ=255)
                      gTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGc                                                            <  2:911380/90‑1 (MQ=255)
                      gTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGc                                                            <  1:555424/90‑1 (MQ=255)
                      gTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGc                                                            <  1:1022145/90‑1 (MQ=255)
                      gTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGc                                                            <  1:1229363/90‑1 (MQ=255)
                        aCGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCgg                                                          >  2:1027752/1‑90 (MQ=255)
                         cgcgGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGa                                                         <  1:930928/90‑1 (MQ=255)
                             gTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTa                                                     >  1:772112/1‑90 (MQ=255)
                                    aaCCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATcgcg                                              <  1:230696/90‑1 (MQ=255)
                                    aaCCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATcgcg                                              <  2:796266/90‑1 (MQ=255)
                                     aCCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGa                                             <  1:663252/90‑1 (MQ=255)
                                       cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCt                                                                                 >  1:1175651/1‑52 (MQ=255)
                                       cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCt                                                                                 <  2:1175651/52‑1 (MQ=255)
                                       cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           <  1:483019/90‑1 (MQ=255)
                                       cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           >  2:929697/1‑90 (MQ=255)
                                       cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           >  2:1433081/1‑90 (MQ=255)
                                       cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           >  1:174328/1‑90 (MQ=255)
                                       cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           >  2:759732/1‑90 (MQ=255)
                                       cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           <  1:184671/90‑1 (MQ=255)
                                       cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           <  1:49113/90‑1 (MQ=255)
                                         gACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTg                                         >  1:1302083/1‑90 (MQ=255)
                                              ggagCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATc                                    <  2:1082265/90‑1 (MQ=255)
                                               gagCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCa                                   >  1:1160367/1‑90 (MQ=255)
                                               gagCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCa                                   >  1:868931/1‑90 (MQ=255)
                                                   cTTTGGGCTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGtc                               >  2:820615/1‑90 (MQ=255)
                                                   cTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGtc                               >  2:650889/1‑90 (MQ=255)
                                                   cTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGtc                               >  2:1569268/1‑90 (MQ=255)
                                                   cTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGtc                               >  1:1280004/1‑90 (MQ=255)
                                                   cTTTGGACTACCATTGTCGCCAGAATTTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGtca                              >  1:296990/1‑90 (MQ=255)
                                                     ttGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGtcat                             <  1:1158899/90‑1 (MQ=255)
                                                     ttGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGtcat                             <  1:572665/90‑1 (MQ=255)
                                                       ggACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGtcatca                           <  1:28461/90‑1 (MQ=255)
                                                          cTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATc                        <  2:1503883/90‑1 (MQ=255)
                                                           tACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           <  2:276997/70‑1 (MQ=255)
                                                           tACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           >  1:276997/1‑70 (MQ=255)
                                                           tACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCt                       >  1:1459471/1‑90 (MQ=255)
                                                           tACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCt                       >  1:295773/1‑90 (MQ=255)
                                                            aCCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTg                      <  2:772112/90‑1 (MQ=255)
                                                             ccATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGc                     <  2:534130/90‑1 (MQ=255)
                                                              cATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCg                    <  1:1238881/90‑1 (MQ=255)
                                                              cATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCg                    <  2:103225/90‑1 (MQ=255)
                                                                    cGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGt              >  1:1476684/1‑90 (MQ=255)
                                                                    cGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGt              >  1:1450610/1‑90 (MQ=255)
                                                                        aGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGAt          <  2:711716/90‑1 (MQ=255)
                                                                        aGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGAt          <  2:905251/90‑1 (MQ=255)
                                                                        aGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGAt          <  1:1554643/90‑1 (MQ=255)
                                                                         gAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATa         >  2:1345333/1‑90 (MQ=255)
                                                                          aaTTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAg        <  2:180379/90‑1 (MQ=255)
                                                                            ttGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGtc      >  2:1001597/1‑90 (MQ=255)
                                                                              gTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGtctc    <  2:1459471/90‑1 (MQ=255)
                                                                                tAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCat  <  1:815380/90‑1 (MQ=255)
                                                                                    |                                                                                     
CATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCAT  >  NZ_CP009273/3318384‑3318553

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