Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,451,094	G→T	11.4%	A166S (GCC→TCC)	gpD →	tail formation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,451,094	0	G	T	11.4%	114.1 / 3.5	35	A166S (GCC→TCC)	gpD	tail formation
Reads supporting (aligned to +/- strand):  ref base G (17/14);  new base T (2/2);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

CTGGGAGGTGTCTTGCCCGGTCGACACCAAGGTCGAGCGCGTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTCCCGCTACCAGTTCCGCCTCGGTGATCGCAACTCGCAGAA  >  NC_002947/3450979‑3451227                                                                                                                    |                                                                                                                                      cTGGGAGGTGTCTTGCCCGGTCGACACCAAGGTCGAGCGCGTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTg                                                                                                                   <  1:273261/136‑1 (MQ=255)                              aGGTCGAGCGCGTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGATAGAAAGATTCTGGTGATTCCGCGCCAAGGCGGGCAAAGcaccta                                                                                      <  2:67101/136‑3 (MQ=255)                                gTCGAGCGCGTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTg                                                                                    <  1:66420/136‑1 (MQ=255)                                gTCGAGCGCGTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTg                                                                                    <  1:397611/136‑1 (MQ=255)                                      cgcgTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAgg                                                                              >  2:232842/1‑136 (MQ=255)                                      cgcgTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAgg                                                                              >  2:236364/1‑136 (MQ=255)                                       gcgTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGgc                                                                             >  1:72638/1‑136 (MQ=255)                                       gcgTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGgc                                                                             >  1:381664/1‑136 (MQ=255)                                       gcgTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGgc                                                                             >  2:243336/1‑136 (MQ=255)                                       gcgTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGgc                                                                             >  1:299223/1‑136 (MQ=255)                                       gcgTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGgc                                                                             >  1:13675/1‑136 (MQ=255)                                         gTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGgcgc                                                                           <  1:93601/136‑1 (MQ=255)                                            gACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTgg                                                                        >  1:303975/1‑136 (MQ=255)                                            gACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTgg                                                                        <  2:348789/136‑1 (MQ=255)                                                    cAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGgtcgtc                                                                <  2:294610/136‑1 (MQ=255)                                                    cAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGgtcgtc                                                                <  2:161483/136‑1 (MQ=255)                                                    cAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGgtcgtc                                                                <  1:131772/136‑1 (MQ=255)                                                    cAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGgtcgtc                                                                <  1:379989/136‑1 (MQ=255)                                                          gTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGTCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCgg                                                                     <  1:35888/125‑1 (MQ=255)                                                          gTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGTCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCgg                                                                     >  2:257944/1‑125 (MQ=255)                                                          gTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGTCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCgg                                                                     <  1:257944/125‑1 (MQ=255)                                                          gTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGTCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCgg                                                                     >  2:35888/1‑125 (MQ=255)                                                          gTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATc                                                          <  1:303376/136‑1 (MQ=255)                                                               aTTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCATGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCTCCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATaa                                                     >  2:74612/1‑136 (MQ=255)                                                                            cACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTcc                                        >  1:71014/1‑136 (MQ=255)                                                                             aCCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTccc                                       <  1:369564/136‑1 (MQ=255)                                                                                              gTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTCCCGCTACCAGTTCCGCCTc                      >  2:384989/1‑136 (MQ=255)                                                                                              gTACGACTGCACCGCCAAGGAGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTCCGGCTACCAGTTCCGCCTc                      >  2:315543/1‑136 (MQ=255)                                                                                                  gACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAg                                                    <  2:111054/102‑1 (MQ=255)                                                                                                  gACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAg                                                    >  1:111054/1‑102 (MQ=255)                                                                                                  gACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTCCCGCTACCAGTTCCGCCTCGGTg                  >  2:303611/1‑136 (MQ=255)                                                                                                    cTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTCCCGCTACCAGTTCCGCCTCGGTGAt                >  1:416723/1‑136 (MQ=255)                                                                                                          cGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTCCCGCTACCAGTTCCGCCTCGGTGATCGCAAc          >  1:41084/1‑136 (MQ=255)                                                                                                               aGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTCCCGCTACCAGTTCCGCCTCGGTGATCGCAACTCGCa     >  2:246536/1‑136 (MQ=255)                                                                                                                  tGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTCCCGCTACCAGTTCCGCCTCGGTGATCGCAACTCGCAGaa  <  2:293655/136‑1 (MQ=255)                                                                                                                    |                                                                                                                                      CTGGGAGGTGTCTTGCCCGGTCGACACCAAGGTCGAGCGCGTCGACCAGCGCAATGAGTCGGATTACAACTTTGTCACCCGCCTGGCCAAGCAGTACGACTGCACCGCCAAGGTGGCCGAAAGCAAGCTGCTGGTGATGCCGCGCCAAGGCGGGCAAAGCACCACTGGCAAGGCGCTGGCGGTCGTCACCATCAATAAGACGTACGTTTCCCGCTACCAGTTCCGCCTCGGTGATCGCAACTCGCAGAA  >  NC_002947/3450979‑3451227

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