Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,737,599 G→C 8.5% intergenic (‑154/+48) 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,5990GC 8.5% 135.4 / 4.1 47intergenic (‑154/+48)hutU/PP_5034urocanate hydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base C (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

TCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACCC  >  NC_002947/5737457‑5737742
                                                                                                                                               |                                                                                                                                               
tctGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcgggcg                                                                                                                                           >  6:267353/1‑149 (MQ=255)
 ctGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACGCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcgggcg                                                                                                                                           >  1:85080/1‑149 (MQ=255)
  tGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcgggcg                                                                                                                                           <  1:76953/147‑1 (MQ=255)
          tcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAACTAATCGCACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                 <  5:173095/149‑1 (MQ=255)
          tcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGCGTATAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                 <  5:167856/149‑1 (MQ=255)
            cgcgGGCGGGTATGTCGTGGTCAATATACAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGTGCGCCGCGCGTGCGGCGCTCGAtctc                                                                                                                               <  8:96592/149‑1 (MQ=255)
                            gTGGTCAATATACAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATc                                                                                                               <  2:85080/149‑1 (MQ=255)
                                tCAATATACAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACa                                                                                                           <  5:267353/149‑1 (MQ=255)
                                    tataCAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGcc                                                                                                        >  4:344421/1‑148 (MQ=255)
                                     ataCAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAg                                                                                                          <  8:258420/145‑1 (MQ=255)
                                     ataCAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGcag                                                                                                       >  7:258420/1‑146 (MQ=255)
                                         aaaCACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAg                                                                                                          >  8:108644/1‑141 (MQ=255)
                                         aaaCACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAg                                                                                                          <  7:108644/141‑1 (MQ=255)
                                         aaaCACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAAc                                                                                                    >  6:230097/1‑147 (MQ=255)
                                         aaaCACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAAc                                                                                                    <  5:230097/147‑1 (MQ=255)
                                          aaCACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCtt                                                                                                 <  4:296413/149‑1 (MQ=255)
                                            cacac‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGTTATAGATCGAGCGCCGAGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCa                                                                                               <  1:336726/149‑1 (MQ=255)
                                             acac‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAg                                                                                              >  5:232949/1‑149 (MQ=255)
                                               ac‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgcgc                                                                                                                       >  1:73767/1‑122 (MQ=255)
                                               ac‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgcgc                                                                                                                       <  2:73767/122‑1 (MQ=255)
                                                       cTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGccc                                                                                         >  3:142177/1‑145 (MQ=255)
                                                          gtatgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGt                                                                                  >  6:20916/1‑149 (MQ=255)
                                                           tatgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATTGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCAAACCTTCAGCGCCCATCAAGtt                                                                                 <  7:214900/149‑1 (MQ=255)
                                                                   atcATGGACAACCAAACTAATCGGCCCATACCGTTACTGGCGTGAGCTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                         <  2:274616/147‑1 (MQ=38)
                                                                   aGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                         <  5:20916/149‑1 (MQ=255)
                                                                             aCCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGcc                                                               >  3:352943/1‑149 (MQ=255)
                                                                              ccAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                             <  5:203753/134‑1 (MQ=255)
                                                                              ccAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                             >  6:203753/1‑134 (MQ=255)
                                                                                                aCCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCAATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAg                                            >  1:284113/1‑149 (MQ=255)
                                                                                                  cGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCGGCACAACCAGCCCCAGc                                           >  6:206337/1‑148 (MQ=255)
                                                                                                       cTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATtc                                      >  7:194033/1‑148 (MQ=255)
                                                                                                         ggCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTACAGCGCCCATCAAGTTCGATCAGGCAGCTCCGCCCCTGCACAACCCGCCCCAGCAATacat                                    >  6:249339/1‑148 (MQ=255)
                                                                                                          gCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATtcatc                                   <  8:194033/148‑1 (MQ=255)
                                                                                                            ttGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAAATCATCAtt                                >  2:3470/1‑149 (MQ=255)
                                                                                                              gagaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCAtt                                <  6:232949/147‑1 (MQ=255)
                                                                                                                 aTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCtt                           >  7:239985/1‑149 (MQ=255)
                                                                                                                 aTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCAATGCCtt                           >  3:250506/1‑149 (MQ=255)
                                                                                                                    tATGTTGGGGTATAGATCGAGCTCCGCGCGGGCGGCGGTCGATATCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCGCATCAAGTTCGATCAGGCAGCAGCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTcc                         <  5:206337/148‑1 (MQ=255)
                                                                                                                             gTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCtgc                          >  6:287821/1‑136 (MQ=255)
                                                                                                                             gTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTc                          <  5:287821/138‑1 (MQ=255)
                                                                                                                             gTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCAACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGc               <  4:250506/149‑1 (MQ=255)
                                                                                                                                taGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAACCGCAAGCAGt            <  8:239985/149‑1 (MQ=255)
                                                                                                                                  gATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTc           >  5:67718/1‑148 (MQ=255)
                                                                                                                                   aTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCg          >  5:154663/1‑148 (MQ=255)
                                                                                                                                     cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCg                                                                                                       <  7:73988/53‑1 (MQ=255)
                                                                                                                                     cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCg                                                                                                       >  8:73988/1‑53 (MQ=255)
                                                                                                                                     cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCGCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaa        <  6:67718/148‑1 (MQ=255)
                                                                                                                                      gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaac       >  3:27198/1‑148 (MQ=255)
                                                                                                                                       aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaacaa     >  1:22861/1‑149 (MQ=255)
                                                                                                                                          gccgcGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAAccc  <  4:352943/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACAC‑CTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACCC  >  NC_002947/5737457‑5737742

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: ‑