Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,783,688 A→G 7.3% L1209P (CTG→CCG)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,783,6880AG 7.3% 158.1 / 2.6 55L1209P (CTG→CCG) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base G (2/2);  total (29/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.50e-01

TGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGG  >  NC_002947/4783541‑4783836
                                                                                                                                                     |                                                                                                                                                    
tGTCGCGCCCATTGGCGCAGGGGTGCCGACGTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAg                                                                                                                                                     >  1:21930/1‑149 (MQ=255)
    gcgcCCATTGGCTCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGt                                                                                                                                                 <  7:97120/149‑1 (MQ=255)
    gcgcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGt                                                                                                                                                 >  2:30336/1‑149 (MQ=255)
    gcgcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGt                                                                                                                                                 >  1:34715/1‑149 (MQ=255)
      gcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTg                                                                                                                                               >  8:3420/1‑149 (MQ=255)
        ccATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCt                                                                                                                                             >  6:132436/1‑149 (MQ=255)
        ccATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCt                                                                                                                                             <  2:190433/149‑1 (MQ=255)
           ttGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGAc                                                                                                                                          <  1:104878/149‑1 (MQ=255)
           ttGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGAc                                                                                                                                          <  2:34715/149‑1 (MQ=255)
              gcgcAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCACTAGGCTAGCCTGATGCTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCAt                                                                                                                                       >  2:38024/1‑149 (MQ=255)
                gcAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCAACg                                                                                                                                     >  7:201659/1‑149 (MQ=255)
                         ccGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCg                                                                                                                            <  2:7940/149‑1 (MQ=255)
                          cGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCg                                                                                                                            >  3:106189/1‑148 (MQ=255)
                              cTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACaa                                                                                                                       >  1:234696/1‑149 (MQ=255)
                                       aGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCt                                                                                                              <  2:234696/149‑1 (MQ=255)
                                          cTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGaa                                                                                                           <  2:40512/149‑1 (MQ=255)
                                             aCCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAAc                                                                                                          <  2:21930/147‑1 (MQ=255)
                                                     cccTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCAc                                                                                                >  5:265000/1‑149 (MQ=255)
                                                      ccTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCTGGCCCCCAGTTCATGCACCGCCAGCAACACCGGTTGATCACAAACGGCCATTTCGAACAATGACCCCCTGAACACCTCGcacc                                                                                                >  1:80990/1‑146 (MQ=255)
                                                           gCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGc                                                                                            >  5:157339/1‑147 (MQ=255)
                                                                 ctgctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGa                                                                                    <  7:3420/149‑1 (MQ=255)
                                                                  tgctgGTAGGCCGTCTGCATGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCATCAACCAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGaa                                                                                                           <  1:191824/126‑1 (MQ=255)
                                                                  tgctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGaa                                                                                                           >  2:191824/1‑125 (MQ=255)
                                                                    ctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAgg                                                                                  <  6:157339/148‑1 (MQ=255)
                                                                     tgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCt                                                                                >  1:72870/1‑149 (MQ=255)
                                                                       gTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCt                                                                                <  6:265000/147‑1 (MQ=255)
                                                                       gTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTg                                                                               <  2:223524/148‑1 (MQ=255)
                                                                                  gCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTc                                                                    <  4:106189/148‑1 (MQ=255)
                                                                                      gTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTgg                                                               <  5:132436/149‑1 (MQ=255)
                                                                                      gTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGAAccaagtttc                                                               >  8:94141/1‑142 (MQ=255)
                                                                                       tCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTggc                                                              <  2:72870/149‑1 (MQ=255)
                                                                                              aCAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAg                                                                                         >  3:210553/1‑115 (MQ=255)
                                                                                              aCAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAg                                                                                         <  4:210553/115‑1 (MQ=255)
                                                                                                   aTGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAg                                                  <  1:30336/149‑1 (MQ=255)
                                                                                                       gcCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCt                                              >  1:182962/1‑149 (MQ=255)
                                                                                                       gcCAAGACACCCCATCCACCACCAGGAGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCt                                              >  7:155879/1‑149 (MQ=255)
                                                                                                        ccAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTg                                             >  2:64434/1‑149 (MQ=255)
                                                                                                             acacCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCAACGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAgg                                        <  8:201659/149‑1 (MQ=255)
                                                                                                               acCCCATCCACGACCAGGTGATGCACCGCCAGCAACAGCCGGTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCGACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGGGGTTTCTAGCGCTGCAAGGCt                                      <  8:155879/149‑1 (MQ=255)
                                                                                                               acCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCt                                      >  6:52809/1‑149 (MQ=255)
                                                                                                                 cccATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGt                                    >  7:153937/1‑149 (MQ=255)
                                                                                                                    aTCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCa                                 <  1:64434/149‑1 (MQ=255)
                                                                                                                        caccaccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCagag                              >  3:163096/1‑149 (MQ=255)
                                                                                                                             ccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGAc                          >  1:36465/1‑148 (MQ=255)
                                                                                                                             ccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACt                         <  2:36465/149‑1 (MQ=255)
                                                                                                                               aGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTAc                       >  2:159890/1‑149 (MQ=255)
                                                                                                                                ggTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAgc                                                 <  3:128257/122‑1 (MQ=255)
                                                                                                                                ggTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAgc                                                 >  4:128257/1‑122 (MQ=255)
                                                                                                                                   gATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTg                    >  4:244273/1‑148 (MQ=255)
                                                                                                                                            tcAGCAACGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGgcggcg                                                          >  6:243603/2‑101 (MQ=255)
                                                                                                                                            tcAGCAACGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGgcggcg                                                          <  5:243603/100‑1 (MQ=255)
                                                                                                                                            tcAGCAACGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCaggagg             <  3:143722/145‑1 (MQ=255)
                                                                                                                                            tcAGCAACGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCaggagg             >  4:143722/2‑146 (MQ=255)
                                                                                                                                              aGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTg        <  7:94141/149‑1 (MQ=255)
                                                                                                                                               gCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTg        >  8:18872/1‑148 (MQ=255)
                                                                                                                                                     aGCCGTTGCTGCCCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGGTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCgg  >  6:141024/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGG  >  NC_002947/4783541‑4783836

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: