Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,982,555 T→C F232S (TTC→TCC)  yqeI → transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,982,5550TC94.7% 217.6 / 4.4 75F232S (TTC→TCC) yqeItranscriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (26/45);  total (28/47)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

ATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCC  >  NZ_CP009273/2982470‑2982641
                                                                                     |                                                                                      
aTTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCct                                                                                    >  2:94116/1‑90 (MQ=255)
aTTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCct                                                                                    <  2:1138752/90‑1 (MQ=255)
   gATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCtt                                                                                 >  2:311559/1‑90 (MQ=255)
   gATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCtt                                                                                 >  1:1018593/1‑90 (MQ=255)
     tGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCctct                                                                                  >  1:1179548/1‑87 (MQ=255)
     tGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCctct                                                                                  <  2:1179548/87‑1 (MQ=255)
     tGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGc                                                                               >  2:1049137/1‑90 (MQ=255)
        gTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCaaa                                                                            >  2:711412/1‑90 (MQ=255)
        gTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCaaa                                                                            >  1:1453039/1‑90 (MQ=255)
          tGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAAt                                                                          >  1:794089/1‑90 (MQ=255)
                agTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAc                                                                    >  2:528704/1‑90 (MQ=255)
                 gTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAcc                                                                   >  2:1169790/1‑90 (MQ=255)
                 gTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAcc                                                                   >  2:526313/1‑90 (MQ=255)
                 gTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAcc                                                                   >  1:1111973/1‑90 (MQ=255)
                 gTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCAcc                                                                   >  1:1335699/1‑90 (MQ=255)
                  ttCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCg                                                                  <  2:1335699/90‑1 (MQ=255)
                  ttCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCg                                                                  <  1:358954/90‑1 (MQ=255)
                      aTAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGga                                                              <  2:7097/90‑1 (MQ=255)
                        aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  1:773430/90‑1 (MQ=255)
                        aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  2:781810/90‑1 (MQ=255)
                        aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  2:1227508/90‑1 (MQ=255)
                        aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  1:140201/90‑1 (MQ=255)
                        aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  2:872902/90‑1 (MQ=255)
                        aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  2:902024/90‑1 (MQ=255)
                        aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  1:730342/90‑1 (MQ=255)
                        aaaTATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatg                                                            <  1:1169790/90‑1 (MQ=255)
                          atatAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGgatgat                                                          <  1:94116/90‑1 (MQ=255)
                          atatAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGgatgat                                                          <  2:909554/90‑1 (MQ=255)
                             taAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGt                                                       <  2:59868/90‑1 (MQ=255)
                              aaaaCACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGtt                                                      >  1:1294655/1‑90 (MQ=255)
                                    cGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCa                                                >  1:364316/1‑90 (MQ=255)
                                    cGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCa                                                >  1:528386/1‑90 (MQ=255)
                                    cGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCa                                                >  1:1449951/1‑90 (MQ=255)
                                    cGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCa                                                >  1:484295/1‑90 (MQ=255)
                                    cGGATACTGGATTCGGGGATTAACTGTAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCAt                                               >  1:219320/1‑90 (MQ=255)
                                      gATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATg                                              <  2:910873/90‑1 (MQ=255)
                                      gATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATg                                              <  1:229823/90‑1 (MQ=255)
                                      gATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATg                                              <  1:213248/90‑1 (MQ=255)
                                      gATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATg                                              <  1:213247/90‑1 (MQ=255)
                                        tACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCa                                            <  2:364316/90‑1 (MQ=255)
                                            ggATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGc                                             <  1:820843/85‑1 (MQ=255)
                                            ggATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGc                                             >  2:820843/1‑85 (MQ=255)
                                            ggATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATaa                                        >  2:91627/1‑90 (MQ=255)
                                            ggATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATaa                                        >  1:1205040/1‑90 (MQ=255)
                                             gATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATaaa                                       <  2:647099/90‑1 (MQ=255)
                                             gATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATaaa                                       <  1:201698/90‑1 (MQ=255)
                                             gATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATaaa                                       <  2:1294655/90‑1 (MQ=255)
                                             gATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATaaa                                       <  1:903720/90‑1 (MQ=255)
                                              aTTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATaaaa                                      >  2:53530/1‑90 (MQ=255)
                                                tCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAAc                                    >  2:900044/1‑90 (MQ=255)
                                                 cGGGGATTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACt                                   >  2:742893/1‑90 (MQ=255)
                                                      aTTAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAAc                              >  2:375088/1‑90 (MQ=255)
                                                       ttAACTGTAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAAcc                             <  1:472693/90‑1 (MQ=255)
                                                              tAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGa                      <  2:653484/90‑1 (MQ=255)
                                                              tAAAAAACATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGa                      <  2:1205040/90‑1 (MQ=255)
                                                                aaaaaCATCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAg                    <  1:528704/90‑1 (MQ=255)
                                                                      aTCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCaa              <  1:1098524/90‑1 (MQ=255)
                                                                      aTCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCaa              <  2:1453039/90‑1 (MQ=255)
                                                                       tCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGGATGCAATAAAAACTATAACCAACATGAAGTGGCaaa             >  2:1561482/1‑90 (MQ=255)
                                                                       tCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCaaa             >  2:448705/1‑90 (MQ=255)
                                                                       tCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCaaa             <  1:262238/90‑1 (MQ=255)
                                                                       tCCGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCaaa             <  1:1529217/90‑1 (MQ=255)
                                                                        ccGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGGATGCAATAAAAACTATAACCAACATGAAGTGGCAAAc            <  1:1561482/90‑1 (MQ=255)
                                                                        ccGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAAc            <  2:57686/90‑1 (MQ=255)
                                                                        ccGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAAc            <  2:19820/90‑1 (MQ=255)
                                                                        ccGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAAc            <  2:146647/90‑1 (MQ=255)
                                                                        ccGTGGTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAAc            <  1:711412/90‑1 (MQ=255)
                                                                             gTTATACTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGttt       <  2:52594/90‑1 (MQ=255)
                                                                                ataCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAc    <  2:1108520/90‑1 (MQ=255)
                                                                                  aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  1:171692/90‑1 (MQ=255)
                                                                                  aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  2:53869/90‑1 (MQ=255)
                                                                                  aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  1:433432/90‑1 (MQ=255)
                                                                                  aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  2:1297709/90‑1 (MQ=255)
                                                                                  aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  2:809008/90‑1 (MQ=255)
                                                                                  aCTCCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAAccc  <  2:1513467/90‑1 (MQ=255)
                                                                                     |                                                                                      
ATTGATGGGTTGAAAGAGTTCGATAAATATAAAACACGGATACTGGATTCGGGGATTAACTGTAAAAAACATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCC  >  NZ_CP009273/2982470‑2982641

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