Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,150,406	A→T	6.9%	S354T (TCG→ACG)	PP_4535 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,150,406	0	A	T	6.9%	216.3 / 3.4	72	S354T (TCG→ACG)	PP_4535	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base T (2/3);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AATTCCTTGCGAAGCTTTTTGAGATCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAACGGGTGTTGAGCGCTTTCTGAATGCTG  >  NC_002947/5150258‑5150551                                                                                                                                                      |                                                                                                                                                   aaTTCCTTGCGAAGCTTTTTGAGATCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGa                                                                                                                                                    >  3:175663/1‑149 (MQ=255)   ttCCTTGCGAAGCTTTTTGAGATCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGAtt                                                                                                                                                  <  5:91698/149‑1 (MQ=255)       ttGCGAAGCTTTTTGAGATCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGc                                                                                                                                                >  2:279563/1‑147 (MQ=255)          cGAAGCTTTTTGAGATCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATTCGCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAg                                                                                                                                             <  1:279563/148‑1 (MQ=255)          cGAAGCTTTTTGACTTCAGGGTCCTCCGCCCGGCGCTCCTCGCCAAAGCCTTGACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACTGTTTGTTCGCTTCATCCTGAGCCACATCCGATTGCGAgg                                                                                                                                            <  2:323005/148‑1 (MQ=255)            aaGCTTTTTGAGATCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGTTTGCgcctcg                                                                                                                                          >  8:351437/1‑143 (MQ=255)                tttttGAGATCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAAc                                                                                                                                    >  2:20522/1‑149 (MQ=255)                   ttGAGATCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGTTTGCGCCTCGC‑AATCgga                                                                                                                                 >  4:3925/1‑148 (MQ=255)                   ttGAGATCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACgg                                                                                                                                  <  1:341932/148‑1 (MQ=255)                     gagaTCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGTTTGCGCCTCGC‑AATCggat                                                                                                                                <  7:236566/148‑3 (MQ=255)                     gagaTCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGTTTGCGCCTCGC‑AATCggat                                                                                                                                <  7:351437/148‑3 (MQ=255)                     gagaTCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACgggg                                                                                                                                <  7:16347/148‑1 (MQ=255)                     gagaTCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACgggg                                                                                                                                >  7:202294/1‑148 (MQ=255)                       gaTCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACgggg                                                                                                                                >  4:192579/1‑146 (MQ=255)                        aTCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGTTTGCGCCTCGC‑AATCggatgtgg                                                                                                                            <  3:3925/149‑7 (MQ=255)                        aTCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGgcg                                                                                                                                          >  7:153508/1‑136 (MQ=255)                        aTCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGa                                                                                                                            >  1:227515/1‑149 (MQ=255)                        aTCAGCGTCCACCGCAGGGCTCGCATCGCCAAAGCCTTCAGGCGCGTAAGTGGAAACTTTGAGCTCGTCTTTAAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTTCGAGgcg                                                                                                                                          <  8:153508/136‑1 (MQ=255)                            gCGTCCACGGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACcgc                                                                                                                        <  6:323757/149‑1 (MQ=255)                            gCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACcgc                                                                                                                        >  8:337355/1‑149 (MQ=255)                                    cGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCaa                                                                                                                <  7:204908/149‑1 (MQ=255)                                     gCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAg                                                                                                               >  5:301137/1‑149 (MQ=255)                                     gCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAg                                                                                                               >  2:330742/1‑149 (MQ=255)                                      cccGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAgg                                                                                                              >  6:244838/1‑149 (MQ=255)                                       ccGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGt                                                                                                             >  5:289331/1‑149 (MQ=255)                                         ggCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTg                                                                                                           >  6:143726/1‑149 (MQ=255)                                         ggCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTg                                                                                                           <  5:244838/149‑1 (MQ=255)                                          gcgcTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGtt                                                                                                            >  1:104549/1‑147 (MQ=255)                                          gcgcTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGtt                                                                                                            <  2:104549/147‑1 (MQ=255)                                           cgcTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTgg                                                                                                          >  8:163756/1‑148 (MQ=255)                                            gcTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCa                                                                                                        <  2:227515/149‑1 (MQ=255)                                                   gCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCg                                                                                                 <  6:289331/149‑1 (MQ=255)                                                   gCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCg                                                                                                 <  7:163756/149‑1 (MQ=255)                                                   gCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCg                                                                                                 <  7:337355/149‑1 (MQ=255)                                                    ccAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCgg                                                                                                <  6:301137/149‑1 (MQ=255)                                                    ccAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCgg                                                                                                <  5:22312/149‑1 (MQ=255)                                                    ccAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCgg                                                                                                >  4:331503/1‑149 (MQ=255)                                                     cAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCggg                                                                                               <  3:192579/149‑1 (MQ=255)                                                     cAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCgg                                                                                                >  2:294674/1‑148 (MQ=255)                                                      aaaGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGa                                                                                              >  8:224546/1‑149 (MQ=255)                                                            ttCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCGGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAgg                                                                                          >  3:358982/1‑147 (MQ=255)                                                                cgcgcgTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTgg                                                                                                          >  1:298138/1‑127 (MQ=255)                                                                cgcgcgTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTgg                                                                                                          <  2:298138/127‑1 (MQ=255)                                                                      tAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGTGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGAc                                                                              <  1:73329/149‑1 (MQ=255)                                                                        aGTGGAAACTTTGAGCTCGTCGGTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGCAAAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGAc                                                                              <  1:294674/148‑1 (MQ=255)                                                                           ggAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGt                                                                         >  2:347705/1‑149 (MQ=255)                                                                            gAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGtt                                                                        >  7:226904/1‑149 (MQ=255)                                                                               aCTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGc                                                                     <  1:20522/149‑1 (MQ=255)                                                                                   tGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAg                                                                 >  5:302651/1‑149 (MQ=255)                                                                                     aGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGc                                                                <  5:143726/148‑1 (MQ=255)                                                                                     aGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGc                                                                <  2:86840/148‑1 (MQ=255)                                                                                      gCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGGTTCTTCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCg                                                              >  5:172540/1‑149 (MQ=255)                                                                                      gCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCg                                                              <  8:202294/149‑1 (MQ=255)                                                                                       cTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACgg                                                                                                                                  <  6:214681/80‑1 (MQ=255)                                                                                       cTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACgg                                                                                                                                  >  5:214681/1‑80 (MQ=255)                                                                                           tcgtTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGc                                                         >  4:319996/1‑149 (MQ=255)                                                                                                    cAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCGTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGt                                                >  3:298644/1‑149 (MQ=255)                                                                                                    cAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGt                                                <  1:347705/149‑1 (MQ=255)                                                                                                      gTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCGTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTgg                                              <  4:298644/149‑1 (MQ=255)                                                                                                      gTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATGCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGTTAGCCAGCGCGGTGACCGAGTTTGCAAAGCCGGTCGCCTGTTCCGTgg                                              <  6:172540/149‑1 (MQ=255)                                                                                                         gATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGgaga                                           <  8:226904/149‑1 (MQ=255)                                                                                                         gATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGgaga                                           >  1:71295/1‑149 (MQ=255)                                                                                                                  ccGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTcag                                    >  6:173486/1‑148 (MQ=255)                                                                                                                   cGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTcag                                    >  1:207765/1‑147 (MQ=255)                                                                                                                         tttGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAAc                            <  3:331503/149‑1 (MQ=255)                                                                                                                           tGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAACg                           <  4:358982/148‑1 (MQ=255)                                                                                                                               cGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAACgg                          <  1:151740/145‑1 (MQ=255)                                                                                                                                    cATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAACGGGTGTTGa                   <  7:224546/147‑1 (MQ=255)                                                                                                                                        cTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAACGGGTGTTGAGCGCtt             >  8:168385/1‑149 (MQ=255)                                                                                                                                         tCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAACGGGTGTTGAGCGCttt            >  1:337417/1‑149 (MQ=255)                                                                                                                                             ccACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAACGGGTGTTGAGCGCTTTCTg         >  4:213211/1‑148 (MQ=255)                                                                                                                                             ccACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAACGGGTGTTGAGCGCTTTCTg         >  3:48583/1‑148 (MQ=255)                                                                                                                                                    cGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGTATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAACGGGTGTTGAGCGCTTTCTGAATGCTg  <  7:168385/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   AATTCCTTGCGAAGCTTTTTGAGATCAGCGTCCACCGCCCGGCGCTCCTCGCCAAAGCCTTCACGCGCGTAAGTGGAAACTTTGAGCTCGTCGTTCAGCCAGTGGATT‑GCTGCCGACAGTTTGTTCGCTTCATCCTCAGCCACATCCGATTGCGAGGCGC‑AAACGGGGCAGACCGCGTGTTCAAGGTTGGCAGATTTCGGGATAGGCAGCGCGGTGACCGAGTTGGCAAAGCCGGTCGCCTGTTCGGTGGAGACTTTCAGCAGATAACGGGTGTTGAGCGCTTTCTGAATGCTG  >  NC_002947/5150258‑5150551

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