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,6490TG95.5% 65.7 / ‑3.6 22H58P (CAC→CCC) mlcDNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (10/11);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGAC  >  NZ_CP009273/1662564‑1662736
                                                                                     |                                                                                       
gTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGc                                                                                     >  1:167876/1‑90 (MQ=255)
gTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGc                                                                                     >  1:187343/1‑90 (MQ=255)
gTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGc                                                                                     >  2:519051/1‑90 (MQ=255)
    cAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCg                                                                                 >  1:241597/1‑90 (MQ=255)
                  aaCCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGAc                                                                   >  1:16207/1‑90 (MQ=255)
                      accaGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAAtt                                                               <  1:355454/90‑1 (MQ=255)
                            cccACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTg                                                         <  1:539786/90‑1 (MQ=255)
                             ccACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGa                                                        <  1:2469/90‑1 (MQ=255)
                                     ggACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCa                                                >  1:489476/1‑90 (MQ=255)
                                      gACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAg                                               <  2:271585/90‑1 (MQ=255)
                                      gACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAg                                               <  2:468086/90‑1 (MQ=255)
                                               ggTTCCCCGCTTCTTTGATTTCCAGCTCATGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGtt                                      >  2:522944/1‑90 (MQ=255)
                                                    cccGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGcc                                 >  1:347733/1‑90 (MQ=255)
                                                                  ttCCAGCTCTTGCACCAGGGGTGCTTTGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATc                   <  1:331738/90‑1 (MQ=255)
                                                                    ccAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGa                 <  2:204251/90‑1 (MQ=255)
                                                                    ccAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGa                 >  2:221984/1‑90 (MQ=255)
                                                                         tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc            <  2:489476/90‑1 (MQ=255)
                                                                         tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc            <  2:505643/90‑1 (MQ=255)
                                                                         tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc            >  1:214077/1‑90 (MQ=255)
                                                                              cACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGAc       <  1:283972/90‑1 (MQ=255)
                                                                              cACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGAc       <  2:342887/90‑1 (MQ=255)
                                                                                   ggTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGAc  <  1:97440/90‑1 (MQ=255)
                                                                                     |                                                                                       
GTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGAC  >  NZ_CP009273/1662564‑1662736

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