Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,183,394 G→C G37G (GGC→GGG thiF ← thiazole biosynthesis adenylyltransferase ThiF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,183,3940GC100.0% 91.5 / NA 27G37G (GGC→GGGthiFthiazole biosynthesis adenylyltransferase ThiF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (15/12);  total (15/12)

GGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATT  >  NZ_CP009273/4183260‑4183521
                                                                                                                                      |                                                                                                                               
gggCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCa                                                                                                                             <  2:556806/139‑1 (MQ=255)
   cgatcgatATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGAc                                                                                                                          >  1:140083/1‑139 (MQ=255)
     atcgatATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCg                                                                                                                        <  1:420003/139‑1 (MQ=255)
                           agGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGa                                                                                                  <  1:69433/139‑1 (MQ=255)
                             gATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAg                                                                                                 >  2:314867/1‑138 (MQ=255)
                                                          caTCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATAtc                                                                   >  2:194028/1‑139 (MQ=255)
                                                          caTCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATAtc                                                                   >  1:435/1‑139 (MQ=255)
                                                            tcgtcgTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATAtcgt                                                                 <  2:575649/139‑1 (MQ=255)
                                                             cgtcgtCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATAtcgtc                                                                >  1:509470/1‑139 (MQ=255)
                                                                 gtcATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAg                                                             <  2:435/138‑1 (MQ=255)
                                                                           aGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTtgctgc                                                                                    <  1:319774/105‑1 (MQ=255)
                                                                           aGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTtgctgc                                                                                    >  2:319774/1‑105 (MQ=255)
                                                                             tACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTggcgg                                                <  1:314867/139‑1 (MQ=255)
                                                                                cAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCa                                                                                                                             >  1:235360/1‑59 (MQ=255)
                                                                                cAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCa                                                                                                                             <  2:235360/59‑1 (MQ=255)
                                                                                   cGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATaa                                          >  1:548927/1‑139 (MQ=255)
                                                                                       ccGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCa                                      >  1:576751/1‑139 (MQ=255)
                                                                                       ccGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCa                                      >  2:308583/1‑139 (MQ=255)
                                                                                                       ccAGGTACAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGtcatt                      <  2:548927/139‑1 (MQ=255)
                                                                                                              cAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGcc               <  2:4247/139‑1 (MQ=255)
                                                                                                              cAGCGCAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGcc               <  1:160183/139‑1 (MQ=255)
                                                                                                                 cgcAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCtat                                            >  2:21061/1‑107 (MQ=255)
                                                                                                                 cgcAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCtat                                            <  1:21061/107‑1 (MQ=255)
                                                                                                                 cgcAGCAGGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAAc            >  1:22702/1‑139 (MQ=255)
                                                                                                                      cagGTGTACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCGGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTg       >  1:155405/1‑139 (MQ=255)
                                                                                                                           gtACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAAtt  >  2:222077/1‑139 (MQ=255)
                                                                                                                           gtACCCAGCCCCCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAAtt  >  1:502135/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                               
GGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATT  >  NZ_CP009273/4183260‑4183521

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