Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,891,803	A→T	8.0%	L13Q (CTG→CAG)	PP_2546 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,891,803	0	A	T	8.0%	139.1 / 2.7	50	L13Q (CTG→CAG)	PP_2546	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (24/22);  new base T (2/2);  total (26/24)
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.46e-01

TTCCCAGGTGGTGAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGGCGC  >  NC_002947/2891658‑2891946                                                                                                                                                  |                                                                                                                                                ttCCCAGGTGGTGAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTgcagc                                                                                                                                               <  3:58933/148‑1 (MQ=255) ttCCCAGGTGGTGAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTgcagc                                                                                                                                               <  5:63867/148‑1 (MQ=255)  tCCCAGGTGGTGAGCGGGCAGCCCAGGGGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACGGGCAAGTGCATCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTgcagca                                                                                                                                              <  7:135319/148‑1 (MQ=255)         tggtgAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAgg                                                                                                                                      <  8:10135/149‑1 (MQ=255)                   gCAGCCCAGGTGCAGGGACTCCACCGCCAGGCCCCAGGCCAGGACCGGCATGTGCAGCAAGAGGGCAGGGGGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTcggcgg                                                                                                                            <  8:352276/149‑1 (MQ=255)                   gCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTcggcgg                                                                                                                            >  3:26951/1‑149 (MQ=255)                           ggTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCgg                                                                                                                    <  7:180932/149‑1 (MQ=255)                            gTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGccg                                                                                                                    >  2:257477/1‑146 (MQ=255)                                      tCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACGAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGAGTCGGCGGCAACGCGGTAGAGCAtgcg                                                                                                         >  1:99705/1‑147 (MQ=255)                                        cACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCAtgggtt                                                                                                       <  8:109999/149‑1 (MQ=255)                                        cACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCAtgggtt                                                                                                       <  8:104819/149‑1 (MQ=255)                                          ccGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCAtgggttgg                                                                                                     <  5:379326/149‑1 (MQ=255)                                            gCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCAtgggttgggt                                                                                                   >  2:59849/1‑149 (MQ=255)                                                       ggccaggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCgg                                                                                             >  6:208957/1‑144 (MQ=255)                                                       ggccaggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCgg                                                                                             <  5:208957/144‑1 (MQ=255)                                                       ggccaggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGAtt                                                                                        >  1:168607/1‑149 (MQ=255)                                                           aggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTgggg                                                                                    <  1:257477/149‑1 (MQ=255)                                                             gACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGtt                                                                                  >  7:104826/1‑149 (MQ=255)                                                                         cagcaACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAg                                                                      <  4:26951/149‑1 (MQ=255)                                                                           gcaACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGt                                                                                   >  3:312530/1‑134 (MQ=255)                                                                           gcaACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGt                                                                                   <  4:312530/134‑1 (MQ=255)                                                                             aaCAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTgg                                                                  <  1:155586/149‑1 (MQ=255)                                                                             aaCAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTgg                                                                  >  7:214098/1‑149 (MQ=255)                                                                                 gggcagggcGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTg                                                              >  7:45343/1‑149 (MQ=255)                                                                                    cagggcGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCCGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCt                                                           >  3:208590/1‑149 (MQ=255)                                                                                    cagggcGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCt                                                           >  6:8043/1‑149 (MQ=255)                                                                                     agggcGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACAAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGCTTATCCCAGCATAGTTGGGTTGCCt                                                           >  1:304607/1‑148 (MQ=255)                                                                                          gcCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTg                                                     >  8:146866/1‑149 (MQ=255)                                                                                           ccAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTgg                                                    <  8:45343/149‑1 (MQ=255)                                                                                                gcagcaCCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCCCCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGAGTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTa                                                >  3:216767/1‑148 (MQ=255)                                                                                                gcagcaCCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTa                                                >  3:94875/1‑148 (MQ=255)                                                                                                  agcaCCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTacac                                             >  6:134770/1‑149 (MQ=255)                                                                                                     aCCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGATGGGTTGCCTGATCTGGCCTACAcccc                                          >  2:279393/1‑149 (MQ=255)                                                                                                        agcagcCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTc                                       >  1:154614/1‑149 (MQ=255)                                                                                                         gcagcCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTcc                                      <  8:214098/149‑1 (MQ=255)                                                                                                          cagcCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCt                                     <  8:104826/149‑1 (MQ=255)                                                                                                             cccGCCGTACATCAGCAGCAGGATGAAGGCCAGGTGGAGCTGGACCAGGGTGTCGGCGGCAAGGGGGTAGTGGATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTTATCTGGGCTACACCCCTTCCTCGa                                  <  2:99705/149‑1 (MQ=255)                                                                                                             cccGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGa                                  <  7:146866/149‑1 (MQ=255)                                                                                                                  cgAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATc                             <  1:59849/149‑1 (MQ=255)                                                                                                                     aCAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGc                          <  2:168607/149‑1 (MQ=255)                                                                                                                     aCAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGc                          <  4:208590/149‑1 (MQ=255)                                                                                                                                  gATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGcca             >  1:205820/1‑149 (MQ=255)                                                                                                                                    tGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGccacc           >  6:109823/1‑149 (MQ=255)                                                                                                                                     gAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCTGTAGAGCATGGGTTGGGTGCCCGGGGATTGGCGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCGTACCTCGACCATCAGCAGGAAAGTCGcaaaca          >  5:107933/1‑144 (MQ=255)                                                                                                                                       aGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTc                                       <  4:219430/118‑1 (MQ=255)                                                                                                                                       aGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTc                                       >  3:219430/1‑118 (MQ=255)                                                                                                                                          cctggTGCTGCTGCACCTGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCAt                              >  7:104587/4‑124 (MQ=255)                                                                                                                                          cctggTGCTGCTGCACCTGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCAt                              <  8:104587/121‑1 (MQ=255)                                                                                                                                           ctggTGCTGCTGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGcc              <  6:54116/137‑1 (MQ=255)                                                                                                                                           ctggTGCTGCTGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGcc              >  5:54116/3‑139 (MQ=255)                                                                                                                                             ggTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTATAGCATGGGTTGGGTGCCCGTGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGgcgc  <  8:301642/149‑1 (MQ=255)                                                                                                                                              gTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGccacca          <  6:134425/140‑1 (MQ=255)                                                                                                                                              gTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGccacca          >  5:134425/1‑140 (MQ=255)                                                                                                                                                  |                                                                                                                                                TTCCCAGGTGGTGAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGGCGC  >  NC_002947/2891658‑2891946

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