Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	346,118	G→A	9.2%	N576N (AAC→AAT) ‡	PP_0287 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	346,118	0	G	A	9.2%	279.6 / 13.8	98	N576N (AAC→AAT) ‡	PP_0287	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/48);  new base A (7/2);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 8.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-01

CGACGCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCG  >  NC_002947/345970‑346266                                                                                                                                                     |                                                                                                                                                     cGACGCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCg                                                                                                                                                      <  3:164650/149‑1 (MQ=255)     gCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGa                                                                                                                                                  >  8:48137/1‑149 (MQ=255)     gCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGa                                                                                                                                                  <  5:296961/149‑1 (MQ=255)       cATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATg                                                                                                                                                >  3:209529/1‑149 (MQ=255)        aatGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATgg                                                                                                                                               >  8:181922/3‑149 (MQ=255)         ttGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGc                                                                                                                                              <  1:263011/149‑1 (MQ=255)            cGTACCCCCCGGCTGCCGCGCCGCCCCTGGCAGGGGGGGTCCTTGCCGCCTGGTTCGCTGCTCACCGTCATGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGaa                                                                                                                                           <  7:232269/149‑1 (MQ=255)            cGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGTGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGaa                                                                                                                                           <  1:253503/149‑1 (MQ=255)            cGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGaa                                                                                                                                           <  7:181922/149‑1 (MQ=255)            cGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGaa                                                                                                                                           <  3:208039/149‑1 (MQ=255)             gTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAg                                                                                                                                          >  8:240790/1‑149 (MQ=255)              tACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGt                                                                                                                                         >  7:4752/1‑149 (MQ=255)                   ccaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgcg                                                                                                                                    >  4:333585/1‑149 (MQ=255)                      ggCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGGGAAGTTTGCGCTg                                                                                                                                 >  7:267491/1‑149 (MQ=255)                      ggCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTg                                                                                                                                 <  3:343020/149‑1 (MQ=255)                         tGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCg                                                                                                                              >  6:258137/1‑149 (MQ=255)                         tGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCg                                                                                                                              >  6:242886/1‑149 (MQ=255)                         tGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCg                                                                                                                              >  5:309518/1‑149 (MQ=255)                           ccGCGCAGCTCCTGTAAGGGGGTGTCCTTGCGGCGTGGTTCGCTTCTCAGCGTCTTGCGATTGAGCGTGTCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCAGGCCATCGTTGATGGCGAAGTTTGCGCTGCCGtt                                                                                                                            <  1:262212/149‑1 (MQ=255)                            cgcgCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCGGCGTCTTGCGATTGAGCGTGGCGATGGATTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCAACGTTGATGGCGAAGTTTGCGGTGCCGtt                                                                                                                            >  5:108821/1‑148 (MQ=255)                              cgcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                         <  1:373316/149‑1 (MQ=255)                              cgcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                         <  6:309518/149‑1 (MQ=255)                              cgcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                         <  3:250930/149‑1 (MQ=255)                               gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                        <  7:240790/149‑1 (MQ=255)                               gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                        <  8:4752/149‑1 (MQ=255)                                    tCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgcgc                                                                                                                                   >  3:254697/1‑133 (MQ=255)                                    tCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgcgc                                                                                                                                   <  4:254697/133‑1 (MQ=255)                                      cTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGa                                                                                                                 >  6:320438/1‑149 (MQ=255)                                       tGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                >  7:344873/1‑149 (MQ=255)                                        ggAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                >  6:24788/1‑148 (MQ=255)                                         gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                <  6:264138/147‑1 (MQ=255)                                         gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                <  8:344873/147‑1 (MQ=255)                                           aGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGt                                                                                                            <  5:258137/149‑1 (MQ=255)                                             ggggTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTg                                                                                                           <  5:252897/148‑1 (MQ=255)                                                   tCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTg                                                                                                    >  7:126545/1‑149 (MQ=255)                                                     cTTGCCGTGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGc                                                                                                  >  4:85906/1‑149 (MQ=255)                                                      ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCt                                                                                                                                  <  3:60730/116‑1 (MQ=255)                                                      ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCt                                                                                                                                  >  4:60730/1‑116 (MQ=255)                                                      ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGc                                                                                                                                >  1:243836/1‑118 (MQ=255)                                                      ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGc                                                                                                                                <  2:243836/118‑1 (MQ=255)                                                      ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCt                                                                                                 >  5:180602/1‑149 (MQ=255)                                                        gCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCt                                                                                                 <  6:103368/147‑1 (MQ=255)                                                           gcgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                        >  2:315764/1‑121 (MQ=255)                                                            cgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                        >  3:45065/1‑120 (MQ=255)                                                            cgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                        <  4:45065/120‑1 (MQ=255)                                                                gTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGTCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTg                                                                                       <  5:42953/149‑1 (MQ=255)                                                                 ttCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTgg                                                                                      >  6:62976/1‑149 (MQ=255)                                                                    gctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGTCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCt                                                                                    <  4:135318/148‑1 (MQ=255)                                                                     ctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                  <  5:242886/149‑1 (MQ=255)                                                                     ctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                  >  7:165169/1‑149 (MQ=255)                                                                     ctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                  <  8:126545/149‑1 (MQ=255)                                                                        ctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                  >  1:241117/1‑146 (MQ=255)                                                                         tCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                  >  1:75996/1‑145 (MQ=255)                                                                           aGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGAt                                                                             <  1:317902/148‑1 (MQ=255)                                                                            gCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCa                                                                           <  8:165169/149‑1 (MQ=255)                                                                               tCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTg                                                                        <  4:294227/149‑1 (MQ=255)                                                                                cTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTgg                                                                       >  1:164339/1‑149 (MQ=255)                                                                                   gCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGCCCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTgg                                                                       <  4:209529/146‑1 (MQ=255)                                                                                   gCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtc                                                                     <  7:252817/148‑1 (MQ=255)                                                                                    cGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtca                                                                    <  5:81379/148‑1 (MQ=255)                                                                                     gATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagg                                                                  <  1:13786/149‑1 (MQ=255)                                                                                          aGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagggtcag                                                             <  3:333585/149‑1 (MQ=255)                                                                                            cGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTcagca                                                           <  2:59788/149‑1 (MQ=255)                                                                                             gTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTCACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTcagcag                                                          <  5:320438/149‑1 (MQ=255)                                                                                              tGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCt                                                                                    <  5:332139/122‑1 (MQ=255)                                                                                              tGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCt                                                                                    >  6:332139/1‑122 (MQ=255)                                                                                                            cacaACTGCTGTTCCAGGTTGGCATAGACCAGCACGCCATCGTTGATGGCGAAGATTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACcggcgg                                           >  4:255415/1‑149 (MQ=255)                                                                                                               aaCTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGgcggcg                                         >  8:116426/1‑148 (MQ=255)                                                                                                                  tgctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTg                                     >  6:2095/1‑149 (MQ=255)                                                                                                                    ctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTg                                     <  2:75996/147‑1 (MQ=255)                                                                                                                       ttCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCt                                >  2:4262/1‑149 (MQ=255)                                                                                                                         ccAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCtt                               <  7:48137/148‑1 (MQ=255)                                                                                                                             gTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCtt                          >  3:179128/1‑149 (MQ=255)                                                                                                                                ggCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCg                        <  2:241117/148‑1 (MQ=255)                                                                                                                                ggCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCg                        <  6:180602/148‑1 (MQ=255)                                                                                                                                ggCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCg                        <  5:62976/148‑1 (MQ=255)                                                                                                                                   aTTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTc                     >  4:276984/1‑148 (MQ=255)                                                                                                                                   aTTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTc                     <  3:276984/148‑1 (MQ=255)                                                                                                                                   aTTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTACAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCt                    >  2:284836/1‑149 (MQ=255)                                                                                                                                     ctacgaGCACGCCATCGTTGCTGGCGAAGTTTGCGCTGCCGTTCGGGGTGGCGACCAGTGCTTTGTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCtt                   <  8:63066/143‑1 (MQ=255)                                                                                                                                       aCCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat                >  6:108406/1‑149 (MQ=255)                                                                                                                                        ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCGGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCCGCGGCGTCTGTTCCTTGCCTTCGGTCtt                   >  7:63066/1‑145 (MQ=255)                                                                                                                                        ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGa                                                  >  5:87216/1‑114 (MQ=255)                                                                                                                                        ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGa                                                  <  6:87216/114‑1 (MQ=255)                                                                                                                                        ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat                >  1:226701/1‑148 (MQ=255)                                                                                                                                        aaacttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGGCGtgat                >  4:263116/7‑148 (MQ=255)                                                                                                                                         aacttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgatg               >  6:128337/6‑148 (MQ=255)                                                                                                                                          aGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGaa             <  7:116426/149‑1 (MQ=255)                                                                                                                                          aGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGaa             <  5:108406/149‑1 (MQ=255)                                                                                                                                           cttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTCCCTTGCCTTCGGTCTTGATGAAg            >  3:273920/4‑149 (MQ=255)                                                                                                                                            ttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagggt                                                                <  6:142767/94‑1 (MQ=255)                                                                                                                                            ttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagggt                                                                >  5:142767/3‑96 (MQ=255)                                                                                                                                            ttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCTCTTTACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgt           >  6:6598/3‑149 (MQ=255)                                                                                                                                             tcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTg                                                                        >  3:283479/2‑87 (MQ=255)                                                                                                                                             tcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTg                                                                        <  4:283479/86‑1 (MQ=255)                                                                                                                                             tcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGACCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTATTGATGAAgtg          >  6:110987/2‑149 (MQ=255)                                                                                                                                              cGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgtgt         <  5:24788/149‑1 (MQ=255)                                                                                                                                               gCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGtt        >  2:106373/1‑149 (MQ=255)                                                                                                                                                ccATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGtt        >  7:97282/1‑148 (MQ=255)                                                                                                                                                     gTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCg  >  8:306521/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CGACGCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCG  >  NC_002947/345970‑346266

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