Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,648,558	C→T	5.8%	P209P (CCG→CCA)	PP_4112 ←	sulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,648,558	0	C	T	5.8%	202.0 / 2.5	69	P209P (CCG→CCA)	PP_4112	sulfate transporter
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base T (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

AGCCACGGCTCTCGGTGATTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATAAACCACCGCCATGG  >  NC_002947/4648416‑4648706                                                                                                                                               |                                                                                                                                                     aGCCACGGCTCTCGGTGATTTCGTCGGTCAGGTTGAGGGTGAGCACCGATGTTATCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGGTTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGc                                                                                                                                               <  1:1573895/150‑1 (MQ=255)    cACGGCTCTCGGTGATTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCgg                                                                                                                                                    <  2:1601877/142‑1 (MQ=255)    cACGGCTCTCGGTGATTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCgg                                                                                                                                                    >  1:1601877/1‑142 (MQ=255)             cGGTGATTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGAtgt                                                                                                                                   <  2:1302014/150‑1 (MQ=255)              ggTGATTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGAtgtg                                                                                                                                  >  1:652588/1‑150 (MQ=255)                tGATTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTgcg                                                                                                                                >  1:1861740/1‑150 (MQ=255)                              aGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAgg                                                                                                                  <  1:1535216/150‑1 (MQ=255)                                  tGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGc                                                                                                              <  2:1320046/150‑1 (MQ=255)                                    aGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCgg                                                                                                            <  2:56851/150‑1 (MQ=255)                                    aGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCgg                                                                                                            <  2:403562/150‑1 (MQ=255)                                      ggTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGt                                                                                                          >  1:1765173/1‑150 (MQ=255)                                       gTGAGCAGGGTTTCCCGAAGGCCGGCCATGGTCATCAGCACTGTGTGGGGGGTGTTTATCGTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTa                                                                                                         <  2:424508/150‑1 (MQ=255)                                       gTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCCCATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCCGCGTGCGGTTa                                                                                                         <  2:1528929/150‑1 (MQ=255)                                          aGCAGTGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGCCGCCCAGAGTGCGGGTAGGc                                                                                                      <  2:573190/150‑1 (MQ=255)                                           gCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCa                                                                                                     >  1:515995/1‑150 (MQ=255)                                                        cagGCCGACCATGGCCATCAGCACTGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTa                                                                                        <  2:23601/150‑1 (MQ=255)                                                        cagGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTa                                                                                        <  2:515995/150‑1 (MQ=255)                                                          ggCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCCCCAGCCCCGGCCGGCCCCCGGAGCTGTACGCCAAGTCGCCTATGTTGCTGGTTGGCAGGCCGAGCACGTaca                                                                                      <  2:805466/150‑1 (MQ=255)                                                           gCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTacac                                                                                     >  1:55367/1‑150 (MQ=255)                                                              gACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTAcaccac                                                                                  >  1:139966/1‑150 (MQ=255)                                                               aCCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCACGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACaccacc                                                                                 >  1:1847794/1‑150 (MQ=255)                                                                 cATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAg                                                                               <  1:1829286/150‑1 (MQ=255)                                                                    ggCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGcc                                                                            <  2:55367/150‑1 (MQ=255)                                                                      ccATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCaa                                                                          >  1:1686743/1‑150 (MQ=255)                                                                       cATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAac                                                                         >  2:354252/1‑150 (MQ=255)                                                                       cATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAac                                                                         <  2:139966/150‑1 (MQ=255)                                                                         tCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAacac                                                                       >  2:783508/1‑150 (MQ=255)                                                                         tCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGGACACCACCAGGCCAacac                                                                       >  1:1554105/1‑150 (MQ=255)                                                                               ccGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAggg                                                                 >  2:1235354/1‑149 (MQ=255)                                                                                   gTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGc                                                             >  1:983250/1‑150 (MQ=255)                                                                                      gggggCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGccac                                                          <  2:1012164/150‑1 (MQ=255)                                                                                        gggCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTc                                                                                                                         <  2:824552/85‑1 (MQ=255)                                                                                        gggCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTc                                                                                                                         >  1:824552/1‑85 (MQ=255)                                                                                            gatgatCTTCAGGGTTTCCAGATTCCAAGGTACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGc                                                    <  2:1672446/150‑1 (MQ=255)                                                                                             atgatCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCg                                                                                                             >  1:438418/1‑92 (MQ=255)                                                                                             atgatCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCg                                                                                                             <  2:438418/92‑1 (MQ=255)                                                                                             atgatCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGcc                                                   >  2:912187/1‑150 (MQ=255)                                                                                                atCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCcggc                                                <  2:983250/150‑1 (MQ=255)                                                                                                    tCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCACCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCa                                            >  2:1568476/1‑150 (MQ=255)                                                                                                       gggTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACcg                                         >  1:1704239/1‑150 (MQ=255)                                                                                                         gTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACcgcg                                       >  2:1490015/1‑150 (MQ=255)                                                                                                              cAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGccc                                                                                                 >  1:211880/1‑87 (MQ=255)                                                                                                              cAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGccc                                                                                                 <  2:211880/87‑1 (MQ=255)                                                                                                              cAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGt                                  <  2:1704239/150‑1 (MQ=255)                                                                                                               aGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTc                                 >  1:1348853/1‑150 (MQ=255)                                                                                                                gATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGCTGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAgg                                                                  >  1:1505292/1‑116 (MQ=255)                                                                                                                gATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGa                                                                                               >  1:644029/1‑87 (MQ=255)                                                                                                                gATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGa                                                                                               <  2:644029/87‑1 (MQ=255)                                                                                                                gATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAgg                                                                  <  2:1505292/116‑1 (MQ=255)                                                                                                                 aTTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAg                               <  2:1861740/150‑1 (MQ=255)                                                                                                                  ttCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGt                              <  2:1220429/150‑1 (MQ=255)                                                                                                                   tCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGtt                             <  1:1490015/150‑1 (MQ=255)                                                                                                                     cAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTc                           <  2:1348853/150‑1 (MQ=255)                                                                                                                     cAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTACGACATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTc                           >  2:23354/1‑150 (MQ=255)                                                                                                                      aaGGCACATCCGGCAGCGCCAGCCCTGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGAt                                                                >  1:1291579/1‑112 (MQ=255)                                                                                                                      aaGGCACATCCGGCAGCGCCAGCCCTGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGAt                                                                <  2:1291579/112‑1 (MQ=255)                                                                                                                      aaGGCACATCCGGCAGCGCCAGCCCTGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGAt                                                                >  1:1291546/1‑112 (MQ=255)                                                                                                                      aaGGCACATCCGGCAGCGCCAGCCCTGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGAt                                                                <  2:1291546/112‑1 (MQ=255)                                                                                                                      aaGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCg                          >  2:886362/1‑150 (MQ=255)                                                                                                                       aGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCgg                         <  2:1765173/150‑1 (MQ=255)                                                                                                                       aGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCgg                         >  2:606440/1‑150 (MQ=255)                                                                                                                       aGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCACAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCgg                         >  1:1809524/1‑150 (MQ=255)                                                                                                                                  gCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGTGGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATaaa              >  1:637672/1‑150 (MQ=255)                                                                                                                                  gCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATaaa              <  1:912187/150‑1 (MQ=255)                                                                                                                                   cAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCCGGCCAACCCCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATAAac             <  2:1469128/150‑1 (MQ=255)                                                                                                                                     gcgcCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGAACTCCACCAATGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATAAacca           >  1:1886666/1‑150 (MQ=255)                                                                                                                                       gcCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATAAaccacc         >  1:1102314/1‑150 (MQ=255)                                                                                                                                        ccAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATAAACCACCg        >  1:1775134/1‑150 (MQ=255)                                                                                                                                         cAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATAAACCACCGc       >  1:1895674/1‑150 (MQ=255)                                                                                                                                            ccccGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATAAACCACCGCCAt    >  1:699293/1‑150 (MQ=255)                                                                                                                                             cccGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATAAACCACCGCCATg   <  1:886362/150‑1 (MQ=255)                                                                                                                                              ccGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCGCCGCGCGGGTCAGTTTCGGCAGCACATAAACCACCGCCATgg  >  1:1059836/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                     AGCCACGGCTCTCGGTGATTTCGTCGGTCAGGTTGAGGGTGAGCAGGGTTTCCAGCAGGCCGACCATGGCCATCAGCACCGCGTAGGGGGCGATGATCTTCAGGGTTTCCAGATTCCAAGGCACATCCGGCAGCGCCAGCCCCGGCAGGCCACCGGCGATGTGCGCCATGTCGCCCAGGGTGCGGGTAGGCAGCCCGAGCACGTACACCACCAGGCCAACACCGAGGATCGCCACCAAGGCCGGCGGCACCGCGCGGGTCAGTTTCGGCAGCACATAAACCACCGCCATGG  >  NC_002947/4648416‑4648706

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