Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,239,106 G→T 13.3% intergenic (‑75/‑300) argF ← / → PP_1080 ornithine carbamoyltransferase/molybdopterin‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,239,1060GT13.3% 69.1 / 3.0 30intergenic (‑75/‑300)argF/PP_1080ornithine carbamoyltransferase/molybdopterin‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (12/14);  new base T (2/2);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

GGATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTT  >  NC_002947/1238968‑1239249
                                                                                                                                          |                                                                                                                                               
ggATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCTTACGCCGCag                                                                                                                                       >  1:364001/1‑148 (MQ=255)
ggATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCTTACGCCGCag                                                                                                                                       <  2:364001/149‑2 (MQ=255)
      ccccGAGCAATTCGTCGGTGGTGAAGGCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCgg                                                                                                                                >  2:610307/1‑150 (MQ=255)
       cccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACgggctgta                                                                                                                               >  2:462596/1‑145 (MQ=255)
       cccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGa                                                                                                                               <  2:145654/150‑1 (MQ=255)
                ttCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAgg                                                                                                                      <  2:680353/150‑1 (MQ=255)
                         ggtgAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGgcgc                                                                                                             <  2:739114/150‑1 (MQ=255)
                             aaGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATa                                                                                                         <  2:175857/150‑1 (MQ=255)
                                  cagcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCg                                                                                                                                 >  1:327189/1‑121 (MQ=255)
                                  cagcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCg                                                                                                                                 <  2:327189/121‑1 (MQ=255)
                                  cagcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACtata                                                                                                    >  2:389141/1‑150 (MQ=255)
                                              aaaGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAAt                                                                                             <  1:35800/145‑1 (MQ=255)
                                              aaaGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAAt                                                                                             >  2:35800/1‑145 (MQ=255)
                                                               tGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTTACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGa                                                                       >  2:488517/1‑150 (MQ=255)
                                                                gAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTTACAACGGGGGTGACCTGCGGTGTACGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAggggggca                                                                      >  2:678192/1‑148 (MQ=255)
                                                                           ttATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGa                                                           >  2:102266/1‑150 (MQ=255)
                                                                            tATCTGCAACGATATGCAAATAGCCCGTATGATTTTTTTTAACAACTGAGGTTACCTTCGGCGTAAGCCGCACGGTGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAt                                                          <  2:271038/150‑1 (MQ=255)
                                                                             aTCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAGCGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAtt                                                         <  1:102266/150‑1 (MQ=255)
                                                                                       gTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtg                                               >  1:419881/1‑150 (MQ=255)
                                                                                         ttgtaGGTCGGCCGTATGGTTTTTTTTGACAACGGGGGCGTGCTGCGGCGTAAGCCGCAGGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgtt                                             <  1:384040/146‑1 (MQ=255)
                                                                                         aTGCCGGTCGTCCATATGTTTTGGTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGATGTGGCAGTAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgtt                                             <  1:610307/150‑1 (MQ=255)
                                                                                         aTGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgtt                                             <  1:486196/150‑1 (MQ=255)
                                                                                                  gACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgttgcaaagggt                                    <  1:389141/150‑1 (MQ=255)
                                                                                                         tGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTa                             >  1:597742/1‑150 (MQ=255)
                                                                                                                ttttGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATaa                                                                                                   >  2:44681/1‑73 (MQ=255)
                                                                                                                ttttGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATaa                                                                                                   <  1:44681/73‑1 (MQ=255)
                                                                                                                ttttGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGAAGAAATAGGGAAAGCCGCAATACTATAAGTAAATGGCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTg                      <  1:462596/150‑1 (MQ=255)
                                                                                                                               ccccgtGCGGCTTACGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATtgcaaagggtgttgcaaagggtgt                                  <  1:546480/118‑1 (MQ=255)
                                                                                                                               ccccgtGCGGCTTACGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATtgcaaagggtgttgcaaagggtgt                                  >  2:546480/6‑123 (MQ=255)
                                                                                                                                    tGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCGGGAATTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTAttt  >  1:228644/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                               
GGATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTT  >  NC_002947/1238968‑1239249

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