Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,654,145 A→G I249T (ATC→ACC)  iscS ← cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,654,1450AG100.0% 116.1 / NA 32I249T (ATC→ACC) iscScysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/15);  total (17/15)

TCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACAGCGCACCGATACCTTTCGGGCCATAGATTTTGTGAC  >  NZ_CP009273/2654010‑2654280
                                                                                                                                       |                                                                                                                                       
tCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCt                                                                                                                                      >  1:431047/1‑139 (MQ=255)
        aCCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCt                                                                                                                                      >  2:533672/1‑131 (MQ=255)
        aCCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCt                                                                                                                                      <  1:533672/131‑1 (MQ=255)
         ccGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAg                                                                                                                             >  1:190431/1‑139 (MQ=255)
         ccGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAg                                                                                                                             >  2:520379/1‑139 (MQ=255)
            ttCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAgg                                                                                                                            >  1:68205/1‑137 (MQ=255)
            ttCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCa                                                                                                                          <  2:387298/139‑1 (MQ=255)
                    aaCTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGcc                                                                                                                   <  2:423301/138‑1 (MQ=255)
                            cGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCa                                                                                                           <  2:622755/138‑1 (MQ=255)
                            cGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCa                                                                                                           <  1:138524/138‑1 (MQ=255)
                               tatCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATAcc                                                                                                       >  2:37184/1‑139 (MQ=255)
                                                            cGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGAccgccgcc                                                                                     >  1:499303/1‑128 (MQ=255)
                                                            cGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGAccgccgcc                                                                                     <  2:499303/128‑1 (MQ=255)
                                                                   cgcgCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGAt                                                                   >  1:232961/1‑139 (MQ=255)
                                                                     cgcAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGc                                                                 >  1:638744/1‑139 (MQ=255)
                                                                                  tctcGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGtt                                                     >  2:482381/1‑138 (MQ=255)
                                                                                      ggTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTac                                                  <  2:431047/137‑1 (MQ=255)
                                                                                      ggTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTac                                                  <  1:162021/137‑1 (MQ=255)
                                                                                      ggTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTac                                                  <  2:232961/137‑1 (MQ=255)
                                                                                      ggTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTac                                                  <  1:449564/137‑1 (MQ=255)
                                                                                      ggTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACAGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTacga                                                >  2:387261/1‑139 (MQ=255)
                                                                                        tCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTacgacg                                              >  2:534846/1‑139 (MQ=255)
                                                                                        tCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTacgacg                                              >  2:602017/1‑139 (MQ=255)
                                                                                        tCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTac                                                  >  1:532529/1‑135 (MQ=255)
                                                                                        tCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTac                                                  <  2:532529/135‑1 (MQ=255)
                                                                                                  ttcttTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGtacatac                                       <  2:565038/136‑1 (MQ=255)
                                                                                                  ttcttTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACAgc                                    >  1:277480/1‑139 (MQ=255)
                                                                                                    cttTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACAgc                                    <  1:602017/137‑1 (MQ=255)
                                                                                                    cttTTGCGATGCGATAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACAg                                     <  1:387261/136‑1 (MQ=255)
                                                                                                         gcgatgcgatAGGCCTCGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACAGCGCACCg                              <  1:216648/138‑1 (MQ=255)
                                                                                                                         cGCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACAGCGCACCGATACCTTTCGGGCCa               >  2:73344/1‑137 (MQ=255)
                                                                                                                          gCCCATTCCGACGGTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACAGCGCACCGATACCTTTCGGGCCa               >  1:638181/1‑136 (MQ=255)
                                                                                                                                    acggTCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACAGCGCACCGATACCTTTCGGGCCATAGATTTTGTGAc  <  1:73344/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGATGCGATAGGCCTCGCCCATTCCGACGATCTGGTGAACAGGCAGAGTGCCGGAACGCATACCGCGCTCGTGACCGCCGCCGTGCATTTGCGCTTCGATGCGTACGCGCGGTTTACGACGTACATACAGCGCACCGATACCTTTCGGGCCATAGATTTTGTGAC  >  NZ_CP009273/2654010‑2654280

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