Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,737,658 G→C 19.2% intergenic (‑157/+45) hutU ← / ← PP_5034 urocanate hydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,737,6580GC19.2% 92.2 / 19.7 44intergenic (‑157/+45)hutU/PP_5034urocanate hydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/22);  new base C (5/3);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-01

GCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAAC  >  NC_002947/5737516‑5737796
                                                                                                                                              |                                                                                                                                          
gCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGTTCGAGCGCCGCGCGGGCGgcgc                                                                                                                                     >  2:148994/1‑150 (MQ=255)
 cGGTCGTCCGCGGGCGGTTATGTCGTAGTCAATATACAACCCCCCCTTGCCTTGTATGTACAAGCATGGACAACGAAAATAATCGGTCCTTACGGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                    <  1:186420/150‑1 (MQ=255)
   gtcgtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCg                                                                                                                                  <  2:1189128/150‑1 (MQ=255)
    tcgtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATCATCGGACCTTACCGTTCCTGCCGTTGGAGTTATGGTTTGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                 <  1:812901/150‑1 (MQ=255)
    tcgtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                 <  2:1014026/150‑1 (MQ=255)
    tcattggcgggcggGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                 <  2:670161/144‑1 (MQ=255)
                               aaTATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATc                                                                                                             >  2:551579/1‑143 (MQ=255)
                               aaTATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATc                                                                                                             <  1:551579/143‑1 (MQ=255)
                                 tataCAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGa                                                                                                    >  1:400810/1‑150 (MQ=255)
                                    aCAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAAcc                                                                                                 <  1:949827/150‑1 (MQ=255)
                                                        atgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCg                                                                             <  1:332137/150‑1 (MQ=255)
                                                          gtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGa                                                                            >  1:1052732/1‑149 (MQ=255)
                                                          gtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGa                                                                            <  2:1052732/149‑1 (MQ=255)
                                                          gtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGACGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                           >  2:155200/1‑150 (MQ=255)
                                                              aaGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCAAGGAGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                       >  1:209210/1‑150 (MQ=255)
                                                               cgCATGGACAACCAAAATAATCGGAGCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGgc                                                                      <  2:515278/149‑1 (MQ=255)
                                                                gCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGACGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGgca                                                                     <  1:155200/150‑1 (MQ=255)
                                                                    ggACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCAc                                                                 >  1:261964/1‑150 (MQ=255)
                                                                          ccAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGcccc                                                           >  1:1055771/1‑150 (MQ=255)
                                                                          ccAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGcccc                                                           <  1:459385/150‑1 (MQ=255)
                                                                                  atCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAAc                                                   >  2:346520/1‑150 (MQ=255)
                                                                                               gTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAAt                                      <  2:261964/150‑1 (MQ=255)
                                                                                                ttCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACaa                                                    >  1:781242/1‑135 (MQ=255)
                                                                                                ttCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACaa                                                    <  2:781242/135‑1 (MQ=255)
                                                                                                 tCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATtc                                    <  2:1130218/150‑1 (MQ=255)
                                                                                                    tGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATtcatc                                 <  2:400810/150‑1 (MQ=255)
                                                                                                       cTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgcgc                                                                                                                     >  2:57332/1‑63 (MQ=255)
                                                                                                       cTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgcgc                                                                                                                     <  1:57332/63‑1 (MQ=255)
                                                                                                           agaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCt                          <  2:209210/150‑1 (MQ=255)
                                                                                                               ttATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCa                      <  2:1183589/150‑1 (MQ=255)
                                                                                                                  tGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGcc                   <  2:1055771/150‑1 (MQ=255)
                                                                                                                   gTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGc                            <  2:1166599/140‑1 (MQ=255)
                                                                                                                   gTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGc                            >  1:1166599/1‑140 (MQ=255)
                                                                                                                           tgaGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGc                            >  2:171432/3‑132 (MQ=255)
                                                                                                                           tgaGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGc                            <  1:171432/130‑1 (MQ=255)
                                                                                                                              gATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTcc                       <  1:673309/134‑1 (MQ=255)
                                                                                                                              gATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTcc                       >  2:673309/1‑134 (MQ=255)
                                                                                                                               aTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaa      >  1:617755/1‑150 (MQ=255)
                                                                                                                                 cGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaaca    >  2:576452/1‑150 (MQ=255)
                                                                                                                                  gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAg                                          >  1:804646/1‑111 (MQ=255)
                                                                                                                                  gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAg                                          <  2:804646/111‑1 (MQ=255)
                                                                                                                                   aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaacaac  >  2:941554/1‑150 (MQ=255)
                                                                                                                                      gccgcGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCa                                           >  1:580627/1‑106 (MQ=255)
                                                                                                                                      gccgcGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCa                                           <  2:580627/106‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                          
GCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAAC  >  NC_002947/5737516‑5737796

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