Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,890,438	T→C	6.6%	R58G (AGG→GGG)	PP_2544 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,890,438	0	T	C	6.6%	210.3 / 2.6	75	R58G (AGG→GGG)	PP_2544	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/38);  new base C (3/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGGGCGTGCG  >  NC_002947/2890291‑2890585                                                                                                                                                    |                                                                                                                                                    cATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCtt                                                                                                                                                    >  5:328455/1‑149 (MQ=255) cATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCtt                                                                                                                                                    >  3:131190/1‑149 (MQ=255)  aTGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTg                                                                                                                                                   <  7:19145/149‑1 (MQ=255)  aTGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTg                                                                                                                                                   >  6:52104/1‑149 (MQ=255)   tGGCCTTCATGCCTGACTACCAGATCCTGTTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGTCGGCTCCCTTgg                                                                                                                                                  <  7:123588/149‑1 (MQ=255)     gCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTggg                                                                                                                                                 >  8:287799/1‑148 (MQ=255)      ccTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGa                                                                                                                                               >  4:10458/1‑149 (MQ=255)       cTTCATGCCTGACTACAAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATCGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCAGGCTCCCTTGGGGAc                                                                                                                                              >  3:330788/1‑149 (MQ=255)        ttCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGAcc                                                                                                                                             >  7:226800/1‑149 (MQ=255)            tGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGc                                                                                                                                         <  7:120196/149‑1 (MQ=255)              ccTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGcc                                                                                                                                        <  3:27949/148‑1 (MQ=255)               cTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGt                                                                                                                                      >  3:53827/1‑149 (MQ=255)               cTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGt                                                                                                                                      >  6:358234/1‑149 (MQ=255)                      ccAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTgg                                                                                                                                                  <  8:184951/130‑1 (MQ=255)                      ccAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTgg                                                                                                                                                  >  7:184951/1‑130 (MQ=255)                      ccAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCgt                                                                                                                               >  4:335621/1‑149 (MQ=255)                             cTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTgg                                                                                                                         <  3:10458/148‑1 (MQ=255)                              tGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGc                                                                                                                       >  1:139283/1‑149 (MQ=255)                              tGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGc                                                                                                                       >  5:140269/1‑149 (MQ=255)                                 tCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTg                                                                                                                    <  5:297578/149‑1 (MQ=255)                                 tCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTg                                                                                                                    >  4:145853/1‑149 (MQ=255)                                      tCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAATGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCt                                                                                                                <  1:294960/148‑1 (MQ=255)                                      tCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCt                                                                                                                <  1:289855/148‑1 (MQ=255)                                         aCCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCttttt                                                                                                            <  7:69518/149‑1 (MQ=255)                                           cGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCATTTCGCGTTGTTGGGCCTGGGCTTTTTCt                                                                                                          >  5:347831/1‑149 (MQ=255)                                             gTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGc                                                                                                        <  6:347831/149‑1 (MQ=255)                                                tCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTGCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCggg                                                                                                     >  5:210947/1‑149 (MQ=255)                                                tCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCgg                                                                                                      <  5:115654/148‑1 (MQ=255)                                                tCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCCGCgg                                                                                                      <  7:287799/148‑1 (MQ=255)                                                 cGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCggg                                                                                                     >  3:381399/1‑148 (MQ=255)                                                  gCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGc                                                                                                    <  4:381399/148‑1 (MQ=255)                                                  gCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCa                                                                                                   >  8:120701/1‑149 (MQ=255)                                                    tGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCACGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCGTTTTCTGCGGGCAtg                                                                                                 >  7:374039/1‑149 (MQ=255)                                                     gcggcgGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCAtgt                                                                                                <  3:94633/149‑1 (MQ=255)                                                      cggcggGCGTGGTGGCGTTGGGTGTGCCGCACCACCTCCGCTCGGGGAGCCAGCAGCACTAAGCCGCACTCGCCTCTTGCAACCGCCGGCTCCCCTTGGGACCCGGCCGGTTCGCGTTGTTGGCCCTGGGCTTTTTCTGCGGGCAcgtg                                                                                               <  4:330788/149‑1 (MQ=255)                                                        gcggTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGa                                                                                             <  7:313706/149‑1 (MQ=255)                                                          ggTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGCCCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGaaa                                                                                           >  8:159952/1‑149 (MQ=255)                                                           gTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAAc                                                                                          <  5:52104/149‑1 (MQ=255)                                                           gTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCACCCCTGGGGCACCGGCCGTTTAGCGGTGTTGGGCCTGGGCCTTTTCTTCGGGCATGTGGgaaa                                                                                          >  5:59407/1‑148 (MQ=255)                                                             cGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGc                                                                                        >  3:276591/1‑149 (MQ=255)                                                              gtggtgGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGAGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCt                                                                                       >  7:227299/1‑149 (MQ=255)                                                                ggtggCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTgg                                                                                     <  4:276591/149‑1 (MQ=255)                                                                    gCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCa                                                                                 <  4:64043/149‑1 (MQ=255)                                                                     cGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCAt                                                                                <  2:139283/149‑1 (MQ=255)                                                                      gATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATc                                                                               <  3:335621/149‑1 (MQ=255)                                                                       aTGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATc                                                                               >  4:354428/1‑148 (MQ=255)                                                                       aTGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCGGGGCGTTTTCTGCGGGCATGTGGAAACGCTGGttactc                                                                               >  4:354433/1‑144 (MQ=255)                                                                         gggTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGcc                                                                                                                                        >  6:26172/1‑89 (MQ=37)                                                                            tgtgGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGt                                                                         <  6:328455/149‑1 (MQ=255)                                                                                 tgCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAAcc                                                                    <  2:77409/149‑1 (MQ=255)                                                                                 tgCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAAcc                                                                    <  5:358234/149‑1 (MQ=255)                                                                                    aCCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGGCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAg                                                                 <  5:112176/149‑1 (MQ=255)                                                                                     ccAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAg                                                                 >  8:116999/1‑148 (MQ=255)                                                                                      cAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCa                                                               <  3:145853/149‑1 (MQ=255)                                                                                        gATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAgg                                                             <  4:53827/149‑1 (MQ=255)                                                                                             gCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTc                                                                                                           >  6:308113/1‑98 (MQ=255)                                                                                             gCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTc                                                                                                           <  5:308113/98‑1 (MQ=255)                                                                                             gCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAg                                                              >  1:300781/1‑143 (MQ=255)                                                                                             gCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAg                                                              <  2:300781/143‑1 (MQ=255)                                                                                             gCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGc                                                        <  4:296533/149‑1 (MQ=255)                                                                                                   gTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCGTTGGGCCCCGGCCGTTTCGCGGTGTTGTGCGTGGGGTTTTTATTCGGGGATGTGGTAACGCTGGTTCATCCACGGTGATACGAGCAGGTGTGCCGCAGc                                                  >  8:359148/1‑149 (MQ=255)                                                                                                      gcaagcaGCACTAAGCGTCCCTCGCCTCTGGAACCGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCTAt                                               <  3:354433/149‑1 (MQ=255)                                                                                                      gcaagcaGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGAt                                               <  3:354428/149‑1 (MQ=255)                                                                                                          gcagcaCTAAGCGGCACTCGCCTCTGGAACCGGCCGGCGCCCTTGGGGACCCGGCCGTTTCGCGGTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAgg                                           <  8:374039/149‑1 (MQ=255)                                                                                                            agcaCTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGc                                         <  8:227299/149‑1 (MQ=255)                                                                                                               aCTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAg                                      <  7:120701/149‑1 (MQ=255)                                                                                                                 tAAGCGTCACTCGCCTCTGGAAACGGCCGGCTACCGAGGGGACCCCGCCCTTTGGCGTTGTTGGGCCTGGGCTTTTTGAGCGGGGATTTGGAAACGCAGGTTAAAGCAAGCTGAACCGAGCCTCGCGTCCGCAGCGCTCACGGGCCGGa                                    >  8:336265/1‑149 (MQ=255)                                                                                                                   tgCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATg                                  <  6:59407/148‑1 (MQ=255)                                                                                                                       tCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCCGCCGTTTCGCGTTGTTGGGCCGGGGCGTTTTCTGCGGGCATGTGGGAACGCTGGTTGATCCACGGTGGACCGAGCAGGGCTGCCGCAGCGCTCAGGGGCAGGATGGc                                >  7:312877/1‑147 (MQ=255)                                                                                                                        cACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGGTTCGCGTTGTTGGGCCGGGGGTTTTTCTGCGGGCATGTGGAAACGCTGGTTGATCCCCGGTGCACCGAGCAGGGCTGCCGCAGCGCTCCGGGCCAGGATGGCGAg                             >  2:210934/1‑149 (MQ=255)                                                                                                                          cTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAggg                           >  2:93432/1‑149 (MQ=255)                                                                                                                             gCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGa                        <  4:131190/149‑1 (MQ=255)                                                                                                                               ctctGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATg                      <  6:37453/149‑1 (MQ=255)                                                                                                                                  gcgAAACGGCCGGGTCCCCAAGGGAGCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAg                   >  2:19018/3‑149 (MQ=39)                                                                                                                                                  ccTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGGGCGt     <  6:210947/147‑1 (MQ=255)                                                                                                                                                  ccTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGGGCGTgc   >  5:147381/1‑149 (MQ=255)                                                                                                                                                   cTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTACGTCTTCACGGGCGTgcg  >  6:194551/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAGATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGCATGTGGAAACGCTGGTTCATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGGGCGTGCG  >  NC_002947/2890291‑2890585

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