Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,826,490 Δ1 bp 12.2% coding (4745/12954 nt) pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,826,4900T.12.2% 113.2 / 11.8 41coding (4745/12954 nt)pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (18/18);  new base . (3/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.41e-01

TGTTGCAGCTCGGTTCGTGGCTCGACATGCTCACGCTGCTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGG  >  NC_002947/4826343‑4826635
                                                                                                                                                   |                                                                                                                                                 
tgttgCAGCTCGGTTCGTGGCTCGACATGCTCACGCTGCTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCa                                                                                                                                                 <  2:584342/150‑1 (MQ=255)
tgttgCAGCTCGGTTCGTGGCTCGACATGCTCACGCTGCTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCa                                                                                                                                                 <  1:508714/150‑1 (MQ=255)
                            gCTCACGCTGCTGCCACACTGGGGCAGGCAGTATGGGTCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATCTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCt                                                                                                                     >  1:1188956/1‑150 (MQ=255)
                             cTCACGCTGCTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTg                                                                                                                    <  1:371820/150‑1 (MQ=255)
                              tCACGCTGCTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTgg                                                                                                                   <  2:756623/150‑1 (MQ=255)
                                   tgcctgcCCCAGTGTGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCCGGGGCCCAATGGCATCTGTGCCAAAGGGCTGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCt                                                                                                              >  2:402924/4‑150 (MQ=255)
                                               tGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATa                                                                                                  >  2:38114/1‑150 (MQ=255)
                                                   gcaggcagTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCg                                                                                              >  1:1098179/1‑150 (MQ=255)
                                                       gcagTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCa                                                                                          >  1:1027066/1‑150 (MQ=255)
                                                       gcagTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCa                                                                                          >  2:757261/1‑150 (MQ=255)
                                                        cagTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGc                                                                                                                                      <  1:1058004/105‑1 (MQ=255)
                                                        cagTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGc                                                                                                                                      >  2:1058004/1‑105 (MQ=255)
                                                              cTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTgcgc                                                                                   <  1:38114/150‑1 (MQ=255)
                                                                      cGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATgg                                                                           <  1:757261/150‑1 (MQ=255)
                                                                           ttGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTcc                                                                      >  2:908039/1‑150 (MQ=255)
                                                                             gCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCtt                                                                    <  1:776614/150‑1 (MQ=255)
                                                                               cGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTtgt                                                                  <  2:417517/150‑1 (MQ=255)
                                                                               cGCTGGGGACCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTtgt                                                                  >  2:380311/1‑150 (MQ=255)
                                                                                   ggggCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTAt                                                              <  1:908039/150‑1 (MQ=255)
                                                                                    gggCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATc                                                             <  2:1027066/150‑1 (MQ=255)
                                                                                    gggCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATc                                                             >  2:228835/1‑150 (MQ=255)
                                                                                              gCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTtg                                                                   >  1:163393/1‑134 (MQ=255)
                                                                                              gCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTtg                                                                   <  2:163393/134‑1 (MQ=255)
                                                                                                aTCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGcc                                                 >  2:28781/1‑150 (MQ=255)
                                                                                                                 gAGTTGCGCCGGCACCAGGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCtg                                <  2:1098179/150‑1 (MQ=255)
                                                                                                                   gTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCtgtg                              >  1:140287/1‑150 (MQ=255)
                                                                                                                    ttGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGc                                                                                                >  1:458210/1‑83 (MQ=255)
                                                                                                                    ttGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGc                                                                                                <  2:458210/83‑1 (MQ=255)
                                                                                                                                    tggctgagCAATTGCCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCa            >  1:995640/8‑150 (MQ=255)
                                                                                                                                    tggctgagCAATTGCCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCa            >  1:164616/8‑150 (MQ=255)
                                                                                                                                     aGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTtgt                                                                  >  2:617340/1‑96 (MQ=255)
                                                                                                                                     aGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTtgt                                                                  <  1:617340/96‑1 (MQ=255)
                                                                                                                                        cTGGCAATTGCTCAGCGAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCCGCAGGCAGGCCTGCACtt         >  1:41545/1‑150 (MQ=255)
                                                                                                                                         tGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTc        <  1:437107/150‑1 (MQ=255)
                                                                                                                                         gagCAATTGCCAGCCAGTACGGCCAGCACCTCGGCAGGCTggg                                                                                                                  >  2:387131/3‑43 (MQ=255)
                                                                                                                                         gagCAATTGCCAGCCAGTACGGCCAGCACCTCGGCAGGCTggg                                                                                                                  <  1:387131/41‑1 (MQ=255)
                                                                                                                                           gCAATTGCCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTc     <  2:164616/150‑1 (MQ=255)
                                                                                                                                             aaTTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCg    >  2:501747/1‑150 (MQ=255)
                                                                                                                                              aTTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTtgt                                                                  >  2:671815/1‑87 (MQ=255)
                                                                                                                                              aTTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTtgt                                                                  <  1:671815/87‑1 (MQ=255)
                                                                                                                                               ttGCTCACCAAGTACGGGTAGCGCTTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTAACATCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCggg  <  1:380311/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGTTGCAGCTCGGTTCGTGGCTCGACATGCTCACGCTGCTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATGTAGGCTGGCAATTGCTCAGCCAGTACGGCCAGCACCTCGGCAGGCTGGGCGCTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGG  >  NC_002947/4826343‑4826635

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