Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,249,474 C→A 12.1% A180A (GCG→GCT)  aruH ← arginine‑‑pyruvate transaminase AruH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,249,4740CA12.1% 98.2 / 3.0 33A180A (GCG→GCT) aruHarginine‑‑pyruvate transaminase AruH
Reads supporting (aligned to +/- strand):  ref base C (22/7);  new base A (2/2);  total (24/9)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

CCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGAAAGCCATTCTCGGGCTTGACCGGTACCGGCACCACCTTGGCGCCACAG  >  NC_002947/4249344‑4249609
                                                                                                                                  |                                                                                                                                       
ccGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGcccc                                                                                                                                    >  1:369196/1‑136 (MQ=255)
            cGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTtgtg                                                                                                                        >  2:643683/1‑136 (MQ=255)
             ggggCTGACATGCCCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGAGCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGc                                                                                                                       <  2:621666/136‑1 (MQ=255)
              gggCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCg                                                                                                                      <  1:643683/136‑1 (MQ=255)
                            gCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGa                                                                                                                                 <  1:14214/111‑1 (MQ=255)
                            gCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGa                                                                                                                                 >  2:14214/1‑111 (MQ=255)
                            gCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGcc                                                                                                        >  2:613346/1‑136 (MQ=255)
                                         agcTCGCTGTAGACCTCGTCAGAGATCAGCCTCAGATCCTGTGCAACGCACTGTTCGGCCAGCGCCTCCCTTGTGTCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTgc                                                                                           <  2:608298/136‑1 (MQ=255)
                                           ctctctGTAGACCTCGTCAGAGATCAACCACAGATCATAGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCTGGCTGTTCAGAGCCAGCGCGCGGGTgcgc                                                                                         >  1:79990/1‑136 (MQ=255)
                                           cTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTgcgc                                                                                         >  1:96758/1‑136 (MQ=255)
                                           cTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGATCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTgcgc                                                                                         >  1:293931/1‑136 (MQ=255)
                                           cTCGCTGTAGACCTCGTCAGAGATAAGCCACAGATCATGGGCAAAGCACAGTTCGGTCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCACCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTgcgc                                                                                         >  1:474792/1‑136 (MQ=255)
                                              gCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTTGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTAGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCggg                                                                                      >  1:36348/1‑136 (MQ=255)
                                              gCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTAGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCggg                                                                                      >  1:7436/1‑136 (MQ=255)
                                                             agagATCAGCCACAGATCATGGGCAACGCACAGTTCTGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCg                                                                       >  1:167354/1‑136 (MQ=255)
                                                             agagATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGcc                                                                       >  1:410125/1‑135 (MQ=255)
                                                             agagATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCg                                                                       >  1:674244/1‑136 (MQ=255)
                                                             agagATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCg                                                                       >  1:673943/1‑136 (MQ=255)
                                                             agagATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCg                                                                       >  2:232094/1‑136 (MQ=255)
                                                             agagATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCg                                                                       >  1:316785/1‑136 (MQ=255)
                                                             agagATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCg                                                                       >  2:5665/1‑136 (MQ=255)
                                                             agagATCAGCCACAGATAATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCg                                                                       >  1:78669/1‑136 (MQ=255)
                                                               agaTCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGAc                                                                     >  2:381671/1‑136 (MQ=255)
                                                                       cacaGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGATCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCAcgcg                                                             >  2:657353/1‑136 (MQ=255)
                                                                                   gCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGaa                                                 >  2:242522/1‑136 (MQ=255)
                                                                                    cAACGCAGAGTGCGGCCAGCGCCGCCCCGGTGGCGCGCGGCTGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGaaa                                                <  1:635812/136‑1 (MQ=255)
                                                                                         cacaGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGAAAGCCAt                                           <  2:137004/136‑1 (MQ=255)
                                                                                         cacaGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGAAAGCCAt                                           >  1:347649/1‑136 (MQ=255)
                                                                                                   ccAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGAAAGCCATTCTCGGGCtt                                 <  2:172202/136‑1 (MQ=255)
                                                                                                    cAGCGCCTCCCAGGTGGCGCGCGGCAGGCTAGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGAAAGCCATTCTCGGGCTTg                                <  2:7436/136‑1 (MQ=255)
                                                                                                    cAGCGCCTCCCAGGTGGCGCGCGGCAGGCTAGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGAAAGCCATTCTCGGGCTTg                                <  2:36348/136‑1 (MQ=255)
                                                                                                         ccTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGAAAGCCATTCTCGGGCTTGACCgg                           >  2:144027/1‑136 (MQ=255)
                                                                                                                       cgcgGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGAAAGCCATTCTCGGGCTTGACCGGTACCGGCaccacca             >  2:263485/1‑135 (MQ=255)
                                                                                                                                  cGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGAAAGCCATTCTCGGGCTTGACCGGTACCGGCACCACCTTGGCGCCACAg  <  1:177979/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
CCGGCAGGCTGCCGGGGCTGACATGCTCGCCTTCATACAGCAGCTCGCTGTAGACCTCGTCAGAGATCAGCCACAGATCATGGGCAACGCACAGTTCGGCCAGCGCCTCCCAGGTGGCGCGCGGCAGGCTCGCCCCTGAAGGGTTGTGCGGGCTGTTCAGGGCCAGCGCGCGGGTGCGCGGGGTAATGCGCTCGGCGACGTCACGCGGGCACACGCGAAAGCCATTCTCGGGCTTGACCGGTACCGGCACCACCTTGGCGCCACAG  >  NC_002947/4249344‑4249609

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑