Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,526,050 A→G T95A (ACC→GCC)  cysK → cysteine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,526,0500AG100.0% 225.1 / NA 67T95A (ACC→GCC) cysKcysteine synthase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/33);  total (34/33)

GTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGG  >  NZ_CP009273/2525963‑2526132
                                                                                       |                                                                                  
gTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGcc                                                                                  <  1:976031/90‑1 (MQ=255)
   gAACCGACCAGCGGTACTACCGGGATTGCCCTGGCTTATGTGGCTGCCGCTCGCGGTTACAAACTCCCCCTGACCCTGCCAGAAGCCATg                                                                               <  1:1209763/90‑1 (MQ=255)
   gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  1:187281/90‑1 (MQ=255)
   gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  1:536628/90‑1 (MQ=255)
   gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  1:484097/90‑1 (MQ=255)
    aaCCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGa                                                                              >  1:781819/1‑90 (MQ=255)
       cGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTa                                                                           >  2:336007/1‑90 (MQ=255)
       cGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTa                                                                           >  1:653326/1‑90 (MQ=255)
         aCCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTAtt                                                                         >  1:1073146/1‑90 (MQ=255)
                  aaTACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAcgccgc                                                                >  1:330147/1‑90 (MQ=255)
                    tACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCaa                                                              <  1:292282/90‑1 (MQ=255)
                       cGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgct                                                           >  1:1115278/1‑90 (MQ=255)
                        gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               >  1:800793/1‑69 (MQ=255)
                        gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               >  1:1500934/1‑69 (MQ=255)
                        gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  2:1500934/69‑1 (MQ=255)
                        gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  2:800793/69‑1 (MQ=255)
                        gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgctg                                                          >  2:997937/1‑90 (MQ=255)
                              gCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGaaa                                                    >  2:796371/1‑90 (MQ=255)
                              gCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGaaa                                                    >  1:1132061/1‑90 (MQ=255)
                                   ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGtt                                               >  1:1430739/1‑90 (MQ=255)
                                   ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGtt                                               <  2:295863/90‑1 (MQ=255)
                                   ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGtt                                               <  2:330147/90‑1 (MQ=255)
                                   ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGtt                                               >  1:1082801/1‑90 (MQ=255)
                                    gCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              <  1:421973/90‑1 (MQ=255)
                                    gCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              <  2:1215641/90‑1 (MQ=255)
                                        aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAc                                                                     <  2:1217532/63‑1 (MQ=255)
                                        aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAc                                                                     >  1:1217532/1‑63 (MQ=255)
                                        aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg                                          >  1:585696/1‑90 (MQ=255)
                                         tGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGc                                         <  2:275228/90‑1 (MQ=255)
                                          gTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCa                                        <  1:389855/90‑1 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               >  1:47717/1‑48 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  2:47717/48‑1 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  1:330924/1‑90 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  2:1322550/1‑90 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGGGCAAAc                                     >  1:949025/1‑90 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGCATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  1:448033/1‑90 (MQ=255)
                                              cTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAcc                                    >  2:1021179/1‑90 (MQ=255)
                                                gCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTg                                  >  1:1031122/1‑90 (MQ=255)
                                                  cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              >  2:133066/1‑76 (MQ=255)
                                                  cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              <  1:133066/76‑1 (MQ=255)
                                                  cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGt                                >  1:569076/1‑90 (MQ=255)
                                                    cTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTTGGTGCAAACCTGGTGc                              <  1:1417869/90‑1 (MQ=255)
                                                    cTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGc                              <  1:62347/90‑1 (MQ=255)
                                                    cTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGc                              <  1:49595/90‑1 (MQ=255)
                                                     tCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCt                             >  2:940718/1‑90 (MQ=255)
                                                            tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa                      <  1:966480/90‑1 (MQ=255)
                                                            tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa                      <  2:177842/90‑1 (MQ=255)
                                                            tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAACCTGGTGCTGACGGAAg                     <  1:924061/90‑1 (MQ=255)
                                                                    cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa              >  2:214313/1‑90 (MQ=255)
                                                                    cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa              >  2:1573147/1‑90 (MQ=255)
                                                                    cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa              <  1:790567/90‑1 (MQ=255)
                                                                    cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa              <  1:760341/90‑1 (MQ=255)
                                                                    cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa              <  2:585696/90‑1 (MQ=255)
                                                                      cccTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              >  1:385453/1‑56 (MQ=255)
                                                                      cccTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              <  2:385453/56‑1 (MQ=255)
                                                                         tGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCa         <  2:1430739/90‑1 (MQ=255)
                                                                          gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgct                                                           >  2:1158468/1‑39 (MQ=39)
                                                                          gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgct                                                           <  1:1158468/39‑1 (MQ=39)
                                                                          gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAt        >  1:1227910/1‑90 (MQ=255)
                                                                          gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAt        >  1:944681/1‑90 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    >  1:1484356/1‑90 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    >  2:980956/1‑90 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    <  1:1005324/90‑1 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    <  2:1006090/90‑1 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    <  1:980956/90‑1 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    <  1:917768/90‑1 (MQ=255)
                                                                                gCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAgg  >  1:396778/1‑90 (MQ=255)
                                                                                       |                                                                                  
GTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGG  >  NZ_CP009273/2525963‑2526132

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