Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,004,326 T→G V256G (GTG→GGG)  yigL → sugar/pyridoxal phosphate phosphatase YigL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,004,3260TG100.0% 137.6 / NA 39V256G (GTG→GGG) yigLsugar/pyridoxal phosphate phosphatase YigL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (12/27);  total (12/27)

AGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATA  >  NZ_CP009273/4004240‑4004410
                                                                                      |                                                                                    
aGGCTGCATTATGGGCAGTGCGCATCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGcc                                                                                   <  2:110614/90‑1 (MQ=255)
     gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCAtt                                                                              >  1:852934/1‑90 (MQ=255)
       aTTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTAt                                                                            <  2:695011/90‑1 (MQ=255)
         tATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCt                                                                          <  2:769762/90‑1 (MQ=255)
         tATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCt                                                                          <  1:576708/90‑1 (MQ=255)
            gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCg                                                                       <  1:502395/90‑1 (MQ=255)
             ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGt                                                                      <  1:665475/90‑1 (MQ=255)
              gCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTa                                                                     <  1:646437/90‑1 (MQ=255)
               cAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTaa                                                                    >  2:559895/1‑90 (MQ=255)
                aGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTaaa                                                                   >  1:848748/1‑90 (MQ=255)
                 gTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAAc                                                                  >  2:264603/1‑90 (MQ=255)
                 gTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAAc                                                                  <  1:65923/90‑1 (MQ=255)
                 gTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAAc                                                                  <  2:821260/90‑1 (MQ=255)
                  tGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAAct                                                                 <  1:563569/90‑1 (MQ=255)
                  tGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAAct                                                                 >  2:57835/1‑90 (MQ=255)
                   gcgcACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAActc                                                                >  2:214377/1‑90 (MQ=255)
                      cACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTAt                                                             >  2:310897/1‑90 (MQ=255)
                           gCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATc                                                        <  1:703711/90‑1 (MQ=255)
                           gCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATc                                                        >  1:455045/1‑90 (MQ=255)
                                tGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAAt                                                   <  1:316565/90‑1 (MQ=255)
                                 gAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                                                  <  1:253844/90‑1 (MQ=255)
                                 gAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                                                  <  2:40142/90‑1 (MQ=255)
                                 gAAAGACCTTCATCCCGAGCTGGAAGTAATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                                                  <  2:326448/90‑1 (MQ=255)
                                     gACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTc                                              <  1:383958/90‑1 (MQ=255)
                                     gACCTTAATCCCGAGCTGGTTGGGATTGGCACTAATGCCGACGACGCGGGGCCGCTTAATCTTCGTAGAGTCTCTTTATCGTAATGGTTc                                              >  2:876883/1‑90 (MQ=255)
                                         ttCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                                                  >  1:309899/1‑82 (MQ=255)
                                         ttCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATc                                                  <  2:309899/82‑1 (MQ=255)
                                                 gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCa                                  <  1:739440/90‑1 (MQ=255)
                                                       gAAGTGATTGGTACTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTc                            <  2:455045/90‑1 (MQ=255)
                                                                   aCTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACtt                <  2:920567/90‑1 (MQ=255)
                                                                    cTAATGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTc               <  2:289949/90‑1 (MQ=255)
                                                                      aaTGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGc             >  2:790224/1‑90 (MQ=255)
                                                                        tGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTa           <  1:457551/90‑1 (MQ=255)
                                                                        tGCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTa           <  1:358831/90‑1 (MQ=255)
                                                                         gCCGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTAc          <  2:154551/90‑1 (MQ=255)
                                                                           cgacgaCGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACaa        >  2:92741/1‑90 (MQ=255)
                                                                                acgCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGAt   <  2:401416/90‑1 (MQ=255)
                                                                                acgCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGAt   <  1:564556/90‑1 (MQ=255)
                                                                                 cgcgGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATAGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATa  <  1:518157/90‑1 (MQ=255)
                                                                                      |                                                                                    
AGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATA  >  NZ_CP009273/4004240‑4004410

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