Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,034,572 T→G Q74H (CAA→CAC mobB ← molybdopterin‑guanine dinucleotide biosynthesis protein MobB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,034,5720TG100.0% 96.4 / NA 30Q74H (CAA→CACmobBmolybdopterin‑guanine dinucleotide biosynthesis protein MobB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/17);  total (13/17)

TCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATG  >  NZ_CP009273/4034486‑4034658
                                                                                      |                                                                                      
tCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgc                                                                                     <  2:217989/90‑1 (MQ=255)
    gCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctggc                                                                                 <  1:193879/90‑1 (MQ=255)
    gCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctggc                                                                                 >  1:413654/1‑90 (MQ=255)
      tttGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTAATCAAGGCCCATCGGTgctggctg                                                                               <  1:18424/90‑1 (MQ=255)
             tATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGa                                                                        <  1:57668/90‑1 (MQ=255)
             tATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGa                                                                        >  2:69834/1‑90 (MQ=255)
              aTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAt                                                                       >  2:504477/1‑90 (MQ=255)
              aTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAt                                                                       >  1:509233/1‑90 (MQ=255)
                   ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGttt                                                                  >  2:389534/1‑90 (MQ=255)
                                aaaCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGc                                                     <  1:504477/90‑1 (MQ=255)
                                 aaCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGcc                                                    >  1:441188/1‑90 (MQ=255)
                                       aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                              <  2:8399/90‑1 (MQ=255)
                                        gATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCa                                             <  2:312363/90‑1 (MQ=255)
                                            cAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTc                                         <  1:278232/90‑1 (MQ=255)
                                             aGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                              >  2:128678/1‑84 (MQ=255)
                                             aGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                              <  1:128678/84‑1 (MQ=255)
                                             aGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCa                                        >  2:401494/1‑90 (MQ=255)
                                                 cttcttCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                    <  2:86516/90‑1 (MQ=255)
                                                   tcttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGa                                                                        >  2:243036/1‑52 (MQ=255)
                                                   tcttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGa                                                                        <  1:243036/52‑1 (MQ=255)
                                                   tcttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTa                                  >  1:503022/1‑90 (MQ=255)
                                                      tcGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCt                               >  2:18383/1‑90 (MQ=255)
                                                       cGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCtt                              <  1:222141/90‑1 (MQ=255)
                                                                     gTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCa                <  2:518544/90‑1 (MQ=255)
                                                                          caAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc           <  2:509233/90‑1 (MQ=255)
                                                                            aGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCa         <  2:147814/90‑1 (MQ=255)
                                                                                 cATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                    <  2:143368/58‑1 (MQ=255)
                                                                                 cATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                    >  1:143368/1‑58 (MQ=255)
                                                                                 cATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCa    >  2:391157/1‑90 (MQ=255)
                                                                                   tCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATg  <  1:25377/90‑1 (MQ=255)
                                                                                      |                                                                                      
TCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATG  >  NZ_CP009273/4034486‑4034658

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