Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,662,649 T→G H58P (CAC→CCC)  mlc ← DNA‑binding transcriptional repressor Mlc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,662,6490TG100.0% 89.1 / NA 27H58P (CAC→CCC) mlcDNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (14/13);  total (14/13)

GTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGC  >  NZ_CP009273/1662564‑1662726
                                                                                     |                                                                             
gTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGc                                                                           >  1:116447/1‑90 (MQ=255)
gTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGc                                                                           >  2:225619/1‑90 (MQ=255)
gTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGc                                                                           >  1:1060462/1‑90 (MQ=255)
gTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGc                                                                           >  1:470765/1‑90 (MQ=255)
gTGCCAAGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCATGCACCAGGGGTGc                                                                           >  1:702059/1‑90 (MQ=255)
  gCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCtt                                                                         <  1:1146641/90‑1 (MQ=255)
                tCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACgg                                                           >  2:720938/1‑90 (MQ=255)
                tCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACgg                                                           <  1:720938/90‑1 (MQ=255)
                  aaCCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGAc                                                         <  1:764367/90‑1 (MQ=255)
                  aaCCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGAc                                                         <  2:982674/90‑1 (MQ=255)
                         aGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTa                                                  >  1:288063/1‑90 (MQ=255)
                          gCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAg                                                 <  2:288063/90‑1 (MQ=255)
                          gCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAg                                                 >  1:1239358/1‑90 (MQ=255)
                          gCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAg                                                 <  1:755282/90‑1 (MQ=255)
                             ccACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGa                                              >  2:258777/1‑90 (MQ=255)
                               aCCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATa                                            <  1:1119920/90‑1 (MQ=255)
                                                   ccccGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTgcgc                        <  2:500273/90‑1 (MQ=255)
                                                    cccGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGcc                       >  2:1008345/1‑90 (MQ=255)
                                                    cccGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGcc                       >  1:665551/1‑90 (MQ=255)
                                                          tcttTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACgg                 >  1:273713/1‑90 (MQ=255)
                                                          tcttTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACgg                 >  2:941394/1‑90 (MQ=255)
                                                             ttGATTTCCAGCTCTTGCCCCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGaa              <  2:779261/90‑1 (MQ=255)
                                                                  ttCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATc         <  2:474386/90‑1 (MQ=255)
                                                                         tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc  <  2:45569/90‑1 (MQ=255)
                                                                         tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc  <  2:470765/90‑1 (MQ=255)
                                                                          ctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATc         >  1:281737/1‑82 (MQ=255)
                                                                          ctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATc         <  2:281737/82‑1 (MQ=255)
                                                                                     |                                                                             
GTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGC  >  NZ_CP009273/1662564‑1662726

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