Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,646,266	A→T	7.0%	D540V (GAC→GTC)	PP_4954 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,646,266	0	A	T	7.0%	195.7 / 3.1	68	D540V (GAC→GTC)	PP_4954	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base T (2/3);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

CACCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGCTGAC  >  NC_002947/5646120‑5646414                                                                                                                                                    |                                                                                                                                                     cacCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACt                                                                                                                                                    <  7:26779/149‑1 (MQ=255)   cccTGGCCACTGCCGCGATCACCCTGATGTTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTgg                                                                                                                                                  <  7:106961/149‑1 (MQ=255)    ccTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCATGACTgg                                                                                                                                                  >  6:177923/1‑148 (MQ=255)    ccTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTACTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTgg                                                                                                                                                  >  2:127498/1‑148 (MQ=255)        gCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAggg                                                                                                                                             <  1:234138/149‑1 (MQ=255)        gCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAg                                                                                                                                               <  7:60443/147‑1 (MQ=255)         ccACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCAGTcaggcaggaa                                                                                                                                            <  5:56346/149‑3 (MQ=255)             tGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGAc                                                                                                                                        >  4:228522/1‑149 (MQ=255)               ccGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCAGTcaggcaggaacagga                                                                                                                                       <  1:221492/148‑8 (MQ=255)                cgcgATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGc                                                                                                                                     <  7:261698/149‑1 (MQ=255)                     tCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAAc                                                                                                                                <  2:19334/149‑1 (MQ=255)                        cccTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAg                                                                                                                             >  6:281912/1‑149 (MQ=255)                        cccTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAg                                                                                                                             <  5:281912/149‑1 (MQ=255)                         ccTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCTGGGCCGACGGCTGAACAAgg                                                                                                                            >  4:140823/1‑149 (MQ=255)                         ccTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTCCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAgg                                                                                                                            >  8:119842/1‑149 (MQ=255)                                        cTGCTTCAACCAGAACGGCGATGGCTACCGCCTGTTCCTGCCGCGCCTGTTCGATACACT‑GGTCGGCAGCCTGATTGCCAATCTAGCAGTGTTACTCTTCCAGCCTGACTGGCAAGGACGACGGCTGAACAAGTCGCTCGCCAAAAca                                                                                                              >  8:210100/1‑147 (MQ=255)                                            ttCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGc                                                                                                                                     <  3:205111/121‑1 (MQ=255)                                            ttCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGc                                                                                                                                     >  4:205111/1‑121 (MQ=255)                                             tCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGaa                                                                                                                                 <  6:169774/124‑1 (MQ=255)                                             tCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGaa                                                                                                                                 >  5:169774/1‑124 (MQ=255)                                              cAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGc                                                                                                         <  1:237818/147‑1 (MQ=255)                                                    gATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATa                                                                                                                >  6:227401/1‑134 (MQ=255)                                                    gATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATa                                                                                                                <  5:227401/134‑1 (MQ=255)                                                      tCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGc                                                                                               <  4:269880/149‑1 (MQ=255)                                                      tCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTCCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGc                                                                                               <  7:119842/149‑1 (MQ=255)                                                         gCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCgtg                                                                                            <  8:88660/149‑1 (MQ=255)                                                                cTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTgc                                                                                     >  6:282788/1‑149 (MQ=255)                                                                cTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTgc                                                                                     <  5:56023/149‑1 (MQ=255)                                                                   cGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGcc                                                                                  <  1:228097/149‑1 (MQ=255)                                                                    ggCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTgcgc                                                                                   >  1:60292/1‑147 (MQ=255)                                                                        tGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCACCAGCGTGTACCTGCGCCAGAt                                                                              >  6:261887/1‑148 (MQ=255)                                                                          ttCCTGCCGCGCCTGTTCGATACCCTGGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCa                                                                            >  1:3864/1‑149 (MQ=255)                                                                          ttCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTg                                                                                                           >  1:51100/1‑117 (MQ=255)                                                                          ttCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTg                                                                                                           <  2:51100/117‑1 (MQ=255)                                                                            ccTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                          <  2:3864/148‑1 (MQ=255)                                                                            ccTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTGGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                          <  3:228522/148‑1 (MQ=255)                                                                                ccGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATgcagca                                                                     <  1:127498/149‑1 (MQ=255)                                                                                 cgcgCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGcagcag                                                                    >  7:174011/1‑149 (MQ=255)                                                                                  gcgcCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGa                                                                               <  8:61622/137‑1 (MQ=255)                                                                                  gcgcCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGa                                                                               >  7:61622/1‑137 (MQ=255)                                                                                      cTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGc                                                               >  3:15416/1‑149 (MQ=255)                                                                                       tGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGcc                                                              >  4:13120/1‑149 (MQ=255)                                                                                         ttCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCa                                                            >  1:275412/1‑149 (MQ=255)                                                                                         ttCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCa                                                            >  4:125900/1‑149 (MQ=255)                                                                                           cGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACg                                                          >  3:101858/1‑149 (MQ=255)                                                                                            gATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACg                                                          <  5:282788/148‑1 (MQ=255)                                                                                                 ccT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTgc                                                                                     <  7:157692/116‑1 (MQ=255)                                                                                                 ccT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTgc                                                                                     >  8:157692/1‑116 (MQ=255)                                                                                                      gTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGcga                                                <  5:177923/149‑1 (MQ=255)                                                                                                          gCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTCACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCACCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCgacgac                                            <  5:261887/149‑1 (MQ=255)                                                                                                                 gATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGcc                                      >  6:4916/1‑148 (MQ=255)                                                                                                                      ccATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGCTCATGCAGCAGTATGCCCACGGCAAGCGCGACGCCCTGGCCTCCCGc                                >  8:120194/1‑149 (MQ=255)                                                                                                                         ttCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTg                             <  2:275412/149‑1 (MQ=255)                                                                                                                            ttGCAGTGTTCCTGTTACTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCAAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTCGCCTACCGCCAGGcc                          >  8:88490/1‑149 (MQ=255)                                                                                                                              ggaGTGGTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGccc                         <  4:173719/146‑1 (MQ=255)                                                                                                                                aGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgc                      <  3:125900/149‑1 (MQ=255)                                                                                                                                 gtgtTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgcg                     <  5:4916/149‑1 (MQ=255)                                                                                                                                 gtgtTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCGGACCAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGGACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgcg                     <  7:120194/149‑1 (MQ=255)                                                                                                                                       ctgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATg                <  6:189382/148‑1 (MQ=255)                                                                                                                                        tgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGcc              >  8:261151/1‑149 (MQ=255)                                                                                                                                          ttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCTcc                               <  8:223160/130‑1 (MQ=255)                                                                                                                                          ttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGcc                               >  7:223160/1‑130 (MQ=255)                                                                                                                                          ttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCca            <  3:13120/149‑1 (MQ=255)                                                                                                                                          gccctgccagTCAGGCAGTACCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCgt                                                                                             <  3:267720/66‑1 (MQ=255)                                                                                                                                           tcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATg                                                                >  2:164171/1‑96 (MQ=255)                                                                                                                                           ccctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCca            >  6:176062/2‑148 (MQ=255)                                                                                                                                             ctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACaa         >  7:114710/1‑149 (MQ=255)                                                                                                                                                 ctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGCt     <  4:101858/149‑1 (MQ=255)                                                                                                                                                    aCTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGCTGAc  <  2:54820/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CACCCTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCT‑GGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGCTGAC  >  NC_002947/5646120‑5646414

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