Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,976,582 T→G Y147S (TAC→TCC)  araH ← arabinose ABC transporter permease AraH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,976,5820TG100.0% 153.4 / NA 40Y147S (TAC→TCC) araHarabinose ABC transporter permease AraH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (20/20);  total (20/20)

AACCAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGTACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGCCACGCCAATCCACA  >  NZ_CP009273/1976457‑1976717
                                                                                                                             |                                                                                                                                       
aaCCAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGATGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTc                                                                                                                             >  2:379610/1‑138 (MQ=255)
  ccAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGaa                                                                                                                          <  2:191482/139‑1 (MQ=255)
           aTCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCa                                                                                                                 <  1:282968/139‑1 (MQ=255)
                    gcgACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCa                                                                                                        <  1:194831/139‑1 (MQ=255)
                      gACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAAt                                                                                                      >  1:510345/1‑139 (MQ=255)
                          gTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCg                                                                                                  >  2:73553/1‑139 (MQ=255)
                              gCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGAt                                                                                              >  1:106344/1‑139 (MQ=255)
                              gCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGAt                                                                                              >  2:322027/1‑139 (MQ=255)
                                        cgcAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCatt                                                                                     >  1:544917/1‑138 (MQ=255)
                                                 aCCGAACCAGATGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGt                                                                           <  1:379610/139‑1 (MQ=255)
                                                         aGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGAt                                                                   <  2:403044/139‑1 (MQ=255)
                                                           ttGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATaa                                                                 >  2:298376/1‑139 (MQ=255)
                                                               cGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGa                                                                                                                           >  2:583450/1‑77 (MQ=255)
                                                               cGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGa                                                                                                                           <  1:583450/77‑1 (MQ=255)
                                                                   aCCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCa                                                         <  2:8295/139‑1 (MQ=255)
                                                                         ggCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGAcc                                                   <  2:634702/139‑1 (MQ=255)
                                                                             aagaagCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGc                                               >  1:414414/1‑139 (MQ=255)
                                                                              agaagCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGa                                                                                                                           >  1:698680/1‑62 (MQ=255)
                                                                              agaagCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGa                                                                                                                           <  2:698680/62‑1 (MQ=255)
                                                                              agaagCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGc                                               >  2:175692/1‑138 (MQ=255)
                                                                              agaagCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGc                                               >  1:54507/1‑138 (MQ=255)
                                                                                       cATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGaa                                      <  2:544917/138‑1 (MQ=255)
                                                                                             cGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGccc                                 <  1:322027/137‑1 (MQ=255)
                                                                                             cGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGccc                                 <  1:297257/137‑1 (MQ=255)
                                                                                              gATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGc                              <  2:54507/139‑1 (MQ=255)
                                                                                                     accgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGAcc                                                   >  1:492062/1‑111 (MQ=255)
                                                                                                     accgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGAcc                                                   <  2:492062/111‑1 (MQ=255)
                                                                                                     accgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAAcc                       >  1:741777/1‑139 (MQ=255)
                                                                                                     accgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAAcc                       >  2:153879/1‑139 (MQ=255)
                                                                                                         cTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACa                                                                                                                        >  1:72369/1‑38 (MQ=255)
                                                                                                         cTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACa                                                                                                                        <  2:72369/38‑1 (MQ=255)
                                                                                                             aCCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGc                <  1:175692/138‑1 (MQ=255)
                                                                                                              ccGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGCCa              <  1:73553/139‑1 (MQ=255)
                                                                                                              ccGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGCCa              <  1:298376/139‑1 (MQ=255)
                                                                                                                 tCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGCCACGc           >  2:413763/1‑139 (MQ=255)
                                                                                                                 tCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGGGATCAGAGCATTTATTTTCAGTATGGCGATAACAAAGCCATTGACCAGGCCACAGATCACGCCCAGCCACAAACCCGCTGCCAAGc           >  2:156796/1‑139 (MQ=255)
                                                                                                                   tGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGAcc                                                   >  1:75309/1‑97 (MQ=255)
                                                                                                                   tGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGAcc                                                   <  2:75309/97‑1 (MQ=255)
                                                                                                                   tGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGCCACGcc          >  1:569219/1‑138 (MQ=255)
                                                                                                                          atgGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGCCACGCCAATCcaca  <  1:413763/139‑1 (MQ=255)
                                                                                                                             |                                                                                                                                       
AACCAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGTACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGCCACGCCAATCCACA  >  NZ_CP009273/1976457‑1976717

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