Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,647,851	C→A	6.9%	D430Y (GAC→TAC) ‡	PP_4112 ←	sulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,647,851	0	C	A	6.9%	227.9 / 4.6	73	D430Y (GAC→TAC) ‡	PP_4112	sulfate transporter
Reads supporting (aligned to +/- strand):  ref base C (29/38);  new base A (3/2);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

ACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGAGC  >  NC_002947/4647704‑4647993                                                                                                                                                    |                                                                                                                                               aCAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGCCCGAATAATCGACAAAGCTCAAATCCTGGAAATCCAGCGTGACCTGGGCCGGGTcc                                                                                                                                               >  6:289387/1‑149 (MQ=255)  cAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTccc                                                                                                                                              >  7:1498/1‑149 (MQ=255)    gCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCg                                                                                                                                            <  3:90507/149‑1 (MQ=255)     cttcttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCgg                                                                                                                                           <  8:1498/149‑1 (MQ=255)     cttcttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCgg                                                                                                                                           <  2:310144/149‑1 (MQ=255)        cttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGcc                                                                                                                                         <  3:213062/148‑1 (MQ=255)           gCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGGACCCGGcccag                                                                                                                                      <  7:184239/148‑4 (MQ=255)            cACCGCTCCGACAAATTCACCACGCTCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTc                                                                                                                                    <  6:234714/149‑1 (MQ=255)            cACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTc                                                                                                                                    >  8:104427/1‑149 (MQ=255)             aCCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTc                                                                                                                                    >  4:275621/1‑148 (MQ=255)             aCCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCg                                                                                                                                   <  1:124823/149‑1 (MQ=255)                gCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAAc                                                                                                                                <  6:321431/149‑1 (MQ=255)                    cGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTgg                                                                                                                             <  4:157007/148‑1 (MQ=255)                    cGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGt                                                                                                                            <  7:104427/149‑1 (MQ=255)                     gACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTGTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGtt                                                                                                                           <  3:275621/149‑1 (MQ=255)                      aCAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTg                                                                                                                          <  6:95337/149‑1 (MQ=255)                          aTGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAgg                                                                                                                       >  3:80684/1‑148 (MQ=255)                             caccacGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAAc                                                                                                                    <  6:232876/148‑1 (MQ=255)                             caccacGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAAc                                                                                                                    >  4:842/1‑148 (MQ=255)                                  cgcgCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTTCCCGGCCGTGTCGAACTGGTTGAGGAACGAggtg                                                                                                              <  6:26120/149‑1 (MQ=255)                                  cgcgCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAggtg                                                                                                              >  5:308508/1‑149 (MQ=255)                                     gcAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTc                                                                                                           >  4:286030/1‑149 (MQ=255)                                       aaaTGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGa                                                                                                         <  5:144850/149‑1 (MQ=255)                                        aaTGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAg                                                                                                        <  6:201251/149‑1 (MQ=255)                                         aTGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAgg                                                                                                       >  1:156338/1‑149 (MQ=255)                                          tGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGc                                                                                                      <  4:177778/149‑1 (MQ=255)                                           gCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGATGAGGAACGAGGTGGTCGAGGc                                                                                                      >  5:147646/1‑148 (MQ=255)                                            cTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCa                                                                                                     <  3:161478/148‑1 (MQ=255)                                            cTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATACTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCa                                                                                                     <  5:313095/148‑1 (MQ=255)                                               gCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCaaa                                                                                                   <  6:154186/147‑1 (MQ=255)                                               gCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGACAAAg                                                                                                  >  5:238490/1‑148 (MQ=255)                                                        ggcctcgcgcTCGCGCCCGGTTTTCCGCGCCGCAACGGCCGAATAATCGCCAAACCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGt                                                                                          <  8:47151/141‑1 (MQ=255)                                                        ggCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTg                                                                                         <  3:286030/148‑1 (MQ=255)                                                        ggCATAGCGCTCGCGCAGGGTTTTCAGCCCCGCAACGGCCGAATAATAGACAAAGCTCAAATGCTGGCAATCCAGCGTGACATGGGCCGGTTCCCCGGCCGAGTCGACCTCGTTGAGGAAAGAGGGGGTCGAGGCAAAGAACAGCGTGc                                                                                        >  6:45812/1‑149 (MQ=255)                                                           aTAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCAt                                                                                     >  7:169223/1‑149 (MQ=255)                                                           aTAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCAt                                                                                     <  4:80684/149‑1 (MQ=255)                                                            tAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATg                                                                                    >  4:54484/1‑149 (MQ=255)                                                              gcgcTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGc                                                                                   <  1:310747/148‑1 (MQ=255)                                                              gcgcTCGCGCAGGGGTGTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGGTAACTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCa                                                                                  <  3:245256/149‑1 (MQ=255)                                                                gcTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGt                                                                                >  4:131025/1‑149 (MQ=255)                                                                gcTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGt                                                                                <  3:842/149‑1 (MQ=255)                                                                    gcgcAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGa                                                                             >  5:29603/1‑148 (MQ=255)                                                                      gcAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATaa                                                                          <  1:246129/149‑1 (MQ=255)                                                                      gcAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATaa                                                                          <  6:106608/149‑1 (MQ=255)                                                                        aGGGTATTCAGCGCCGCAACGGCCCAATAATCGACAAAGCTCAAATCCTGGCAATCAAGCGTGACCTGGGCCGGGTCCCCGGCCGAGTCGAAATGGTTGAGGAAACAGGTGGTCGAGGCAAAGAACAGCTTGCCATGCAGATGATAac                                                                         >  8:289135/1‑148 (MQ=255)                                                                         gggTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAacac                                                                       >  3:107656/1‑149 (MQ=255)                                                                                cAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTgccg                                                                 <  3:160427/148‑1 (MQ=255)                                                                                   cgccgcAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCt                                                             >  5:281347/1‑149 (MQ=255)                                                                                   cgccgcAACGGCCGAATAATCGACAAAGCACAAATGCTGGCAATCCAGCGTGCCCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGATGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCt                                                             >  6:59407/1‑149 (MQ=255)                                                                                        cAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCcgcg                                                        <  6:29603/149‑1 (MQ=255)                                                                                           cGGCCGAATAATCGACAAAGCTCAAATGCTGGAAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTCAGGAAAGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTcc                                                     >  7:219410/1‑149 (MQ=255)                                                                                                   taatCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATgg                                               >  1:306600/1‑147 (MQ=255)                                                                                                   taatCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATgg                                               <  2:306600/147‑1 (MQ=255)                                                                                                   taatCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGgt                                             >  6:321041/1‑149 (MQ=255)                                                                                                    aatCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTTGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGgtg                                            >  3:179553/1‑149 (MQ=255)                                                                                                       cGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCg                                         >  2:287766/1‑149 (MQ=255)                                                                                                            aaGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCGCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATa                                    <  6:147646/149‑1 (MQ=255)                                                                                                              gCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACa                                  <  3:54484/149‑1 (MQ=255)                                                                                                                tCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCaaa                                                                                                   >  2:308667/1‑82 (MQ=255)                                                                                                                tCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCaaa                                                                                                   <  1:308667/82‑1 (MQ=255)                                                                                                                  aaaTGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTc                              >  4:94944/1‑149 (MQ=255)                                                                                                                  aaaTGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTc                              >  3:19442/1‑149 (MQ=255)                                                                                                                  aaaTGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTCCCATGCAGTTGATAACACTTCCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTc                              >  1:284447/1‑149 (MQ=255)                                                                                                                          gCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTTGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCa                       <  4:179553/148‑1 (MQ=255)                                                                                                                          gCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCa                       <  3:94944/148‑1 (MQ=255)                                                                                                                          gCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCa                       <  2:156338/148‑1 (MQ=255)                                                                                                                               ccAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCAtgctgc                 >  6:228635/1‑149 (MQ=255)                                                                                                                                  gCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAg              <  3:131025/149‑1 (MQ=255)                                                                                                                                       gacacggCCGGGGACCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGGCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCa           >  2:234482/6‑147 (MQ=255)                                                                                                                                        ccTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAg        >  3:282940/1‑149 (MQ=255)                                                                                                                                        acacggCCGGGGACCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAg        >  5:65673/5‑149 (MQ=255)                                                                                                                                        acacggCCGGGGACCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAg        >  4:189186/5‑149 (MQ=255)                                                                                                                                            gtgctggTCCCCGGCAGTGTCGAACTGGTTGAGGTACTAGGTGGTCGAGGCAAAGACCAGCGTGCCATGCAGTTGATAACACTTGCCGACCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGa    <  5:45812/144‑1 (MQ=255)                                                                                                                                              ccGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGAAAAGTTGAgc  <  6:281347/149‑1 (MQ=255)                                                                                                                                              ccGGGGACCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGAgc  <  6:65673/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               ACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGAGC  >  NC_002947/4647704‑4647993

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