Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 463,678 G→A 7.8% F473F (TTC→TTT)  davB ← lysine 2‑monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947463,6780GA 7.8% 145.7 / 2.2 51F473F (TTC→TTT) davBlysine 2‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (23/24);  new base A (2/2);  total (25/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 = 5.69e-01

TCTGCACCGCGCCTTCAACCCAGGCGGGGGTCCATGACACGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCA  >  NC_002947/463530‑463823
                                                                                                                                                    |                                                                                                                                                 
tCTGCACCGCGCCTTCAACCCAGGCGGGGGTCCATGACACGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTg                                                                                                                                                   <  8:259121/149‑1 (MQ=255)
    cACCGCGCCTTCAACCCAGGCGGGGGTCCATGACACGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAgg                                                                                                                                               <  3:67760/149‑1 (MQ=255)
       cgcgCCTTCAACCCAGGCGGGGGTCCATGACACGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTCAAAGgcgc                                                                                                                                            <  2:67008/149‑1 (MQ=255)
            cTTCAACCCAGGCGGGGGTCCATGACACGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCGGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTCAAAGgcgcgcttc                                                                                                                                       <  5:6288/149‑6 (MQ=255)
                    cAGGCGGGGGTCCATGACACGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACgg                                                                                                                               >  3:60178/1‑149 (MQ=255)
                     acgCGGGGGTCCATGACACGTCGTCACCGGCATTGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACggg                                                                                                                              <  1:285966/147‑1 (MQ=255)
                     aGGCGGGGGTCCATGACACGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACggg                                                                                                                              >  2:303018/1‑149 (MQ=255)
                        cGGGGGTCCATGACACGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCg                                                                                                                           <  2:215057/149‑1 (MQ=255)
                                    acacGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGa                                                                                                                <  1:185506/148‑1 (MQ=255)
                                     cacGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGaaa                                                                                                              <  1:330770/149‑1 (MQ=255)
                                      acGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAAc                                                                                                             >  3:70314/1‑149 (MQ=255)
                                       cgtcgtCACCTGCAATGAAGATACGGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAATGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAAcc                                                                                                            <  2:185762/149‑1 (MQ=255)
                                        gtcgtcACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCg                                                                                                           <  5:172871/149‑1 (MQ=255)
                                         tcgtcACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTATGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGt                                                                                                          >  1:110277/1‑149 (MQ=255)
                                              aCCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCg                                                                                                     >  1:43940/1‑149 (MQ=255)
                                               ccGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCGAAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCgg                                                                                                    <  6:250274/149‑1 (MQ=255)
                                               ccGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCgg                                                                                                    <  4:70314/149‑1 (MQ=255)
                                                              ccGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGAtgt                                                                                      <  4:60178/148‑1 (MQ=255)
                                                               cgcgCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGaa                                                                                                                                     <  8:307781/100‑1 (MQ=255)
                                                                    tgctCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGc                                                                               >  3:40527/1‑149 (MQ=255)
                                                                       tCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGa                                                                             >  8:268040/1‑148 (MQ=255)
                                                                          gccgGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATc                                                                         <  5:158509/149‑1 (MQ=255)
                                                                           ccgGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCg                                                                        <  4:40527/149‑1 (MQ=255)
                                                                            cgGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCgg                                                                       <  6:172542/149‑1 (MQ=255)
                                                                            cgGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCgg                                                                       <  1:303018/149‑1 (MQ=255)
                                                                             ggCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAg                                                                                                                  >  1:222364/1‑105 (MQ=255)
                                                                             ggCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAg                                                                                                                  <  2:222364/105‑1 (MQ=255)
                                                                                   tCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCgg                                                                >  6:323172/1‑149 (MQ=255)
                                                                                    ccTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCggg                                                               >  8:138991/1‑149 (MQ=255)
                                                                                         tgcATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGAGGATGTGCCCGGCGATATCGGTCTTCGGGTAGAt                                                          >  8:244644/1‑149 (MQ=255)
                                                                                           cATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGAtct                                                        >  6:313567/1‑149 (MQ=255)
                                                                                             tGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGAtcttc                                                      <  5:313567/149‑1 (MQ=255)
                                                                                                 gTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCa                                                  <  2:43940/149‑1 (MQ=255)
                                                                                                      cgTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCa                                             <  5:323172/149‑1 (MQ=255)
                                                                                                       gTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCAt                                            >  1:129873/1‑149 (MQ=255)
                                                                                                       gTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATAGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCa                                             >  3:8674/1‑148 (MQ=255)
                                                                                                           aTGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAg                                        >  5:58578/1‑149 (MQ=255)
                                                                                                             gcgcTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAggg                                      >  6:106022/1‑149 (MQ=255)
                                                                                                               gcTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGc                                     >  4:137821/1‑148 (MQ=255)
                                                                                                                   gTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGtt                                <  2:110277/149‑1 (MQ=255)
                                                                                                                      gTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTa                             >  1:225933/1‑149 (MQ=255)
                                                                                                                       tAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTAc                            <  7:268040/149‑1 (MQ=255)
                                                                                                                               aTGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGACGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCtttt                    <  4:8674/149‑1 (MQ=255)
                                                                                                                                tGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTc                   >  2:215136/1‑149 (MQ=255)
                                                                                                                                 ggCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTcc                  >  5:317049/1‑149 (MQ=255)
                                                                                                                                    ccGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCAcc               >  6:107774/1‑149 (MQ=255)
                                                                                                                                         tcggcgcCTTCAAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCTGGTg          >  8:114832/4‑149 (MQ=255)
                                                                                                                                         tcggcgcCTTCAAAGGCCCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTg          >  6:277095/4‑149 (MQ=255)
                                                                                                                                          gcgcgcCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGc         <  7:138991/149‑1 (MQ=255)
                                                                                                                                           cgcgcCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCg        <  6:317049/149‑1 (MQ=255)
                                                                                                                                              gcCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCgg       >  8:124847/1‑147 (MQ=255)
                                                                                                                                                 tttGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGACTCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGgcagca  >  4:293024/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCTGCACCGCGCCTTCAACCCAGGCGGGGGTCCATGACACGTCGTCACCGGCAATGAAGATACCGCGCTGCTCCGCCGGCATGTCCTGCTGCATGAAGTGCGCGTACATGCGCTGGTTGTAGCGGTAATGGCCCGGAAGCGCGCCTTTGAAGGCGCCGAGGAAGTACGGGTCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCA  >  NC_002947/463530‑463823

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