Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,379,928 T→C 8.1% F505L (TTC→CTC)  recN → DNA repair/recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,379,9280TC 8.1% 135.2 / 2.8 49F505L (TTC→CTC) recNDNA repair/recombination protein
Reads supporting (aligned to +/- strand):  ref base T (21/24);  new base C (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

CACCCTGGTGTTCGACGAAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGG  >  NC_002947/5379780‑5380072
                                                                                                                                                     |                                                                                                                                                
cACCCTGGTGTTCGACGAAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGt                                                                                                                                                  <  3:261218/149‑1 (MQ=255)
 aCCCTGGTGTTCGACGAAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGtt                                                                                                                                                 <  7:149989/149‑1 (MQ=255)
    cTGGTGTTCGACGAAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGCTCGt                                                                                                                                              >  4:175140/1‑149 (MQ=255)
    cTGGTGTTCGACGAAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGCTCGt                                                                                                                                              >  4:169146/1‑149 (MQ=255)
       gtgtTCGACGAAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGTCCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTg                                                                                                                                             <  1:82759/147‑1 (MQ=255)
          ttCGACGAAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTg                                                                                                                                             >  1:184392/1‑144 (MQ=255)
          ttCGACGAAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTg                                                                                                                                             <  2:184392/144‑1 (MQ=255)
            cgacgaAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAg                                                                                                                                       <  8:6078/148‑1 (MQ=255)
             gacgaAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGt                                                                                                                                     <  4:361618/149‑1 (MQ=255)
                gaAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTgcg                                                                                                                                  >  1:301275/1‑149 (MQ=255)
                      gACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAg                                                                                                                            <  5:169592/149‑1 (MQ=255)
                                        ccAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCCGTGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGacac                                                                                                                      >  1:257928/1‑138 (MQ=255)
                                        ccAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGacac                                                                                                                      <  2:257928/137‑1 (MQ=255)
                                         cAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACAccgcc                                                                                                         >  5:156667/1‑149 (MQ=255)
                                         cAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACAccgcc                                                                                                         <  6:42869/149‑1 (MQ=255)
                                                       gtcgGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTg                                                                                            <  5:346613/148‑1 (MQ=255)
                                                              aaCTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGt                                                                                    <  6:156667/149‑1 (MQ=255)
                                                                    tgcGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAgcgc                                                                              >  2:256746/1‑149 (MQ=255)
                                                                     gcggcgCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGcgcg                                                                             >  7:35052/1‑149 (MQ=255)
                                                                     gcggcgCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCATGGCAAGCGTGAGcccg                                                                             >  4:227460/1‑149 (MQ=255)
                                                                       ggcgCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTg                                                                           >  6:350507/1‑149 (MQ=255)
                                                                          gcC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGgaa                                                                        <  3:48873/149‑1 (MQ=255)
                                                                               gggCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGaagaag                                                                    >  6:302159/1‑149 (MQ=255)
                                                                                ggCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGgaagaa                                                                     <  6:133543/147‑1 (MQ=255)
                                                                                ggCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTACGTGCACAAGGTGCGCTACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGaagaag                                                                    <  1:120589/148‑1 (MQ=255)
                                                                                               gTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGgaa                                                                        <  3:249597/129‑1 (MQ=255)
                                                                                               gTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGgaa                                                                        >  4:249597/1‑129 (MQ=255)
                                                                                               gTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTcgg                                                    >  5:279270/1‑149 (MQ=255)
                                                                                                    gACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTgg                                                                          <  6:188256/122‑1 (MQ=255)
                                                                                                     aCCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGa                                               <  1:87551/148‑1 (MQ=255)
                                                                                                     aCCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGCTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGa                                               <  3:169146/148‑1 (MQ=255)
                                                                                                     aCCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGCTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGa                                               <  3:175140/148‑1 (MQ=255)
                                                                                                      ccGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATc                                             <  5:350507/149‑1 (MQ=255)
                                                                                                        gTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGa                                           >  4:278112/1‑149 (MQ=255)
                                                                                                         tGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGAc                                          >  2:198812/1‑149 (MQ=255)
                                                                                                           aCCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCt                                        >  7:123671/1‑149 (MQ=255)
                                                                                                            cccACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCt                                        <  1:263998/148‑1 (MQ=255)
                                                                                                             ccACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACGGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTGGGGGGGATCGCCCttc                                      >  6:261011/1‑147 (MQ=255)
                                                                                                             ccACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGa                                      >  5:134598/1‑149 (MQ=255)
                                                                                                             ccACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGa                                      >  4:325110/1‑149 (MQ=255)
                                                                                                             ccACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGa                                      <  6:279270/149‑1 (MQ=255)
                                                                                                             ccACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGa                                      >  4:162410/1‑149 (MQ=255)
                                                                                                                     cGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGGaa                              <  8:123671/149‑1 (MQ=255)
                                                                                                                         ggTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGGAATccc                          <  6:134598/149‑1 (MQ=255)
                                                                                                                             gcggcgCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTgg                       <  1:256746/148‑1 (MQ=255)
                                                                                                                               ggcgCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGTAATCCCTGGccc                    <  5:15721/149‑1 (MQ=255)
                                                                                                                               ggcgCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTGGcc                     >  3:261799/1‑148 (MQ=255)
                                                                                                                                            cACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCCCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGAt       >  6:106880/1‑149 (MQ=255)
                                                                                                                                             aCCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTTGGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGAtg      <  2:301275/149‑1 (MQ=255)
                                                                                                                                                aTCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTTGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGAtggtg   <  5:302159/149‑1 (MQ=255)
                                                                                                                                                 tCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGAtggtgg  >  5:240971/1‑149 (MQ=255)
                                                                                                                                                     ttCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGAtggtgg  <  4:261799/145‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CACCCTGGTGTTCGACGAAGTCGACGTGGGTATTGGCGGGCCAACGGCAGAAATCGTCGGCCAACTGCTGCGGCGCC‑TGGGCGAACGCGGCCAGGTACTGACCGTGACCCACCTGCCGCAGGTGGCGGCGCAGGGGCATCACCATCTGTTCGTGCACAAGGTGCGCAACAGCGACACCACCCACACCGCCGTGGCTAGCCTTGGCAAGCGTGAGCGCGTGGAAGAAGTGGCGCGGATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGG  >  NC_002947/5379780‑5380072

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