Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,006,479	C→A	13.4%	D1020E (GAC→GAA)	PP_2627 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,006,479	0	C	A	13.4%	120.4 / 6.5	43	D1020E (GAC→GAA)	PP_2627	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base A (4/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAG  >  NC_002947/3006343‑3006620                                                                                                                                         |                                                                                                                                                  aGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGcc                                                                                                                                       >  6:335360/1‑149 (MQ=255)  gCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCt                                                                                                                                      <  7:81329/149‑1 (MQ=255)      cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTg                                                                                                                                    >  2:71360/1‑147 (MQ=255)      cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCa                                                                                                                                  <  2:231447/149‑1 (MQ=255)      cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCa                                                                                                                                  <  5:210242/149‑1 (MQ=255)            gCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCg                                                                                                                             <  3:87054/148‑1 (MQ=255)                   gCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCAc                                                                                                                     >  3:270696/1‑149 (MQ=255)                       ccccTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACg‑g                                                                                                                  <  5:6435/147‑1 (MQ=255)                          cTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCt                                                                                                                                      <  6:185644/125‑1 (MQ=255)                          cTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCt                                                                                                                                      >  5:185644/1‑125 (MQ=255)                              gtggtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATcc‑‑c                                                                                                        <  4:270696/148‑1 (MQ=255)                                          cccGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGa                                                                                           >  8:8392/1‑149 (MQ=255)                                           ccGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGAc                                                                                          >  4:148559/1‑149 (MQ=255)                                           ccGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGAc                                                                                          <  5:335360/149‑1 (MQ=255)                                             gATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGt                                                                                        >  7:56411/1‑149 (MQ=255)                                               tGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGa                                                                                      <  8:56411/149‑1 (MQ=255)                                                         gATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCt                                                                            >  5:63794/1‑149 (MQ=255)                                                           tGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCg                                                                          <  4:167063/149‑1 (MQ=255)                                                                  gAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGCGGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTgg                                                                    <  1:71360/149‑1 (MQ=255)                                                                       gTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAg‑c                                                             <  3:148559/149‑1 (MQ=255)                                                                                  cGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTTGAGCTGGCGGAG‑CGCATCCGCGa                                                   >  6:25318/1‑148 (MQ=255)                                                                                   gAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGTGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGaaa                                                 >  7:285463/1‑149 (MQ=255)                                                                                         gcgcCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGTGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACcttct                                           <  8:285463/149‑1 (MQ=255)                                                                                               gcagcaGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGAt                                      >  7:104175/1‑148 (MQ=255)                                                                                                     gTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTg                               <  7:8392/149‑1 (MQ=255)                                                                                                      tttAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAg‑c                                                             >  8:270513/1‑118 (MQ=255)                                                                                                      tttAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAg‑c                                                             <  7:270513/118‑1 (MQ=255)                                                                                                           tGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGc                         <  8:50280/149‑1 (MQ=255)                                                                                                                 gTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCCGACTGATTCGTGCCGAGGTGTGCGGCCGGCAGGAGAT‑GCGGCGCCGCTTCCGCGAACCCTTCTTCGATCACGGTGGAATCCTAAg                     >  1:276040/1‑149 (MQ=255)                                                                                                                    tctcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg                 <  8:104175/148‑1 (MQ=255)                                                                                                                     ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGa                                                                         <  8:141181/92‑1 (MQ=38)                                                                                                                     ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGa                                                                         >  7:141181/1‑92 (MQ=38)                                                                                                                     ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg                 >  6:83343/1‑147 (MQ=255)                                                                                                                     ctcAAGGTTGTCCTTCAGGAAAACCTTGAGACCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGAGCCGACGAGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCCCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg                 >  7:256601/1‑147 (MQ=255)                                                                                                                      tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg              >  1:133696/1‑149 (MQ=255)                                                                                                                       cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCg                                                                          <  7:56179/89‑1 (MQ=255)                                                                                                                       cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCg                                                                          >  8:56179/1‑89 (MQ=255)                                                                                                                       cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt             >  4:99122/1‑149 (MQ=255)                                                                                                                       cAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGTCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt             >  6:246671/1‑149 (MQ=255)                                                                                                                          ggTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGtt            >  5:7966/1‑147 (MQ=255)                                                                                                                            ttttCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGa                                                                                  >  8:30641/1‑76 (MQ=255)                                                                                                                               tccttcaGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCAt      <  5:246671/142‑1 (MQ=255)                                                                                                                                  tGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAg  >  8:192999/1‑149 (MQ=255)                                                                                                                                         |                                                                                                                                                  AGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAG  >  NC_002947/3006343‑3006620

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