Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,976,083 T→G A313A (GCA→GCC 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,0830TG97.4% 137.9 / ‑5.9 38A313A (GCA→GCCaraHarabinose ABC transporter permease AraH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (10/27);  minor base C (1/0);  total (11/27)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTT  >  NZ_CP009273/1975950‑1976214
                                                                                                                                     |                                                                                                                                   
gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATTGTGGGGAAAAAAGCCAAATTTTTGCAGAAAAAAACATCAGACAGTGCGTTTCGCCTTTTGCTTTTTACGGGTGGAGATCACCCCCGCcag                                                                                                                                >  1:593993/1‑139 (MQ=255)
gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCcag                                                                                                                                <  1:180476/139‑1 (MQ=255)
gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCcag                                                                                                                                <  1:202706/139‑1 (MQ=255)
              cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTTGAAGATCACCGCGGCCAGCAGGATTAAGCcgc                                                                                                                  <  1:665496/139‑1 (MQ=255)
              cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCcgc                                                                                                                  <  1:757393/139‑1 (MQ=255)
              cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCcgc                                                                                                                  <  2:443005/139‑1 (MQ=255)
              cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCcgc                                                                                                                  <  2:195700/139‑1 (MQ=255)
              cGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCcgc                                                                                                                  <  1:266210/139‑1 (MQ=255)
                aagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGa                                                                                                                <  2:29380/139‑1 (MQ=255)
                   aaCAGCCAACCGTTAGGGAGGTTGGCTATGGTGGGAAGAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAAcc                                                                                                             <  1:600907/139‑1 (MQ=255)
                                  tGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa                                                                                               <  1:618585/138‑1 (MQ=255)
                                  tGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa                                                                                               <  1:547087/138‑1 (MQ=255)
                                      ggTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgag                                                                                          <  2:94140/139‑1 (MQ=255)
                                           gCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAata                                                                                     <  1:376181/139‑1 (MQ=255)
                                                  tggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCa                                                                              <  2:314575/139‑1 (MQ=255)
                                                      aaaaaaCGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGtt                                                                          >  2:723049/1‑139 (MQ=255)
                                                         aaaCGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGa                                                                                                                           <  1:22602/87‑1 (MQ=255)
                                                         aaaCGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGa                                                                                                                           >  2:22602/1‑87 (MQ=255)
                                                          aaCGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATg                                                                      <  1:342424/139‑1 (MQ=255)
                                                           aCGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgag                                                                                          >  1:161789/1‑118 (MQ=255)
                                                           aCGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgag                                                                                          <  2:161789/118‑1 (MQ=255)
                                                                      ttgCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACg                                                          >  2:663511/1‑139 (MQ=255)
                                                                           gAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGcc                                                     <  2:75819/139‑1 (MQ=255)
                                                                            aaaaaaGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCt                                                    <  1:98554/139‑1 (MQ=255)
                                                                                       aGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGa                                          <  2:114271/138‑1 (MQ=255)
                                                                                           aGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgag                                                                                          >  1:378970/1‑86 (MQ=255)
                                                                                           aGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgag                                                                                          <  2:378970/86‑1 (MQ=255)
                                                                                             tGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAGGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGc                                   >  2:357122/1‑139 (MQ=255)
                                                                                             tGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCcgc                                                                                                                  >  1:599307/1‑60 (MQ=255)
                                                                                             tGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCcgc                                                                                                                  <  2:599307/60‑1 (MQ=255)
                                                                                             tGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa                                                                                               <  2:411941/79‑1 (MQ=255)
                                                                                             tGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa                                                                                               >  1:411941/1‑79 (MQ=255)
                                                                                             tGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGc                                   >  1:534143/1‑139 (MQ=255)
                                                                                             tGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGc                                   <  1:271377/139‑1 (MQ=255)
                                                                                                 tttCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGccacc                               <  2:593993/139‑1 (MQ=255)
                                                                                                                       tCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGa                                                                                                                <  2:186607/36‑1 (MQ=255)
                                                                                                                       tCGAAGATCACCGCGGCCAGCAGGATTAAGCCGCGa                                                                                                                >  1:186607/1‑36 (MQ=255)
                                                                                                                              tCACCGCGGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCttt  <  1:663511/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
GTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTT  >  NZ_CP009273/1975950‑1976214

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