Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	618,358	A→G	5.9%	S71P (TCG→CCG)	PP_0532 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	618,358	0	A	G	5.9%	250.7 / 4.1	84	S71P (TCG→CCG)	PP_0532	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (45/34);  new base G (3/2);  total (48/36)
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.29e-01

GGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTCATGCTCGGGCTGG  >  NC_002947/618211‑618506                                                                                                                                                    |                                                                                                                                                     ggccggccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGAt                                                                                                                                                     <  1:250756/149‑1 (MQ=255)      gccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTgg                                                                                                                                                <  4:350302/149‑1 (MQ=255)      gccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTgg                                                                                                                                                <  3:354593/149‑1 (MQ=255)      gccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTgg                                                                                                                                                >  6:113849/1‑149 (MQ=255)       cccTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCATCGGGGTACCCGGTGCTgtg                                                                                                                                               >  8:146572/1‑147 (MQ=255)       cccTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGgc                                                                                                                                               <  2:230347/149‑1 (MQ=255)        ccTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCATCGGGGTACCCGGTGCTgtgg                                                                                                                                              <  7:146572/149‑4 (MQ=255)            tCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTg                                                                                                                                           >  6:148726/1‑148 (MQ=255)            tCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTg                                                                                                                                           >  6:154984/1‑148 (MQ=255)            tCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTg                                                                                                                                           <  6:266050/148‑1 (MQ=255)              gATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGa                                                                                                                                         <  4:333542/148‑1 (MQ=255)              gATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGAc                                                                                                                                        >  4:249555/1‑149 (MQ=255)                tGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGa                                                                                                                                      >  5:128257/1‑149 (MQ=255)                           gTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTgg                                                                                                                                                >  2:370013/1‑128 (MQ=255)                           gTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTgg                                                                                                                                                <  1:370013/128‑1 (MQ=255)                           gTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGcc                                                                                                                           <  5:389202/149‑1 (MQ=255)                           gTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGc                                                                                                                            >  4:165112/1‑148 (MQ=255)                           gTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCACCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGc                                                                                                                            >  8:10658/1‑148 (MQ=255)                               ggTAACCCACCGCCACCAGCAGGGGCACCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAAcc                                                                                                                       >  5:360569/1‑149 (MQ=255)                                 tAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCca                                                                                                                     <  3:165112/149‑1 (MQ=255)                                     ccACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCCACCCACAgg                                                                                                                 >  2:173891/1‑149 (MQ=255)                                     ccACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAgg                                                                                                                 >  8:29881/1‑149 (MQ=255)                                     ccACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATAGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAgg                                                                                                                 >  4:37037/1‑149 (MQ=255)                                      cACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCAATCTAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGc                                                                                                                >  6:258696/1‑149 (MQ=255)                                       aCCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCg                                                                                                               >  3:82064/1‑149 (MQ=255)                                        ccGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCgg                                                                                                              <  3:249555/149‑1 (MQ=255)                                              ccaGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCCACCCACAGGCGGTAGTGa                                                                                                        >  6:260615/1‑149 (MQ=255)                                               cagcagGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGAc                                                                                                       <  7:123052/149‑1 (MQ=255)                                               cagcagGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGAc                                                                                                       <  5:346365/149‑1 (MQ=255)                                                    ggggCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGACTTGCCAACCCACAGGCGGTAGTGACGCAgg                                                                                                  <  6:333283/149‑1 (MQ=255)                                                     gggCACCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGc                                                                                                 <  6:360569/149‑1 (MQ=255)                                                              gCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGcc                                                                                         <  5:113849/148‑1 (MQ=255)                                                              gCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGc                                                                                          <  6:128257/147‑1 (MQ=255)                                                              gCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGc                                                                                          <  6:37725/147‑1 (MQ=255)                                                                 cACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGt                                                                                     >  6:71052/1‑149 (MQ=255)                                                                      cGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGa                                                                                <  7:10658/149‑1 (MQ=255)                                                                        aGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATc                                                                              <  4:82064/149‑1 (MQ=255)                                                                             gTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCg                                                                                            <  7:203389/130‑1 (MQ=255)                                                                             gTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCg                                                                                            >  8:203389/1‑130 (MQ=255)                                                                              tCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTg                                                                        >  4:168212/1‑149 (MQ=255)                                                                                cAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTggg                                                                      >  3:88595/1‑149 (MQ=255)                                                                                cAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTgg                                                                       <  5:154984/148‑1 (MQ=255)                                                                                cAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTgg                                                                       <  5:148726/148‑1 (MQ=255)                                                                                   ttGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTg                                                                    >  4:218017/1‑148 (MQ=255)                                                                                     gctgctGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAg                                                                 >  4:12018/1‑149 (MQ=255)                                                                                      ctgctgATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAg                                                                 <  1:205890/148‑1 (MQ=255)                                                                                        gctgATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGa                                                              >  4:168227/1‑149 (MQ=255)                                                                                          tgATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGa                                                              >  7:176042/1‑147 (MQ=255)                                                                                            aTCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAgg                                                          >  3:204854/1‑149 (MQ=255)                                                                                            aTCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAg                                                           <  7:127200/148‑1 (MQ=255)                                                                                                aGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAg                                                      >  2:122216/1‑149 (MQ=255)                                                                                                    tCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTcc                                                  >  1:129758/1‑149 (MQ=255)                                                                                                     ccccAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCAGATTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCt                                                 >  1:123499/1‑149 (MQ=255)                                                                                                       ccAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCCGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTgg                                               >  5:268857/1‑149 (MQ=255)                                                                                                       ccACCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTAGCGCTGGACCATACGGCTTGCCAACCCCCAGGCGGTAGTGACGCAGCCCATCGGCCGAGTCCACATCAAGCTGGGTGAAGAGATAGGCCAGGACCTgg                                               >  4:322383/1‑149 (MQ=255)                                                                                                           cgcAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCt                                           <  3:218017/149‑1 (MQ=255)                                                                                                             cagcaTTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTg                                          >  5:58702/1‑148 (MQ=255)                                                                                                                 aTTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAaca                                     >  2:64056/1‑149 (MQ=255)                                                                                                                   tCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACAcc                                   >  5:37480/1‑149 (MQ=255)                                                                                                                     ccATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGa                                 >  3:284758/1‑149 (MQ=255)                                                                                                                      cATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGAc                                >  5:134401/1‑149 (MQ=255)                                                                                                                      cATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGAc                                <  8:176042/149‑1 (MQ=255)                                                                                                                      cATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGAc                                >  6:148692/1‑149 (MQ=255)                                                                                                                      cATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGAc                                >  3:256659/1‑149 (MQ=255)                                                                                                                       aTCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCa                                                                                                    >  3:254595/1‑80 (MQ=255)                                                                                                                       aTCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCa                                                                                                    <  4:254595/80‑1 (MQ=255)                                                                                                                        tCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGATTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGAGCTGTAACACCGACgt                              >  6:160545/1‑149 (MQ=255)                                                                                                                          gAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACgtgt                            <  3:168212/149‑1 (MQ=255)                                                                                                                             cATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCa                         >  4:330114/1‑149 (MQ=255)                                                                                                                                  ccaGCATCGGGGTACCCGGTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtcttc                    >  4:164715/2‑149 (MQ=255)                                                                                                                                  acaGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGAt                                                                               >  8:367993/1‑90 (MQ=255)                                                                                                                                  acaGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGAt                                                                               <  7:367993/90‑1 (MQ=255)                                                                                                                                  acaGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACCCCGACGTGTCCAtcttc                    >  7:152519/1‑149 (MQ=255)                                                                                                                                   caGCATCGGGGTACCCGGTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGCTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtcttct                   >  3:40502/1‑149 (MQ=255)                                                                                                                                    aGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTg                  <  3:37037/149‑1 (MQ=255)                                                                                                                                    aGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTg                  <  2:129758/149‑1 (MQ=255)                                                                                                                                       aCCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTc               <  6:58702/149‑1 (MQ=255)                                                                                                                                       aCCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTc               <  6:134401/149‑1 (MQ=255)                                                                                                                                                cccGGTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTCATGCTCggg      <  3:164715/149‑1 (MQ=255)                                                                                                                                                cccGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTCATGCTCggg      >  6:172332/1‑149 (MQ=255)                                                                                                                                                cccGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTCATGCTCggg      >  6:140238/1‑149 (MQ=255)                                                                                                                                                 ccGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGAc                                <  7:49728/122‑1 (MQ=255)                                                                                                                                                 ccGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGAc                                >  8:49728/1‑122 (MQ=255)                                                                                                                                                 ccGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACTTGTCCATCTGATGGGCAGGCTCgg       >  1:15731/1‑147 (MQ=255)                                                                                                                                                    aTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTCATGCTCGGGCTgg  <  5:258696/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTCATGCTCGGGCTGG  >  NC_002947/618211‑618506

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