Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,780,984 G→A 8.7% C2110C (TGC→TGT)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,780,9840GA 8.7% 128.7 / 2.6 48C2110C (TGC→TGT) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base A (2/2);  total (23/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 = 5.66e-01

GTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACGCTCACCGAG  >  NC_002947/4780837‑4781127
                                                                                                                                                    |                                                                                                                                               
gTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCg                                                                                                                                                 >  6:124255/1‑148 (MQ=255)
                  ggcATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCt                                                                                                                              >  8:226469/1‑149 (MQ=255)
                  ggcATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCt                                                                                                                              <  3:182648/149‑1 (MQ=255)
                  ggcATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCt                                                                                                                              >  3:344822/1‑149 (MQ=255)
                     aTTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCTCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGc                                                                                                                           >  7:209677/1‑149 (MQ=255)
                        gATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTc                                                                                                                         >  4:294109/1‑148 (MQ=255)
                         aTGGGCAGGTGCAGCATCTGTGACGGCGCGGCGTGTCACTGCACCTGCTCACAGACCAACAAGGTAGGCGGCGCTAGGAAGCCTTCCAGCATCGC‑TTTCGCCAAGCCCAGGCAATGCTCATCGAACACCAGCACCCCCGCCTCGCTCa                                                                                                                        >  2:88082/1‑148 (MQ=255)
                          tGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAg                                                                                                                      <  7:226469/149‑1 (MQ=255)
                           gggCAGGTGCATCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAATGTCGGCGGCGCTTGGAAGCCTTCGAGGATCGCTTTGTGCCAAGGCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCaa                                                                                                                       <  1:88082/148‑1 (MQ=255)
                             gCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAgcgc                                                                                                                   <  1:15448/149‑1 (MQ=255)
                               aGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCgcag                                                                                                                 <  5:124255/149‑1 (MQ=255)
                                     agcaTCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATc                                                                                                           <  6:188938/149‑1 (MQ=255)
                                      gcaTCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑CTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCg                                                                                                          <  3:94412/149‑1 (MQ=255)
                                       caTCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCgg                                                                                                         >  3:65861/1‑149 (MQ=255)
                                           tgtgCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAg                                                                                                     >  5:70543/1‑149 (MQ=255)
                                            gtgCCGGCGCGGCGTGTGCCTGCACCTGCGCGCAGACGAAGAAGGGCGGGGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACTCCAGCACGCGCGCCTGGCTCAAGCGCAGCACATCGGCCAg                                                                                                     <  4:239744/148‑1 (MQ=255)
                                               ccGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGgcgc                                                                                                 >  5:54918/1‑149 (MQ=255)
                                                 ggcgcggcgTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGgcgcgc                                                                                               <  2:292596/149‑1 (MQ=255)
                                                 ggcgcggcgTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCACACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCg                                                                                                          >  6:23585/1‑138 (MQ=255)
                                                 ggcgcggcgTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCACACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCg                                                                                                          <  5:23585/138‑1 (MQ=255)
                                                    gcggcgTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCTCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCCGGCGCGCCAc                                                                                            <  8:209677/149‑1 (MQ=255)
                                                            tgCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAg                                                                                     <  6:70543/148‑1 (MQ=255)
                                                              ccTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCCCATCGGCCAGGCGCGCCCCCGGCAAGGAc                                                                                  >  5:235307/1‑149 (MQ=255)
                                                               cTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACg                                                                                 <  3:294109/149‑1 (MQ=255)
                                                                   aCCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATc                                                                             <  4:65861/149‑1 (MQ=255)
                                                                        cTCCCAGACCAACACGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAg                                                                          <  7:199150/147‑1 (MQ=255)
                                                                           ccAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCgg                                                                                         >  5:109591/1‑129 (MQ=255)
                                                                           ccAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCgg                                                                                         <  6:109591/129‑1 (MQ=255)
                                                                                    catggggggcggcTCTTGCAAGCCTTCGAGCATCGT‑TTGGGCAACCCCCAGCAAATGCTAATCCCACACCACCCCCCGCGCCTGGCAAACGCCCACCACTTGGGCCCGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCgg                                                            <  4:344822/142‑1 (MQ=255)
                                                                                            gcggcgCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGGACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGcc                                                    >  5:109606/1‑149 (MQ=255)
                                                                                            gcggcgCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGcc                                                    >  8:66808/1‑149 (MQ=255)
                                                                                                      gAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACAccc                                          >  6:109825/1‑149 (MQ=255)
                                                                                                             atgAGCATTGC‑CTGGGCTTGGCCCAAGCGATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCAGCCTGGAACACCCCGaccac                                   >  1:141995/3‑149 (MQ=38)
                                                                                                                 gCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATg                               >  3:206034/1‑149 (MQ=255)
                                                                                                                      gC‑TTGGGCCAAGCCCAGGCAATGCTCATCGGACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCa                           <  6:109606/148‑1 (MQ=255)
                                                                                                                      gC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCa                           <  7:66808/148‑1 (MQ=255)
                                                                                                                      gC‑TTGGGCCAAGCCCAGGCAATGCTCATCACACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACAcc                                           <  6:266442/132‑1 (MQ=255)
                                                                                                                      gC‑TTGGGCCAAGCCCAGGCAATGCTCATCACACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACAcc                                           >  5:266442/1‑132 (MQ=255)
                                                                                                                             gCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGa                    <  5:109825/149‑1 (MQ=255)
                                                                                                                               cAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGa                    >  5:157798/1‑147 (MQ=255)
                                                                                                                                 aGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTgg                                                 <  7:38596/116‑1 (MQ=255)
                                                                                                                                 aGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTgg                                                 >  8:38596/1‑116 (MQ=255)
                                                                                                                                  gCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGt               >  3:54416/1‑149 (MQ=255)
                                                                                                                                        gCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACGCt         >  5:103648/1‑149 (MQ=255)
                                                                                                                                        ccAATGCTAAGCGCACACACGCCCCCGCCCGGGGCTAACCCGCATCAAAGCGCCCAGGCCCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACGc          <  4:336808/147‑1 (MQ=255)
                                                                                                                                           aTGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACGCTCAc      >  2:289361/1‑149 (MQ=255)
                                                                                                                                            tccTCACCGCACCCCCCCCCCCGCGCCTGGCTCACCCCCACCACATCGGCCCGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACGCTCAcc     <  4:54416/147‑1 (MQ=255)
                                                                                                                                              cTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACGCTCACCGa   <  6:157798/149‑1 (MQ=255)
                                                                                                                                               tCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACGCTCACCGAg  >  5:210893/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGC‑TTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACGCTCACCGAG  >  NC_002947/4780837‑4781127

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