Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	58,435	G→C	11.9%	intergenic (‑346/+342)	PP_5424 ← / ← PP_0050	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	58,435	0	G	C	11.9%	113.6 / 8.0	42	intergenic (‑346/+342)	PP_5424/PP_0050	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base C (3/2);  total (22/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 = 3.96e-01

ACCTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGC  >  NC_002947/58289‑58581                                                                                                                                                   |                                                                                                                                                   aCCTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGc                                                                                                                                                   <  2:9120/148‑1 (MQ=255)            gccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgc                                                                                                                                       <  7:225421/149‑2 (MQ=255)              cggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgctg                                                                                                                                     <  2:212214/149‑4 (MQ=255)                          tGCATCAGATCTTGATTGACTCGCGGGCCGCTTTGGATAGCCTTGCTATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACGAGCATGAGCCCGCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcggcgg                                                                                                                         <  1:308088/149‑1 (MQ=255)                           gCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCggg                                                                                                                        >  4:178873/1‑149 (MQ=255)                           gCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCggg                                                                                                                        >  6:45325/1‑149 (MQ=255)                           gCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCggg                                                                                                                        >  6:289514/1‑149 (MQ=255)                            cATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCggg                                                                                                                        >  5:204616/1‑148 (MQ=255)                             aTCAGATCTTGATTGACACGCGTGCCGCATTGCATAGCCTTGCAATTGAGGGCTGCAACTACGCATTTCTCGGGCCAGGGTACCAGCATGAGCACCCGCAATATTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcgtccgcca                                                                                                                      >  2:186277/1‑142 (MQ=255)                                aGATCTTGATTGACACGCGTGCCGCGTTGGATAGCCTTGCAATTGAGGTCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTg                                                                                                                   <  5:139524/149‑1 (MQ=255)                                     ttgattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGt                                                                                                              >  3:47840/1‑149 (MQ=255)                                     ttgattgaCACGCGTGACGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGt                                                                                                              >  2:14967/1‑149 (MQ=255)                                       gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTc                                                                                                             <  1:41323/148‑1 (MQ=255)                                       gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTc                                                                                                             <  3:178873/148‑1 (MQ=255)                                                     ccGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCaa                                                                                              <  4:47840/149‑1 (MQ=255)                                                     ccGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTACGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGTGGGCATTGCCAGTCGGATTCTGAACCCaa                                                                                              <  1:186277/149‑1 (MQ=255)                                                                  ccTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGAc                                                                                   >  3:129070/1‑147 (MQ=255)                                                                  ccTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTa                                                                                 <  5:289514/149‑1 (MQ=255)                                                                               ggCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGAt                                                                    <  4:95891/149‑1 (MQ=255)                                                                                     aaCTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGc                                                              <  1:14967/149‑1 (MQ=255)                                                                                        tGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGac                                                           >  8:34941/1‑149 (MQ=255)                                                                                            cATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGactact                                                       >  8:162529/1‑149 (MQ=255)                                                                                                  tCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGAc                                                                                   >  6:242173/1‑115 (MQ=255)                                                                                                  tCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGAc                                                                                   <  5:242173/115‑1 (MQ=255)                                                                                                   cGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGAGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATc                                                >  4:14441/1‑149 (MQ=255)                                                                                                        cTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTa                                           <  5:45325/149‑1 (MQ=255)                                                                                                              aCCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGac                                                           <  2:75367/127‑1 (MQ=255)                                                                                                              aCCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGac                                                           >  1:75367/1‑127 (MQ=255)                                                                                                                  gCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGt                                  >  2:40764/1‑148 (MQ=255)                                                                                                                    aTGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                               <  3:14441/149‑1 (MQ=255)                                                                                                                     tGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAcg                              >  3:107235/1‑149 (MQ=255)                                                                                                                        gCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAcgcg                            <  4:129070/148‑1 (MQ=255)                                                                                                                            ccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcggt                       >  8:169507/1‑149 (MQ=255)                                                                                                                                  ttttGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTaaa                 >  3:184655/1‑149 (MQ=255)                                                                                                                                      tacttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGc             >  6:14233/6‑149 (MQ=255)                                                                                                                                      ggCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACgg              <  4:184655/148‑1 (MQ=255)                                                                                                                                           cagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCttt        >  7:45849/1‑149 (MQ=255)                                                                                                                                             gagGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc       >  7:178605/1‑148 (MQ=255)                                                                                                                                                  aGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAg                                                                                             <  2:226624/57‑1 (MQ=255)                                                                                                                                                  aGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAg                                                                                             >  1:226624/1‑57 (MQ=255)                                                                                                                                                  aGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGc  >  2:141055/1‑148 (MQ=255)                                                                                                                                                  aGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGc  <  1:141055/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   ACCTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGC  >  NC_002947/58289‑58581

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