Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,402,723	A→G	5.1%	L271P (CTG→CCG)	PP_1227 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,402,723	0	A	G	5.1%	231.7 / 2.0	80	L271P (CTG→CCG)	PP_1227	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base G (2/2);  total (39/41)
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.46e-01

CTCGGCCCATTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGAAGATTTCAGACAGCGCCTGAC  >  NC_002947/1402578‑1402870                                                                                                                                                     |                                                                                                                                                    cTCGGCCCATTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCt                                                                                                                                                  >  4:117857/1‑149 (MQ=255)  tCGGCCCATTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACGGCTc                                                                                                                                                 >  3:106982/1‑149 (MQ=255)  tCGGCCCATTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTc                                                                                                                                                 >  2:161727/1‑149 (MQ=255)    ggCCCATTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCGATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGa                                                                                                                                               <  4:219158/149‑1 (MQ=255)    ggCCCATTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGa                                                                                                                                               >  7:3468/1‑149 (MQ=255)    ggCCCATTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGa                                                                                                                                               <  1:42897/149‑1 (MQ=255)      cccATTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGAtc                                                                                                                                             <  1:161727/149‑1 (MQ=255)      cccATTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGAtc                                                                                                                                             <  8:3468/149‑1 (MQ=255)      cccATTCCGGGGCGCCAGTGAACACCACTGTCAGCCAGCTGTCCGGCCCTTCTCCGGCTACCGCGTCGCCAATGTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGAtc                                                                                                                                             <  1:327426/149‑1 (MQ=255)       ccATTCCGGGTCGCCAGTGAACGACACTTTCAGCCAGCGGTCCGGCCGTTCTCCGCATACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCAAT‑‑‑TAATCGATGGTCTTCGCCGCCATGGTTTTGGGCCGGAGGAAATACCGCGCGaacg                                                                                                                                            >  7:203436/1‑146 (MQ=255)         aTTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGGGCCAGATATCCCATTAGTGAGCGATGGTCTGCGCCGCCATGGTTTTGGGCACGATAAAATACACCACGAt                                                                                                                                              >  1:105367/1‑148 (MQ=255)         aTTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTcc                                                                                                                                          <  7:177962/149‑1 (MQ=255)             cGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCcgc                                                                                                                                       >  2:210724/1‑148 (MQ=255)              gggTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATGTCCcgcg                                                                                                                                      >  7:39288/1‑148 (MQ=255)              gggTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCcgcg                                                                                                                                      >  2:145760/1‑148 (MQ=255)              gggTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCcgcg                                                                                                                                      >  7:44273/1‑148 (MQ=255)               ggTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGaa                                                                                                                                    <  8:203436/149‑1 (MQ=255)                           actcctcTGTCATCCAGCGGTCGGGCCCTTCTCCGGCTACCGCGTCGCCAATCTCGTTGCGCCAGGCTTCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAATTACAGCTGGATCTCGCGCGAACGCCCGACCTTg                                                                                                                        <  5:290002/143‑1 (MQ=255)                            caccacTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTgg                                                                                                                       >  3:274435/1‑149 (MQ=255)                            caccacTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTgg                                                                                                                       <  3:117857/149‑1 (MQ=255)                             accacTCGCAGCCAGCGGGCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCGCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGc                                                                                                                      <  8:238192/149‑1 (MQ=255)                              ccacTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGc                                                                                                                      >  7:21553/1‑148 (MQ=255)                               cacTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGa                                                                                                                    >  1:214979/1‑149 (MQ=255)                               cacTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGa                                                                                                                    <  8:144441/149‑1 (MQ=255)                                    tCAGCCTGCGGTACGGCACTTCTCCGACTACCGCGTCGCAAATCTCGTTGATCCAGGCATCCCAT‑‑‑TCATAGATGTTCTTCGCCGCCTTGGTTTTGGGCACGATGAAATACAGCTCGTTCTCCCGCGAAATCCAGACCTTGTCGACATAg                                                                                                               <  2:214137/149‑1 (MQ=255)                                    tCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAg                                                                                                               >  8:295844/1‑149 (MQ=255)                                    tCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGATCGATCTCCCGCGAACGCCCGACCTTGGCGCCATAg                                                                                                               >  7:23882/1‑149 (MQ=255)                                       gccagcGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGgcg                                                                                                            >  6:182301/1‑149 (MQ=255)                                       gccagcGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGgcg                                                                                                            >  6:283163/1‑149 (MQ=255)                                       gccagcGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGAATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGgcg                                                                                                            >  6:329517/1‑149 (MQ=255)                                         cagcGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGCAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGcgcg                                                                                                          >  7:40679/1‑149 (MQ=255)                                          agcGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCgg                                                                                                         <  1:191736/149‑1 (MQ=255)                                             ggTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGt                                                                                                        >  2:41318/1‑147 (MQ=255)                                               tCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGt                                                                                                        >  7:276508/1‑145 (MQ=255)                                                  ggCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACGGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTg                                                                                                 >  3:103703/1‑149 (MQ=255)                                                    cccTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCa                                                                                               <  8:276508/149‑1 (MQ=255)                                                    cccTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCa                                                                                               <  4:274435/149‑1 (MQ=255)                                                    cccTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCa                                                                                               <  5:182301/149‑1 (MQ=255)                                                      cTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAgg                                                                                             >  8:191022/1‑149 (MQ=255)                                                      cTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAgg                                                                                             <  7:295844/149‑1 (MQ=255)                                                      cTTCTACGCCTACCTCGTCGCCACTCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGGTTTGGGCACGATGAAATACAGCTCGATCTCCCTCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAgg                                                                                             <  2:105367/149‑1 (MQ=255)                                                         ctcCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGCCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGc                                                                                            <  8:269848/147‑1 (MQ=255)                                                            cGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAg                                                                                              >  8:279156/1‑142 (MQ=255)                                                            cGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAg                                                                                              <  7:279156/142‑1 (MQ=255)                                                                   cgcgTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGcga                                                                                <  1:145760/149‑1 (MQ=255)                                                                         gCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAg                                                                           <  1:70353/148‑1 (MQ=255)                                                                         gCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCAGGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAgg                                                                          >  6:270787/1‑149 (MQ=255)                                                                           cAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATACGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGTcc                                                                        <  8:21553/149‑1 (MQ=255)                                                                            aaTCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCt                                                                       >  7:161214/1‑149 (MQ=255)                                                                             aTCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCtt                                                                      >  7:229704/1‑149 (MQ=255)                                                                                tcGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGc                                                                   >  4:121743/1‑149 (MQ=255)                                                                                  gTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACGGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTgg                                                                    <  4:106982/146‑1 (MQ=255)                                                                                  gTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTgg                                                                    <  2:214979/146‑1 (MQ=255)                                                                                  ctggcgcCTTGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGTCGGGCGACGAAGGCCTTg                                                                     <  8:244597/142‑1 (MQ=255)                                                                                   ttGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCAc                                                                <  8:161214/149‑1 (MQ=255)                                                                                     gcgcCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACGGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCAcc                                                               <  4:103703/148‑1 (MQ=255)                                                                                      cgcCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTc                                                             <  1:210724/149‑1 (MQ=255)                                                                                       gcCAGGCATCCCAT‑‑‑TCATCGAAGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACCGCTCGATCTCCCGCGACCGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCa                                                            >  2:233883/1‑149 (MQ=255)                                                                                        ccAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCAt                                                           <  1:41318/149‑1 (MQ=255)                                                                                           ggCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCAGCTa                                                        >  8:73016/1‑149 (MQ=255)                                                                                                cccAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGc                                                   <  3:121743/149‑1 (MQ=255)                                                                                                          aTCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTg                                             <  5:329517/148‑1 (MQ=255)                                                                                                           tCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCatggggcg                                           >  6:301700/1‑142 (MQ=255)                                                                                                           tCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAg                                           <  8:40679/149‑1 (MQ=255)                                                                                                             gATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGt                                         >  3:311796/1‑149 (MQ=255)                                                                                                                     tCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAggg                                 >  8:128286/1‑149 (MQ=255)                                                                                                                       gccgccATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGt                               <  8:229704/149‑1 (MQ=255)                                                                                                                            cATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGt                                                     >  5:269305/1‑122 (MQ=255)                                                                                                                            cATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGt                                                     <  6:269305/122‑1 (MQ=255)                                                                                                                            cATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCa                          <  2:55495/149‑1 (MQ=255)                                                                                                                              tGGGTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAgg                        <  2:222833/149‑1 (MQ=255)                                                                                                                                gTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGaa                      <  7:191022/149‑1 (MQ=255)                                                                                                                                      ggCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGAAGAttt                 >  2:279466/1‑148 (MQ=255)                                                                                                                                           gATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGAAGATTTCAGa             >  1:44429/1‑147 (MQ=255)                                                                                                                                                 aTACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGAAGATTTCAGACAGCGCCt     >  8:294073/1‑149 (MQ=255)                                                                                                                                                  tACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGAAGATTTCAGACAGCGCCt     >  4:277216/1‑148 (MQ=255)                                                                                                                                                   aCAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGAAGATTTCAGACAGCGCCTGa   >  3:315039/1‑149 (MQ=255)                                                                                                                                                   aCAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGAAGATTTCAGACAGCGCCTGa   <  8:39288/149‑1 (MQ=255)                                                                                                                                                   aCAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGAAGATTTCAGACAGCGCCTGa   <  8:44273/149‑1 (MQ=255)                                                                                                                                                    cAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGAAGATTTCAGACAGCGCCTGAc  <  5:270787/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CTCGGCCCATTCCGGGTCGCCAGTGAACACCACTGTCAGCCAGCGGTCCGGCCCTTCTCCGCCTACCGCGTCGCCAATCTCGTTGCGCCAGGCATCCCAT‑‑‑TCATCGATGGTCTTCGCCGCCATGGTTTTGGGCACGATGAAATACAGCTCGATCTCCCGCGAACGCCCGACCTTGGCGACATAGGCGCGGTAAGATTGCAGGCCATGGCGGGCGACGAAGGCCTTGGCCACCTCATCTACGTGCAGCTTGAGGTCGCCAGGGGTGACCAGGAAGATTTCAGACAGCGCCTGAC  >  NC_002947/1402578‑1402870

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