Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	583,595	A→G	V939A (GTC→GCC)	nfrA ←	bacteriophage adsorption protein NfrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	583,595	0	A	G	95.5%	140.7 / ‑0.8	44	V939A (GTC→GCC)	nfrA	bacteriophage adsorption protein NfrA
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (23/19);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTG  >  NZ_CP009273/583460‑583728                                                                                                                                        |                                                                                                                                      agaAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTc                                                                                                                                    >  1:563521/1‑139 (MQ=255)         gttgttGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGc                                                                                                                            >  2:684568/1‑139 (MQ=255)                          tACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGcc                                                                                                           >  2:510429/1‑139 (MQ=255)                             gTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCttt                                                                                                        <  1:337962/139‑1 (MQ=255)                                 aTTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAg                                                                                                               <  1:530049/128‑1 (MQ=255)                                 aTTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAg                                                                                                               >  2:530049/1‑128 (MQ=255)                                 aTTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCa                                                                                                    >  2:527160/1‑139 (MQ=255)                                          cTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCg                                                                                           <  2:649839/139‑1 (MQ=255)                                           ttAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTAt                                                                                                      >  1:341089/1‑127 (MQ=255)                                           ttAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTAt                                                                                                      <  2:341089/127‑1 (MQ=255)                                                    tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGa                                                                                                                <  1:50026/108‑1 (MQ=255)                                                    tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGa                                                                                                                >  2:50026/1‑108 (MQ=255)                                                    tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGAt                                                                                                >  2:591188/1‑124 (MQ=255)                                                    tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGAt                                                                                                <  1:591188/124‑1 (MQ=255)                                                    tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGaa                                                                                 >  2:576455/1‑139 (MQ=255)                                                    tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGaa                                                                                 >  2:230271/1‑139 (MQ=255)                                                    tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGaa                                                                                 >  2:11754/1‑139 (MQ=255)                                                               gACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTa                                                                       <  1:510429/138‑1 (MQ=255)                                                                aCGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAg                                                                     >  1:544982/1‑139 (MQ=255)                                                                 cGCCGAGACTGACTTTGGGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCCCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGCAAgg                                                                    >  2:497742/1‑139 (MQ=255)                                                                             ctTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATACTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTg                                                        <  2:544982/139‑1 (MQ=255)                                                                                 gtgCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCg                                                    >  1:634659/1‑139 (MQ=255)                                                                                      gCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGc                                               >  2:269625/1‑139 (MQ=255)                                                                                           gCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCt                                          >  2:156055/1‑139 (MQ=255)                                                                                                gtAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGa                                     >  1:339804/1‑139 (MQ=255)                                                                                                     gCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCg                                                                                           <  1:79767/80‑1 (MQ=255)                                                                                                     gCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCg                                                                                           >  2:79767/1‑80 (MQ=255)                                                                                                      cGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAg                               >  1:217072/1‑139 (MQ=255)                                                                                                        tctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCg                                                                                           >  2:139910/1‑77 (MQ=255)                                                                                                        tctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCg                                                                                           <  1:139910/77‑1 (MQ=255)                                                                                                        tctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCt                             >  2:493381/1‑139 (MQ=255)                                                                                                        tctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCt                             <  1:527160/139‑1 (MQ=255)                                                                                                         ctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGt                                                                              >  1:765761/1‑89 (MQ=255)                                                                                                         ctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGt                                                                              <  2:765761/89‑1 (MQ=255)                                                                                                            gCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGAcg                         <  1:684568/139‑1 (MQ=255)                                                                                                            gCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGAcg                         <  2:560305/139‑1 (MQ=255)                                                                                                            gCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGAcg                         <  1:450996/139‑1 (MQ=255)                                                                                                            gCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGACAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGa                           <  1:657835/137‑1 (MQ=255)                                                                                                                            cAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTc          <  2:353134/138‑1 (MQ=255)                                                                                                                            cAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTc          <  1:558349/138‑1 (MQ=255)                                                                                                                               cGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTa                                                                       <  2:511092/74‑1 (MQ=255)                                                                                                                               cGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTa                                                                       >  1:511092/1‑74 (MQ=255)                                                                                                                                    cccGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTg  >  2:350058/1‑138 (MQ=255)                                                                                                                                    cccGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTg  >  1:555906/1‑138 (MQ=255)                                                                                                                                        |                                                                                                                                      AGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTG  >  NZ_CP009273/583460‑583728

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