Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,484 T→C 12.0% intergenic (+49/+22) hutF → / ← PP_5037 formiminoglutamate deiminase/lipocalin family lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,740,4840TC12.0% 122.8 / 8.4 50intergenic (+49/+22)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base C (3/3);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

ATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTG  >  NC_002947/5740337‑5740626
                                                                                                                                                    |                                                                                                                                              
aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                               <  3:254551/149‑1 (MQ=255)
aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGt                                                                                                                                                <  6:123274/148‑1 (MQ=255)
 tGGTGAATGGGCAATGGGGGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg                                                                                                                                              <  6:290079/149‑1 (MQ=38)
                    ggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGCGTTGGGGCCCTAGGGGGGGACGCCCCGCACGGGTTGGACCGCCACCCCcaaaa                                                                                                                             >  1:125982/1‑148 (MQ=255)
                       gcgcCTGGGGGGGCTTGCCAAGTATTATTATCGGGGGCGGGCTTTTCCAATGGTGTTGGGCAAGCTTTTGGTCTGTGTTTAGGGGGTTT‑GGGTTCCAGCCCTTGGCGTGCAAGCCGGCTATGGGTTCAAACGCAACCCCAAAAAACTCa                                                                                                                        <  4:139346/149‑1 (MQ=255)
                           cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                              <  8:12027/123‑1 (MQ=255)
                           cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                              >  7:12027/1‑123 (MQ=255)
                                      tGCCCAGGAATATGACAGCGGGCGGGCGTTTGCAAAGGTGTTGCGCCAGCTTTTGGGCTGCGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAt                                                                                                         <  5:294683/149‑1 (MQ=255)
                                           aGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTAGCGTGGAAGCCCGCGAAGGGTTTGAACGCAACCCCAAAAAAATCAGTGCATTTTTACTAttttgt                                                                                                    >  5:61143/1‑146 (MQ=255)
                                                  tGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGATTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCATCCCTTCGCTTGAAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGa                                                                                             <  5:181877/149‑1 (MQ=255)
                                                         gggcgggcGTTTGCCTAGGGGTTGCTCCTGCTTTTGGGCTGAGCTTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCACGCCGGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgc                                                                                      <  6:126466/149‑1 (MQ=255)
                                                                 gTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTTCAGCCCTTAGCGTGCAAGCCCCCGAAGGGTTGGAACGCAACCCCAAAAAACAGAGGGGATTTTTAATATTTGCTGAAAGGACGCCCCCCCGATCaa                                                                              >  6:158771/1‑149 (MQ=255)
                                                                       cAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGGTTCAGCCCATAGCGGGGAAGCCCGCGAACGGTTTGAACCCAAACCCCAAAAAAAAACGGCATTTTTTATATTTTGTGAAAGGCCCCCCGCCACAGAAAGCGGcg                                                                        >  2:208472/1‑148 (MQ=255)
                                                                        aaGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                       <  1:208472/149‑1 (MQ=255)
                                                                             gtTGCGCCAGCTTTTGGGCTGAGCGTACCGNGGTT‑GCGTTGCAGCCCTTCGCGTGGAAGCCCGCGAAGGGTTGGGACGCAACCCCAAAAAACTCAGTGCATTTTTAGTATTTGGTGATGGGGCGCCCGCCCGATGAAGCGGgtggtgg                                                                   >  1:27501/1‑147 (MQ=255)
                                                                                 cgcCTGCTTTTGGGGTCCGCGACCCGGGGTG‑GCGTGGCAGCCTTCCGCGTCCACGCCCGATAAGGGTTCAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc                                                              <  2:325968/149‑1 (MQ=255)
                                                                                 cgcCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc                                                              <  2:100901/149‑1 (MQ=255)
                                                                                   ccAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                                            <  5:158771/149‑1 (MQ=255)
                                                                                                     gTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                          <  8:208215/149‑1 (MQ=255)
                                                                                                              tt‑gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGtt                                 >  5:78056/1‑149 (MQ=255)
                                                                                                                   gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                             <  6:78056/149‑1 (MQ=255)
                                                                                                                   gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                             >  7:137137/1‑149 (MQ=38)
                                                                                                                   gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                             >  6:187726/1‑149 (MQ=38)
                                                                                                                    ttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                            >  4:174085/1‑149 (MQ=255)
                                                                                                                      gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcag                                       <  4:312109/136‑1 (MQ=255)
                                                                                                                      gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcag                                       >  3:312109/1‑136 (MQ=255)
                                                                                                                       caacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                            >  1:214526/4‑146 (MQ=39)
                                                                                                                       caacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                            <  2:214526/143‑1 (MQ=39)
                                                                                                                        aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGa                                                                                  >  5:165934/3‑91 (MQ=21)
                                                                                                                        aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGa                                                                                  <  6:165934/89‑1 (MQ=21)
                                                                                                                        aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                                   <  5:206733/138‑1 (MQ=255)
                                                                                                                        aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                                   >  6:206733/1‑138 (MQ=255)
                                                                                                                        aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc                        >  6:13583/1‑149 (MQ=255)
                                                                                                                         gCCCTTCGCGTGCAAGCCCGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                       <  1:324200/149‑1 (MQ=255)
                                                                                                                          cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCCACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCTGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCGcc                       >  6:39569/1‑148 (MQ=255)
                                                                                                                          cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                            <  3:296459/127‑1 (MQ=255)
                                                                                                                          cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                            >  4:296459/1‑127 (MQ=255)
                                                                                                                            cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                    <  8:137137/149‑1 (MQ=39)
                                                                                                                                 cgTGCATCCCCGCTATGGGTTCAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg               <  4:175931/149‑1 (MQ=255)
                                                                                                                                    gCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAg            >  5:123401/1‑149 (MQ=255)
                                                                                                                                    gCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAg            >  6:74912/1‑149 (MQ=255)
                                                                                                                                      aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgctg                                                      >  7:302256/1‑105 (MQ=255)
                                                                                                                                      aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgctg                                                      <  8:302256/105‑1 (MQ=255)
                                                                                                                                      aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCa             <  6:71257/146‑1 (MQ=255)
                                                                                                                                      aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCa             >  5:71257/1‑146 (MQ=255)
                                                                                                                                            gcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt         >  3:156841/1‑144 (MQ=255)
                                                                                                                                            gcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt         <  4:156841/144‑1 (MQ=255)
                                                                                                                                             cgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtgt   >  7:211219/1‑149 (MQ=255)
                                                                                                                                             cgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtgt   <  6:123401/149‑1 (MQ=255)
                                                                                                                                              gAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGtgtg  >  4:134084/1‑149 (MQ=255)
                                                                                                                                                   gTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg       <  7:160352/139‑1 (MQ=255)
                                                                                                                                                   gTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg       >  8:160352/1‑139 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTG  >  NC_002947/5740337‑5740626

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