Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,737,655 G→C 12.7% 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,6550GC12.7% 160.7 / 14.3 57intergenic (‑154/+48)hutU/PP_5034urocanate hydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/22);  new base C (5/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

CTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACCCCAAC  >  NC_002947/5737512‑5737802
                                                                                                                                               |                                                                                                                                                   
ctctGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCGTGGACAACCAAAATAATCGGACCTTACCGTTCTTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGAgcgggcg                                                                                                                                               >  1:370605/1‑150 (MQ=255)
         cgtcCGCGGGCGGGTATATCGTCGTCAATATACACTGACACCTGTACTTGTATGTACAAGCATGGACTCTTCAGATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                      <  1:283375/150‑1 (MQ=255)
               cgggcgggTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                <  2:1307474/150‑1 (MQ=255)
                           tCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAAt                                                                                                                    >  2:1473599/1‑150 (MQ=255)
                           tCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAAt                                                                                                                    >  2:533607/1‑150 (MQ=255)
                              tGGTCAATATACAAACACACCTTTACTTGTATGTACAATCATGTACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTa                                                                                                                 <  2:370605/150‑1 (MQ=255)
                              tGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTa                                                                                                                 >  1:841409/1‑150 (MQ=255)
                                           aaCACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTc                                                                                                    >  2:215212/1‑150 (MQ=255)
                                              acacCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGTGCTCGa                                                                                                                                       >  1:1427939/1‑112 (MQ=255)
                                              acacCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGTGCTCGa                                                                                                                                       <  2:1427939/112‑1 (MQ=255)
                                              acacCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAgc                                                                                                 >  2:745443/1‑150 (MQ=255)
                                                    ttACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCAGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCAt                                                                                           >  1:1376471/1‑150 (MQ=255)
                                                       cTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCaa                                                                                        <  1:602722/150‑1 (MQ=255)
                                                          gtatgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGtt                                                                                     <  1:745443/150‑1 (MQ=255)
                                                             tgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGa                                                                                  >  2:1637173/1‑150 (MQ=255)
                                                              gtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                                 >  2:1553483/1‑150 (MQ=255)
                                                              gtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                                 >  1:485060/1‑150 (MQ=255)
                                                               taCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTAc                                                                                                                <  2:1290416/118‑1 (MQ=255)
                                                               taCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTAc                                                                                                                >  1:1290416/1‑118 (MQ=255)
                                                                aCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAg                                                                                                              >  1:377388/1‑119 (MQ=255)
                                                                aCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAg                                                                                                              <  2:377382/119‑1 (MQ=255)
                                                                aCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAg                                                                                                              <  2:377388/119‑1 (MQ=255)
                                                                aCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAg                                                                                                              >  1:377382/1‑119 (MQ=255)
                                                                aCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCCCGcgcggc                                                                                                                                             <  1:1168612/88‑1 (MQ=255)
                                                                aCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCCCGcgcggc                                                                                                                                             >  2:1168612/1‑88 (MQ=255)
                                                                 cAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCACCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAg                                                                              >  1:1004577/1‑150 (MQ=255)
                                                                     cATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                             <  1:1304814/147‑1 (MQ=255)
                                                                     cATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                             >  2:1304814/1‑147 (MQ=255)
                                                                          aCAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCg                                                                     >  2:1432439/1‑150 (MQ=255)
                                                                               cAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcgggcg                                                                                                                                               >  1:508957/1‑71 (MQ=255)
                                                                               cAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcgggcg                                                                                                                                               <  2:508957/71‑1 (MQ=255)
                                                                                aaaaTAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAAt                                                                                                                    <  2:1455480/97‑1 (MQ=255)
                                                                                aaaaTAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAAt                                                                                                                    >  1:1455480/1‑97 (MQ=255)
                                                                                          gACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCAGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGc                                                     <  2:1376471/150‑1 (MQ=255)
                                                                                               tACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAg                                                >  1:403552/1‑150 (MQ=255)
                                                                                                aCCGTTCCTGGCTTGAGATTTATGTTGTGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGc                                               >  2:965399/1‑150 (MQ=255)
                                                                                                            ttGAGATTTATGTTGGGGTGTAGATCTCGCGCGGCGCCGTATCATCTCAATCACACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCCTTg                                   >  2:805127/1‑150 (MQ=255)
                                                                                                            ttGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAg                                                                                                  >  2:856279/1‑87 (MQ=255)
                                                                                                            ttGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAg                                                                                                  <  1:856279/87‑1 (MQ=255)
                                                                                                                gaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAGGTTCGATCAgg                                                                             >  1:1566134/1‑104 (MQ=255)
                                                                                                                gaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                             <  2:1566134/104‑1 (MQ=255)
                                                                                                                  tttATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGACTTcc                             >  2:997288/1‑150 (MQ=255)
                                                                                                                       gTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAAt                                                                                                                    <  2:272710/58‑1 (MQ=255)
                                                                                                                       gTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAAt                                                                                                                    >  1:272710/1‑58 (MQ=255)
                                                                                                                       gTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCt                                >  2:368888/1‑142 (MQ=255)
                                                                                                                       gTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCt                                <  1:368888/142‑1 (MQ=255)
                                                                                                                          ggggTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCaa                     <  1:533607/150‑1 (MQ=255)
                                                                                                                          ggggTATAGATCGAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCaa                     <  2:841409/150‑1 (MQ=255)
                                                                                                                            ctgtgaGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGc                   >  1:181592/6‑150 (MQ=255)
                                                                                                                             gTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAg                                                                                                              <  1:978124/58‑1 (MQ=255)
                                                                                                                             gTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAg                                                                                                              >  2:978124/1‑58 (MQ=255)
                                                                                                                               tgaGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACTACCAGCCCCAGCAATTCATCAtt                                    <  1:877462/128‑1 (MQ=38)
                                                                                                                               tgaGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACTACCAGCCCCAGCAATTCATCAtt                                    >  2:877462/3‑130 (MQ=38)
                                                                                                                               ataGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGt                <  1:215212/150‑1 (MQ=255)
                                                                                                                                  gATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGa             <  2:403552/150‑1 (MQ=255)
                                                                                                                                  gATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGa             >  1:1463072/1‑150 (MQ=255)
                                                                                                                                  gATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGa             >  1:1227665/1‑150 (MQ=255)
                                                                                                                                   aTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaa            <  1:1473599/150‑1 (MQ=255)
                                                                                                                                    tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCTATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaac           <  1:1432439/150‑1 (MQ=255)
                                                                                                                                       aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaacaac        >  1:1421535/1‑150 (MQ=255)
                                                                                                                                             gcgcgGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaacaac        >  2:1362214/1‑144 (MQ=255)
                                                                                                                                             gcgcgGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaacaac        <  1:1362214/144‑1 (MQ=255)
                                                                                                                                             gcgcgGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACCCCAAc  <  2:181592/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACCCCAAC  >  NC_002947/5737512‑5737802

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