Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	58,439	G→C	8.1%	intergenic (‑350/+338)	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,439	0	G	C	8.1%	166.9 / 4.1	54	intergenic (‑350/+338)	PP_5424/PP_0050	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base C (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAAT  >  NC_002947/58291‑58585                                                                                                                                                      |                                                                                                                                                   cTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTg‑g                                                                                                                                                    >  7:123666/1‑149 (MQ=255)  tCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑Gt                                                                                                                                                   >  7:311407/1‑149 (MQ=255)      tGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCct                                                                                                                                               >  1:212944/1‑149 (MQ=255)          gccggccATCGTCATGAATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAG‑CTTCCTCTgc                                                                                                                                           <  4:142398/149‑2 (MQ=255)           ccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGaa                                                                                                                                          <  3:317005/149‑1 (MQ=255)              gccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAagta                                                                                                                                       <  8:311407/149‑1 (MQ=255)                cATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAagtagt                                                                                                                                     <  7:295198/149‑1 (MQ=255)                          cATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCggg                                                                                                                            >  4:151096/1‑148 (MQ=255)                          cATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCggg                                                                                                                            <  8:195591/148‑1 (MQ=255)                               gATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCGCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGc                                                                                                                      >  5:114232/1‑149 (MQ=255)                                      attgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCg                                                                                                                >  1:229912/1‑148 (MQ=255)                                      attgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCg                                                                                                                >  2:236920/1‑148 (MQ=255)                                          acacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGAtt                                                                                                            <  2:229912/148‑1 (MQ=255)                                          acacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTc                                                                                                           >  6:130396/1‑149 (MQ=255)                                           cacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCAtt                                                                                                                        <  5:243520/135‑1 (MQ=255)                                           cacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCAtt                                                                                                                        >  6:243520/1‑135 (MQ=255)                                              gcgTGCCGCTTTGGATAGCCTTTCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCGTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                       >  7:141683/1‑149 (MQ=255)                                              gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                       <  8:123666/149‑1 (MQ=255)                                                   ccGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAcccc                                                                                                    >  2:220566/1‑147 (MQ=255)                                                    cGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCaa                                                                                                  <  4:180955/148‑1 (MQ=255)                                                     gCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAgg                                                                                                >  7:108500/1‑149 (MQ=255)                                                      cTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGt                                                                                               >  8:304299/1‑149 (MQ=255)                                                          ggATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGc                                                                                            <  5:69087/148‑1 (MQ=255)                                                           gATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGATCCCAAGGTCGc                                                                                            >  8:225987/1‑147 (MQ=255)                                                           gATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCa                                                                                          <  3:151096/149‑1 (MQ=255)                                                           gATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCa                                                                                          <  8:194097/149‑1 (MQ=255)                                                               gCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGAc                                                                                       <  5:130396/148‑1 (MQ=255)                                                                 cTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAg                                                                                    >  2:293609/1‑149 (MQ=255)                                                                  ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                       <  3:161204/129‑1 (MQ=255)                                                                  ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                       >  4:161204/1‑129 (MQ=255)                                                                  ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAg                                                                                    >  5:313725/1‑148 (MQ=255)                                                                  ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAg                                                                                    >  3:152340/1‑148 (MQ=255)                                                                    gCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAgg                                                                                   <  6:114232/147‑1 (MQ=255)                                                                    gCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAgg                                                                                   <  7:75966/147‑1 (MQ=255)                                                                          gAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTaaa                                                                            <  4:155076/148‑1 (MQ=255)                                                                           aGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCGATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAg                                                                           >  6:46182/1‑148 (MQ=255)                                                                                  cAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGc                                                                     <  6:313725/147‑1 (MQ=255)                                                                                        cgcATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGcc                                                                                                                     <  8:19604/93‑1 (MQ=255)                                                                                         gcATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGactac                                                            >  1:306700/1‑149 (MQ=255)                                                                                           aTTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGAATTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGactact                                                           >  1:97348/1‑149 (MQ=255)                                                                                           aTTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActc                                                          >  2:68112/1‑149 (MQ=255)                                                                                                tCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCAt                                                     >  2:281985/1‑149 (MQ=255)                                                                                                  gggCATGGGTACTAGCATGAGCCCCGGCTATTTTGGGAGCTGAGGAAGCTG‑GTTCCTCTGAAGTAGTGAAGGGGGGGATTGCCTGTCGGATTTTGATCACTAGGTCGACAGAGTAGGAGCTAAAGGATCGATGCGAGTACTCCTGAtgc                                                   <  8:141683/149‑3 (MQ=255)                                                                                                        gggTACCAGCATGAGCCCCCGCGATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTAcc                                             <  5:46182/149‑1 (MQ=255)                                                                                                            aCCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGaa                                          >  2:227174/1‑148 (MQ=255)                                                                                                              cAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAg                                       <  2:306700/149‑1 (MQ=255)                                                                                                                  ataaGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                                   <  2:197728/146‑1 (MQ=255)                                                                                                                       cccccGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtc                              <  8:108500/149‑1 (MQ=255)                                                                                                                         cccGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcgg                            >  8:9733/1‑149 (MQ=255)                                                                                                                                       ttcagAGGAACCAG‑CTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTc                                     <  6:33795/124‑1 (MQ=255)                                                                                                                                       ttcagAGGAACCAG‑CTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTc                                     >  5:33795/3‑126 (MQ=255)                                                                                                                                       ttcagAGGAACCAG‑CTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCt              >  8:9015/3‑149 (MQ=255)                                                                                                                                           gagGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc           >  5:277426/1‑148 (MQ=255)                                                                                                                                             ggAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGTCCTTTCCAg        <  1:236920/149‑1 (MQ=255)                                                                                                                                                   tG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAAt  >  7:62660/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTG‑GTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAAT  >  NC_002947/58291‑58585

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