Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	735,621	C→T	11.1%	intergenic (‑44/+89)	PP_0627 ← / ← PP_0628	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	735,621	0	C	T	11.1%	115.6 / 6.3	45	intergenic (‑44/+89)	PP_0627/PP_0628	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base T (2/3);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

AAGGAGCAGCCACGGCCGGAAGGGCCTGCGCTCGACTTGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGAAGACGATGCTGCGGCCCAGGTTGCCGCGCAG  >  NC_002947/735473‑735763                                                                                                                                                     |                                                                                                                                               aaGGAGCAGCCACGGCCGGAAGGGCCTGCGCTCGACTTGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACa                                                                                                                                               >  1:581943/1‑150 (MQ=255)              ggCCGGAAGGGCCTGCGCTCGACTTGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCGAACGCTGTAGCGTTTGCGCAgt                                                                                                                                  <  2:474620/150‑2 (MQ=255)                 cGGAAGGGCCTGCGCTCGACTTGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGgcgc                                                                                                                               <  2:455462/150‑1 (MQ=255)                   gAAGGGCCTGCGCTCGACTTGGTGCTGGGGGGCTCAGCCTGACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCtt                                                                                                                             <  1:640138/150‑1 (MQ=255)                        gCCTGCGCTCGACTTGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAAc                                                                                                                        >  2:638003/1‑150 (MQ=255)                          cTGCGCTCGACTTGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACaa                                                                                                                      >  2:475121/1‑150 (MQ=255)                           tGCGCTCGACTTGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACaaa                                                                                                                     <  1:128417/150‑1 (MQ=255)                                   aCTTGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACaaa                                                                                                                     <  2:715991/142‑1 (MQ=255)                                   aCTTGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACaaa                                                                                                                     >  1:715991/1‑142 (MQ=255)                                     ttGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCtt                                                                                                                             >  1:184032/1‑132 (MQ=255)                                     ttGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCtt                                                                                                                             <  2:184032/132‑1 (MQ=255)                                        gTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTc                                                                                                        <  2:570918/150‑1 (MQ=255)                                             ggggggCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATcc                                                                                                   >  1:604930/1‑150 (MQ=255)                                                    tCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACaaa                                                                                                                     >  2:643261/1‑125 (MQ=255)                                                    tCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACaaa                                                                                                                     <  1:643261/125‑1 (MQ=255)                                                    tCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGgcgc                                                                                            >  1:424212/1‑150 (MQ=255)                                                        aGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGCCTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAg                                                                                        >  1:347732/1‑150 (MQ=255)                                                              tCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACg                                                                                  >  1:458556/1‑150 (MQ=255)                                                                aTCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGcc                                                                                >  1:122610/1‑150 (MQ=255)                                                                 tCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCa                                                                               >  1:712291/1‑150 (MQ=255)                                                                  cGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCAt                                                                              <  2:424212/150‑1 (MQ=255)                                                                              gATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCg                                                                  <  2:604930/150‑1 (MQ=255)                                                                                  gtCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCGATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGc                                                              >  2:546052/2‑150 (MQ=255)                                                                                    ccGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCAc                                                            <  1:638003/150‑1 (MQ=255)                                                                                       ggTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTgg                                                         <  2:581943/150‑1 (MQ=255)                                                                                         ttCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAGATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGCCGCTGCACTGGAt                                                       >  2:280852/1‑150 (MQ=255)                                                                                             ggCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCaaa                                                   <  1:616417/150‑1 (MQ=255)                                                                                               cGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGa                                                 <  2:472830/150‑1 (MQ=255)                                                                                                        aTGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCt                                        <  1:546052/150‑1 (MQ=255)                                                                                                          ggtggGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTAc                                      <  2:347732/150‑1 (MQ=255)                                                                                                               gCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGa                                 >  1:95191/1‑150 (MQ=255)                                                                                                                  tCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGaa                                >  1:642103/1‑148 (MQ=255)                                                                                                                  tCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGaa                                <  2:642103/148‑1 (MQ=255)                                                                                                                        ggTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGAAGACGAtgc                        >  1:328392/1‑150 (MQ=255)                                                                                                                          tAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGAAGACGAtgctg                      >  2:459536/1‑150 (MQ=255)                                                                                                                                tttGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGAAGACGATGCTGCGGccc                >  2:262776/1‑150 (MQ=255)                                                                                                                                    tgACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGc                                                                                 <  1:582207/81‑1 (MQ=255)                                                                                                                                    tgACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGc                                                                                 >  2:582207/1‑81 (MQ=255)                                                                                                                                      aCTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGAAGACGATGCTGCGGCCCAGGTTg          >  2:319441/1‑150 (MQ=255)                                                                                                                                       cTGCGCGAACGCTGTAGCGTTTGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGa                                             <  2:436921/114‑1 (MQ=39)                                                                                                                                       cTGCGCGAACGCTGTAGCGTTTGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGa                                             >  1:436921/1‑114 (MQ=39)                                                                                                                                       cTGCGCGAACGCTGTAGCGTTTGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCg                                  >  1:76334/1‑125 (MQ=255)                                                                                                                                       cTGCGCGAACGCTGTAGCGTTTGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCg                                  <  2:76334/125‑1 (MQ=255)                                                                                                                                             aaaCGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGAAGACGATGCTGCGGCCCAGGTTGCCGCGCa   <  1:319441/150‑1 (MQ=255)                                                                                                                                              aaCGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGAAGACGATGCTGCGGCCCAGGTTGCCGCGCAg  <  1:262776/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               AAGGAGCAGCCACGGCCGGAAGGGCCTGCGCTCGACTTGGTGCTGGGGGGCTCGCAGATACTCATCGACACGACGCTGATCTTCCGGGTTCAGGCGACTGGTCATGGTGGGCCTCGTCAGGTAGACGTTTGTGACTGCGCAAACGCTACAGCGTTCGCGCAGGCGCTTGAAACAAATGATAATGCTTTCTATCCCTGGCGCCGAGTGTACGCCATCGCAAACAGTCGCTGCACTGGATCAAAGACTGATGCCTACGTCGAAGACGATGCTGCGGCCCAGGTTGCCGCGCAG  >  NC_002947/735473‑735763

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