Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	436,905	A→T	7.0%	L128Q (CTG→CAG)	PP_0358 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	436,905	0	A	T	7.0%	169.9 / 3.1	57	L128Q (CTG→CAG)	PP_0358	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base T (2/2);  total (29/28)
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.76e-01

AGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGC  >  NC_002947/436761‑437047                                                                                                                                                 |                                                                                                                                               agcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAggcg                                                                                                                                            <  3:162555/149‑1 (MQ=255)                     cTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                         >  2:23107/1‑148 (MQ=255)                      ttACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggtg                                                                                                                                         >  4:320467/1‑129 (MQ=255)                      ttACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggtg                                                                                                                                         <  3:320467/131‑3 (MQ=255)                      ttACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGaa                                                                                                                        >  1:49484/1‑148 (MQ=255)                          ggATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAACCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                         >  6:333557/1‑143 (MQ=255)                          ggATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAACCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                         <  5:333557/143‑1 (MQ=255)                                gCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                         >  6:108494/1‑137 (MQ=255)                                gCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                         <  5:108494/137‑1 (MQ=255)                                  ccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                           >  5:1381/1‑149 (MQ=255)                                    aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                           >  5:96768/1‑147 (MQ=255)                                      aGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCCGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAAt                                                                                                                       >  5:67111/1‑133 (MQ=255)                                      aGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCCGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAAt                                                                                                                       <  6:67111/133‑1 (MQ=255)                                        aaGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                           >  5:340993/1‑143 (MQ=255)                                        aaGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                           <  6:340993/143‑1 (MQ=255)                                          gCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCTCACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACAt                                                                                                   <  6:1381/149‑1 (MQ=255)                                            cGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACAtt                                                                                                  >  3:36745/1‑148 (MQ=255)                                              aaGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACCCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGAt                                                                                               >  7:85949/1‑149 (MQ=255)                                                  tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGGAGGCGCTCACCAAGAATCa                                                                                                                     <  8:261664/123‑1 (MQ=255)                                                  tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                           >  4:207337/1‑149 (MQ=255)                                                     aCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggt                                                                                                                                          <  8:72675/99‑2 (MQ=255)                                                     aCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggt                                                                                                                                          >  7:72675/1‑98 (MQ=255)                                                          aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAt                                                                                   <  4:36745/149‑1 (MQ=255)                                                           gCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATg                                                                                  <  1:23107/149‑1 (MQ=255)                                                                 gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCg                                                                                          >  5:334023/1‑135 (MQ=255)                                                                 gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCg                                                                                          <  6:334023/135‑1 (MQ=255)                                                                            ggCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGc                                                                 <  3:207337/149‑1 (MQ=255)                                                                                 gCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGa                                                            <  5:339208/149‑1 (MQ=255)                                                                                 gCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGa                                                            <  2:49484/149‑1 (MQ=255)                                                                                  cTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGGAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGTCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTGGat                                                           <  8:19590/149‑2 (MQ=255)                                                                                      ccACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATa                                                       >  3:163722/1‑149 (MQ=255)                                                                                       cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAgg                                                                                                                                     >  7:23550/1‑70 (MQ=255)                                                                                       cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAgg                                                                                                                                     <  8:23550/70‑1 (MQ=255)                                                                                       cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGc                                                                 >  7:151567/1‑138 (MQ=255)                                                                                       cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGc                                                                 <  8:151567/138‑1 (MQ=255)                                                                                            gCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGc                                                   <  6:96768/147‑1 (MQ=255)                                                                                              cacaCCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAg                                                                                                                    >  4:114446/1‑80 (MQ=255)                                                                                              cacaCCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgtt                                               <  8:85949/149‑1 (MQ=255)                                                                                               acacCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgttg                                              >  5:64383/1‑149 (MQ=255)                                                                                                            cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGAc                                     >  4:164344/1‑145 (MQ=255)                                                                                                            cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGAc                                     <  3:164344/145‑1 (MQ=255)                                                                                                                 ccACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAc                                         >  1:159840/1‑136 (MQ=255)                                                                                                                 ccACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAc                                         <  2:159840/136‑1 (MQ=255)                                                                                                                 ccACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAg                            >  3:143512/1‑149 (MQ=255)                                                                                                                  cACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAGGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAgg                           >  5:268134/1‑149 (MQ=255)                                                                                                                         gCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGc                      >  1:131619/1‑147 (MQ=255)                                                                                                                         gCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGc                      <  2:131619/147‑1 (MQ=255)                                                                                                                               cccACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCTAGCATACGGCCGTTGTTCACAGACCACACGTAGGTCAGCCAGTTGcc              >  8:122502/1‑149 (MQ=255)                                                                                                                                 cACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCgg            >  6:38366/1‑149 (MQ=255)                                                                                                                                 cACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCgg            <  6:64383/149‑1 (MQ=255)                                                                                                                                    cGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGa         <  4:143512/149‑1 (MQ=255)                                                                                                                                      ccagccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAt        >  5:264694/1‑148 (MQ=255)                                                                                                                                       cagccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGTCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAt        >  4:275184/1‑147 (MQ=255)                                                                                                                                       cagccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACAGACCACACGTAGGTCAGCCAGTTGCCGGCGATCa      <  7:122502/149‑1 (MQ=255)                                                                                                                                        agccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAg     <  5:38366/149‑1 (MQ=255)                                                                                                                                           caTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGc  <  2:280563/149‑1 (MQ=255)                                                                                                                                           caTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGc  >  1:280563/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                               AGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGC  >  NC_002947/436761‑437047

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