Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,549,653	T→A	8.7%	intergenic (+18/‑154)	PP_t47 → / → PP_4881	tRNA‑Leu/iron ABC transporter substrate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,549,653	0	T	A	8.7%	130.6 / 4.7	46	intergenic (+18/‑154)	PP_t47/PP_4881	tRNA‑Leu/iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base A (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

CGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGGAA  >  NC_002947/5549511‑5549798                                                                                                                                               |                                                                                                                                                  cGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGt                                                                                                                                            <  1:140402/150‑1 (MQ=255)   aaCGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTgg                                                                                                                                          >  2:604743/1‑150 (MQ=255)       aTAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTa                                                                                                                                              >  2:195351/1‑142 (MQ=255)       aTAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTa                                                                                                                                              <  1:195351/142‑1 (MQ=255)          gTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGtttttt                                                                                                                                   >  2:50768/1‑150 (MQ=255)                            tgTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTAtgt                                                                                                                 >  1:303599/1‑150 (MQ=255)                            tgTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTAt                                                                                                                   >  2:766242/1‑148 (MQ=255)                            tgTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTAt                                                                                                                   <  1:766242/148‑1 (MQ=255)                              tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGc                                                                                                                                                <  2:232250/117‑1 (MQ=255)                              tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTAAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGc                                                                                                                                                >  1:232250/1‑117 (MQ=255)                                   cTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGa                                                                                                                        >  1:313816/1‑136 (MQ=255)                                   cTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGa                                                                                                                        <  2:313816/136‑1 (MQ=255)                                    tttGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTgg                                                                                                                                          >  1:324548/1‑117 (MQ=255)                                    tttGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTgg                                                                                                                                          <  2:324548/117‑1 (MQ=255)                                        cccAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTccccc                                                                                                     >  2:756509/1‑150 (MQ=255)                                          cAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTccccccc                                                                                                   <  1:610192/150‑1 (MQ=255)                                           aGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCt                                                                                                  <  1:50768/150‑1 (MQ=255)                                                        ggTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCt                                                                                     <  1:56974/150‑1 (MQ=255)                                                                  cTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAAt                                                                           <  1:604743/150‑1 (MQ=255)                                                                  cTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAAt                                                                           >  1:227736/1‑150 (MQ=255)                                                                      tttCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCAt                                                                               <  2:267/142‑1 (MQ=255)                                                                      tttCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCAt                                                                               >  1:267/1‑142 (MQ=255)                                                                       ttCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTtctc                                                                      >  2:552464/1‑150 (MQ=255)                                                                          aGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGc                                                                                                               <  1:682537/106‑1 (MQ=255)                                                                          aGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGc                                                                                                               >  2:682537/1‑106 (MQ=255)                                                                                  gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGtt                                                                                                                                       >  2:179057/1‑74 (MQ=255)                                                                                  gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGtt                                                                                                                                       <  1:179057/74‑1 (MQ=255)                                                                                      ccGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACtt                                                       <  2:786238/150‑1 (MQ=255)                                                                                      ccGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACtt                                                       <  1:756509/150‑1 (MQ=255)                                                                                          aaGTTCTTGAAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTATCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTgaga                                                   <  1:532338/150‑1 (MQ=255)                                                                                                       ttCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGttt                                      >  1:510616/1‑150 (MQ=255)                                                                                                             tcttctCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGc                                                                                                               >  1:361241/1‑71 (MQ=255)                                                                                                             tcttctCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGc                                                                                                               <  2:361241/71‑1 (MQ=255)                                                                                                             tcttctCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATg                                                                                         <  1:47246/93‑1 (MQ=255)                                                                                                             tcttctCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATg                                                                                         >  2:47246/1‑93 (MQ=255)                                                                                                                           ccATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGtt                                       <  1:352334/129‑1 (MQ=255)                                                                                                                           ccATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGtt                                       >  2:352334/1‑129 (MQ=255)                                                                                                                            cATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAAAAACAATTTCGTTTACCATTGGTTCCAAATATGTa                 >  1:701633/1‑150 (MQ=255)                                                                                                                                   aaaCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTaa          >  1:399646/1‑150 (MQ=255)                                                                                                                                      cccACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTtctc                                                                      <  2:690149/87‑1 (MQ=255)                                                                                                                                      cccACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTtctc                                                                      >  1:690149/1‑87 (MQ=255)                                                                                                                                      cccACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCg       >  1:489719/1‑150 (MQ=255)                                                                                                                                       ccACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGa      <  2:489719/150‑1 (MQ=255)                                                                                                                                         aCTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCTCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAgg    <  1:552464/150‑1 (MQ=255)                                                                                                                                         aCTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAgg    >  2:440613/1‑150 (MQ=255)                                                                                                                                           tAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCGCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGGaa  <  2:303599/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                  CGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGGAA  >  NC_002947/5549511‑5549798

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