Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,038,395 T→C N279N (AAT→AAC yihF → YdgA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,038,3950TC93.8% 212.1 / 2.1 65N279N (AAT→AACyihFYdgA family protein
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (26/35);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

TTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCAAACGATG  >  NZ_CP009273/4038265‑4038528
                                                                                                                                  |                                                                                                                                     
ttATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACt                                                                                                                               <  1:883063/139‑1 (MQ=255)
ttATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACt                                                                                                                               <  2:365417/139‑1 (MQ=255)
  atcaatcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCg                                                                                                                             >  2:1244765/1‑139 (MQ=255)
   tcaatcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCgg                                                                                                                            >  2:582737/1‑139 (MQ=255)
   tcaatcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCgg                                                                                                                            <  1:264638/139‑1 (MQ=255)
   tcaatcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCg                                                                                                                             <  1:669189/138‑1 (MQ=255)
   tcaatcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCg                                                                                                                             >  2:669189/1‑138 (MQ=255)
     aatcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTa                                                                                                                          <  2:486958/139‑1 (MQ=255)
     aatcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTa                                                                                                                          <  2:871257/139‑1 (MQ=255)
           ggggAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAgtggtg                                                                                                                    <  1:1284074/139‑1 (MQ=255)
             ggAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAg                                                                                                                  >  1:1191871/1‑139 (MQ=255)
                     ttGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCatt                                                                                                          >  1:807474/1‑139 (MQ=255)
                       gCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAgt                                                                                                                        <  2:1241604/123‑1 (MQ=255)
                       gCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTgg                                                                                                               >  1:103791/1‑132 (MQ=255)
                       gCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTgg                                                                                                               <  2:103791/132‑1 (MQ=255)
                       gCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCCCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAgt                                                                                                                        >  1:1241604/1‑123 (MQ=255)
                           aaaTCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGc                                                                                                    <  1:1244765/139‑1 (MQ=255)
                                       tCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTAtt                                                                                          <  2:684765/137‑1 (MQ=255)
                                       tCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTAtt                                                                                          <  2:39824/137‑1 (MQ=255)
                                        cGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGa                                                                                        <  2:1138232/138‑1 (MQ=255)
                                        cGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGa                                                                                        <  1:555015/138‑1 (MQ=255)
                                        cGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGa                                                                                        <  1:377598/138‑1 (MQ=255)
                                        cGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGa                                                                                        <  1:143682/138‑1 (MQ=255)
                                        cGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGa                                                                                        <  1:660615/138‑1 (MQ=255)
                                            gCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCt                                                                                   >  1:1210010/1‑139 (MQ=255)
                                             cAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCt                                                                                   <  1:582737/138‑1 (MQ=255)
                                                       gCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTgcgc                                                                        >  2:1207483/1‑139 (MQ=255)
                                                            aTCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCAACTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAAtt                                                                   >  2:140160/1‑139 (MQ=255)
                                                                cgTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTat                                                                >  1:329103/1‑138 (MQ=255)
                                                                 gTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatca                                                              <  1:1080032/139‑1 (MQ=255)
                                                                 gTTAAAGACAAAGATTATATCAATGTCAATCTCAACTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatca                                                              <  1:140160/139‑1 (MQ=255)
                                                                  ttAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatcat                                                             <  2:807474/139‑1 (MQ=255)
                                                                    aaaGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatcatc                                                            <  1:552393/138‑1 (MQ=255)
                                                                    aaaGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatcatc                                                            <  1:401471/138‑1 (MQ=255)
                                                                    aaaGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatcatc                                                            <  2:1191871/138‑1 (MQ=255)
                                                                     aaGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCaa                                                                                                                                 >  1:235443/1‑68 (MQ=255)
                                                                     aaGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCaa                                                                                                                                 <  2:235443/68‑1 (MQ=255)
                                                                     aaGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatca                                                              <  1:814663/135‑1 (MQ=255)
                                                                     aaGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatca                                                              >  2:814663/1‑135 (MQ=255)
                                                                             gATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTAGAACa                                                  <  2:1141171/139‑1 (MQ=255)
                                                                                  tatCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGc                                              >  2:915854/1‑138 (MQ=255)
                                                                                        tGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGa                                                                                                                   >  2:95189/1‑63 (MQ=255)
                                                                                        tGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGa                                                                                                                   <  1:95189/63‑1 (MQ=255)
                                                                                        tGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTg                                               <  1:306752/131‑1 (MQ=255)
                                                                                        tGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTg                                               >  2:306752/1‑131 (MQ=255)
                                                                                            aaTCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGa                                                                                        >  2:596608/1‑86 (MQ=255)
                                                                                            aaTCTCACCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGa                                                                                        <  1:596608/86‑1 (MQ=255)
                                                                                                  aCCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAAt                                                                                               >  2:802384/1‑73 (MQ=255)
                                                                                                  aCCTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAAt                                                                                               <  1:802384/73‑1 (MQ=255)
                                                                                                   ccTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGa                                                                                        <  2:161573/79‑1 (MQ=255)
                                                                                                   ccTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGa                                                                                        >  1:161573/1‑79 (MQ=255)
                                                                                                    cTACGAACTTGATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTtgct                           >  1:411450/1‑139 (MQ=255)
                                                                                                              gATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatcatc                                                            <  2:1045044/96‑1 (MQ=255)
                                                                                                              gATAAGTTAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatcatc                                                            >  1:1045044/1‑96 (MQ=255)
                                                                                                                    ttAACAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAg           <  1:278715/139‑1 (MQ=255)
                                                                                                                      aaCAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatcatc                                                            <  1:806099/88‑1 (MQ=255)
                                                                                                                      aaCAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatcatc                                                            >  2:806099/1‑88 (MQ=255)
                                                                                                                      aaCAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGc          >  2:846929/1‑138 (MQ=255)
                                                                                                                       aCAAAAGGGAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCa         <  1:1207483/138‑1 (MQ=255)
                                                                                                                             gggAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCAAACGAt   >  2:20032/1‑138 (MQ=255)
                                                                                                                              ggAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatcatc                                                            <  1:602413/80‑1 (MQ=255)
                                                                                                                              ggAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatcatc                                                            >  2:602413/1‑80 (MQ=255)
                                                                                                                              ggAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGca                        >  2:473373/1‑116 (MQ=255)
                                                                                                                              ggAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGca                        <  1:473373/116‑1 (MQ=255)
                                                                                                                              ggAACCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCAAACGATg  >  2:290605/1‑138 (MQ=255)
                                                                                                                                  |                                                                                                                                     
TTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCAAACGATG  >  NZ_CP009273/4038265‑4038528

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