Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,823,341	C→G	8.2%	G2617R (GGC→CGC) ‡	pvdL ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,823,341	0	C	G	8.2%	136.8 / 4.1	49	G2617R (GGC→CGC) ‡	pvdL	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base G (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

TCGAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTACC  >  NC_002947/4823196‑4823488                                                                                                                                                  |                                                                                                                                                    tCGAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGccc                                                                                                                                                    <  6:332666/147‑1 (MQ=255) tCGAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCgg                                                                                                                                                  >  2:165558/1‑149 (MQ=255)   gAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGc                                                                                                                                                 <  8:356122/148‑1 (MQ=255)   gAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGc                                                                                                                                                 >  7:356122/1‑148 (MQ=255)     gcagcaCCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGGGCGGCAAg                                                                                                                                              <  6:7645/149‑1 (MQ=255)       agcaCGAGATGCGCAGGCAGCATGTAGTGGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAACTGCTTGCCGCCCGGCAAGTc                                                                                                                                            <  2:312574/149‑1 (MQ=255)          aCCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAg                                                                                                                                         <  2:107568/149‑1 (MQ=255)              gATGCGCAGGCACCATGTAGTCGGGCACGTGCAGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCACTTGTTGCACCTCGTTACCGGCGGTGGCATGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGAAAGTCCAgcgc                                                                                                                                      >  1:44942/1‑148 (MQ=255)                tGCGCAGGCACCGTGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAg                                                                                                                                   <  5:59979/149‑1 (MQ=255)                tGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAg                                                                                                                                   <  3:177093/149‑1 (MQ=255)                             tGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAg                                                                                                                                              <  2:207877/125‑1 (MQ=255)                                  tCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGAc                                                                                                                 <  7:36502/149‑1 (MQ=255)                                     ggCAGGTGCTGGCGGGCATCGGCTTTCCCCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATc                                                                                                              >  7:197827/1‑149 (MQ=255)                                                   ggCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCa                                                                                                <  1:228314/149‑1 (MQ=255)                                                      aTCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAggc                                                                                             >  6:48192/1‑149 (MQ=255)                                                         ggCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAggcggc                                                                                          >  3:97176/1‑149 (MQ=255)                                                          gCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGct                                                                                         >  8:272399/1‑149 (MQ=255)                                                            tttCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGctct                                                                                       >  2:76494/1‑149 (MQ=255)                                                            tttCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCATGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGctct                                                                                       >  1:258733/1‑149 (MQ=255)                                                                         cgcAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAg                                                                          >  3:365675/1‑149 (MQ=255)                                                                          gcagcaGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGt                                                                         >  7:242373/1‑149 (MQ=255)                                                                              cagTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGa                                                                     >  2:211590/1‑149 (MQ=255)                                                                               agTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGGCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCTGGCGGCTCTCGATTTCGCCCAGTTCGAt                                                                    <  2:44942/149‑1 (MQ=255)                                                                               agTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGCCATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGAt                                                                    >  7:82528/1‑149 (MQ=255)                                                                               agTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGAt                                                                    <  8:197827/149‑1 (MQ=255)                                                                                 ttgttgCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGc                                                                  <  1:76494/149‑1 (MQ=255)                                                                                          gtcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCcgc                                                         <  8:242373/149‑1 (MQ=255)                                                                                          gtcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCcgc                                                         <  3:204810/149‑1 (MQ=255)                                                                                          gtcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCAGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCcgc                                                         <  2:209087/149‑1 (MQ=255)                                                                                          gtagtTACGGGCGGTGGCCTGCGTGCAAAGCTGATAGCCCACCAGCTGGTTGCGGCCCGGCAAGTCGAGCGCCAGCACGAGGTCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCGCAGTTCGATGGGGAAGCcg                                                          <  4:365675/145‑1 (MQ=255)                                                                                             gtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGa                                                      >  8:268616/1‑149 (MQ=255)                                                                                                aCCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATc                                                    <  7:268616/148‑1 (MQ=255)                                                                                                  cggcggTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTc                                                 <  1:165558/149‑1 (MQ=255)                                                                                                    gcggTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCAc                                               <  5:48192/149‑1 (MQ=255)                                                                                                            ctgcctgcACACCAGATAGCCCACCATCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACCGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCTATATCTCCCAGTTCGATGAGGAAGCCGCGGATCTTCACCTGCTGGt                                       >  3:62778/1‑149 (MQ=255)                                                                                                             tgcctgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTc                                      >  5:334420/1‑149 (MQ=255)                                                                                                                    acacCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACgg                               >  6:267548/1‑149 (MQ=255)                                                                                                                        cAGATAGACCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCCTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATGTCGCCCAGTACGATGCGGAAGCCGCGGATCTACACCTGCTGGACGGCCCGGCCg                            >  5:260431/1‑148 (MQ=255)                                                                                                                          gATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACa                         <  4:97176/149‑1 (MQ=255)                                                                                                                             aGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATa                       >  2:322429/1‑148 (MQ=255)                                                                                                                                 caccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTcca                  >  6:46348/1‑149 (MQ=255)                                                                                                                                  accaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTccac                 >  1:148450/1‑149 (MQ=255)                                                                                                                                      gctggatGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACTGCCGACATACTCCACCAgt             <  8:82528/143‑2 (MQ=255)                                                                                                                                        tgctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGc            >  3:171228/1‑148 (MQ=255)                                                                                                                                        gactTGCCGGGCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAg              <  6:88484/144‑1 (MQ=255)                                                                                                                                        gactTGCCGGGCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAg              >  5:88484/3‑146 (MQ=255)                                                                                                                                        gactTGCCGGGCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGc            >  7:57127/3‑148 (MQ=255)                                                                                                                                         gctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATAc                      <  7:241760/137‑1 (MQ=255)                                                                                                                                         gctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATAc                      >  8:241760/1‑137 (MQ=255)                                                                                                                                         gctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCa          >  3:27871/1‑149 (MQ=255)                                                                                                                                                gccCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTAc   >  1:241984/1‑149 (MQ=255)                                                                                                                                                  ccGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTAcc  >  4:12555/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TCGAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTACC  >  NC_002947/4823196‑4823488

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