Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,424,895	C→G	9.1%	intergenic (+483/‑29)	PP_3039 → / → PP_3040	hypothetical protein/holin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,424,895	0	C	G	9.1%	120.9 / 4.3	44	intergenic (+483/‑29)	PP_3039/PP_3040	hypothetical protein/holin
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base G (2/2);  total (24/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 = 6.92e-01

CAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGCTG  >  NC_002947/3424748‑3425042                                                                                                                                                      |                                                                                                                                                         cAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGgcg                                                                                                                                                         >  1:38117/1‑149 (MQ=255)    aCATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCAGCCAAACACTGAAAACCCCGGCCCTGGCGCCg                                                                                                                                                      >  5:60423/1‑149 (MQ=255)          cccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGtt                                                                                                                                                 <  8:252072/148‑1 (MQ=255)             ggACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTg                                                                                                                                             >  3:142108/1‑149 (MQ=255)               aCGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGcc                                                                                                                                           <  1:112622/149‑1 (MQ=255)                     aCCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                                      <  4:142108/148‑1 (MQ=255)                       cACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCGGGCGCCGGGGttt                                                                                                                                                >  4:91200/1‑136 (MQ=255)                        aCTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGTAGGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgggg                                                                                                                                    >  1:218756/1‑149 (MQ=255)                        aCTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCTCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                                              <  8:250024/137‑1 (MQ=255)                        aCTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                                              >  7:250024/1‑137 (MQ=255)                                  gggTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAAc                                                                                                                        <  6:60423/149‑1 (MQ=255)                                         aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTg                                                                                                                                             >  2:203657/1‑121 (MQ=255)                                         aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTg                                                                                                                                             <  1:203657/121‑1 (MQ=255)                                         aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCttt                                                                                                                                        <  2:35387/126‑1 (MQ=255)                                         aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCttt                                                                                                                                        >  1:35387/1‑126 (MQ=255)                                            aaaaCTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCCCAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCgggg                                                                                                                                                   >  2:152615/1‑112 (MQ=255)                                            aaaaCTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCgggg                                                                                                                                                   <  1:152615/112‑1 (MQ=255)                                             aaaCTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGa‑‑‑‑‑gcagc                                                                                                         <  6:208046/148‑1 (MQ=255)                                              aaCTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                            <  5:184824/133‑1 (MQ=255)                                                   gCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGTCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCg                                                                                                    <  5:177785/147‑1 (MQ=255)                                                        ctgcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGtt                                                                                              >  6:70518/1‑148 (MQ=255)                                                           ccACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTgaga                                                                                          <  2:38117/149‑1 (MQ=255)                                                                          ggTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCtggttgg                                                                           >  2:89855/2‑149 (MQ=255)                                                                            tAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCTA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTg                                                                         <  2:218756/149‑1 (MQ=255)                                                                            tAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCtggttggt                                                                          <  5:94979/148‑1 (MQ=255)                                                                              aaaaGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGa                                                                        >  6:30220/1‑148 (MQ=255)                                                                                     cccATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCtggt                                                                              <  8:97834/135‑1 (MQ=255)                                                                                     cccATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCtggt                                                                              >  7:97834/1‑135 (MQ=255)                                                                                        aTCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTcc                                                             >  7:52794/1‑149 (MQ=255)                                                                                            tCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                                     <  8:69763/78‑1 (MQ=255)                                                                                            tCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                                     >  7:69763/1‑78 (MQ=255)                                                                                                            cgcACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGgtgt                                             >  6:23909/1‑147 (MQ=255)                                                                                                                 aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCga                                      <  5:23909/149‑1 (MQ=255)                                                                                                                        cccGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTgg                               >  2:260102/1‑149 (MQ=255)                                                                                                                          cGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCg                             >  5:277523/1‑149 (MQ=255)                                                                                                                             cAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGc                          <  5:70518/149‑1 (MQ=255)                                                                                                                              aaaCACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCtggttg                                                                            <  4:132425/98‑1 (MQ=255)                                                                                                                              aaaCACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCtggttg                                                                            >  3:132425/1‑98 (MQ=255)                                                                                                                                    tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCGTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACAt                    >  3:195578/1‑148 (MQ=255)                                                                                                                                    tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCgggg                                                >  2:2814/1‑120 (MQ=255)                                                                                                                                    tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATg                   >  7:93467/1‑149 (MQ=255)                                                                                                                                     aaaaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATgg                  >  5:47102/2‑149 (MQ=255)                                                                                                                                       aaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGc                 >  7:251184/1‑148 (MQ=255)                                                                                                                                       aaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGc                 <  8:251184/148‑1 (MQ=255)                                                                                                                                          ccccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGc             >  4:126028/1‑149 (MQ=255)                                                                                                                                                     gcgcCGGGGTTTTTGCCTTTGTGGGATTCGATTAACACCGACCTCCCCCGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTa         <  1:89855/147‑1 (MQ=255)                                                                                                                                                     gcgcCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGCTg  <  1:260102/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                         CAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCG‑A‑GACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGA‑‑‑‑‑GCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCGGCTG  >  NC_002947/3424748‑3425042

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