Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,032,761	C→G	6.3%	R94G (CGA→GGA)	PP_4435 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,032,761	0	C	G	6.3%	198.8 / 2.1	63	R94G (CGA→GGA)	PP_4435	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/33);  new base G (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

GCGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCACCGCTA  >  NC_002947/5032614‑5032907                                                                                                                                                    |                                                                                                                                                   gCGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTgcgc                                                                                                                                                    >  3:371905/1‑148 (MQ=255) gCGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCACACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTgcgc                                                                                                                                                    <  1:331418/148‑1 (MQ=255)  cGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGa                                                                                                                                                  <  4:20534/149‑1 (MQ=255)  cGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGa                                                                                                                                                  <  3:192183/149‑1 (MQ=255)    aaCATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGc                                                                                                                                                <  1:85288/149‑1 (MQ=255)      cATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCt                                                                                                                                              >  5:141351/1‑149 (MQ=255)          gACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCtgtg                                                                                                                                           <  8:154140/148‑1 (MQ=255)          gACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTgtgt                                                                                                                                          >  6:305688/1‑149 (MQ=255)           aCAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGtt                                                                                                                                         >  8:92365/1‑149 (MQ=255)            cAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTg                                                                                                                                        <  1:161442/149‑1 (MQ=255)             aGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGa                                                                                                                                       <  3:387535/149‑1 (MQ=255)                      gaggagGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCt                                                                                                                                <  6:127797/147‑1 (MQ=255)                      gaggagGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCt                                                                                                                                <  4:167609/147‑1 (MQ=255)                         gaggCAGGTAGTCACTAGTCGTAAATCCTCACAGTTACTATCTATCTTAACGACATCACCTCCGACTTCGCTGCCATGAACGATGTTTGGGAAGAGTGTGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGt                                                                                                                            <  2:187370/148‑1 (MQ=255)                          aggCCGGTATTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTc                                                                                                                          <  7:185773/149‑1 (MQ=255)                            gCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaa                                                                                                                        <  1:341868/149‑1 (MQ=255)                                    tCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaa                                                                                                                        <  6:183869/141‑1 (MQ=255)                                    tCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaa                                                                                                                        >  5:183869/1‑141 (MQ=255)                                          gTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCg                                                                                                          >  5:25549/1‑149 (MQ=255)                                           tCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGt                                                                                                         >  7:384253/1‑149 (MQ=255)                                             gAAAATCCTCACAGTTACTATCTATCTGAACGAACTCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGtt                                                                                                        <  8:100907/148‑1 (MQ=255)                                              aaaaTCCTCACAGTGACTATCTATCTGAAAGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTTGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGtt                                                                                                        >  6:409965/1‑147 (MQ=255)                                                       acaGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTcc                                                                                               <  7:92365/147‑1 (MQ=255)                                                       acaGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTcc                                                                                               <  5:305688/147‑1 (MQ=255)                                                        caGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCgta                                                                                            <  6:263427/149‑1 (MQ=255)                                                        caGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCgta                                                                                            <  6:25549/149‑1 (MQ=255)                                                              ctatctatCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTagca                                                                                      <  4:249891/149‑1 (MQ=255)                                                              ctatctatCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTagca                                                                                      <  6:242605/149‑1 (MQ=255)                                                                atctatctGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTagcag                                                                                     >  3:51139/1‑148 (MQ=255)                                                                       tGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTg                                                                                                       >  3:342346/1‑123 (MQ=255)                                                                       tGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTg                                                                                                       <  4:342346/123‑1 (MQ=255)                                                                         aaCGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCt                                                                            <  5:409965/148‑1 (MQ=255)                                                                         aaCGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCt                                                                            >  3:212139/1‑148 (MQ=255)                                                                                 caccacCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTa                                                                    >  6:29809/1‑148 (MQ=255)                                                                                   ccaccGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGc                                                                              <  1:228953/136‑1 (MQ=255)                                                                                   ccaccGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGc                                                                              >  2:228953/1‑136 (MQ=255)                                                                                     accGACTTCGCTGCCATGAACGAGGGTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGACGCGGCTATGTTCAAGCCAACACTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCATAAAAgc                                                               >  5:139508/1‑149 (MQ=255)                                                                                      ccGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAgcg                                                              >  2:188093/1‑149 (MQ=255)                                                                                       cGACTTCGCTGCCATGAACGATGTTTGGGAAGAGTGGGTAGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGAAATCGCCTGACGCTTAAAAgag                                                              >  5:58672/1‑146 (MQ=255)                                                                                           ttCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCAccc                                                         >  1:381226/1‑149 (MQ=255)                                                                                                    aTGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGgcgc                                                <  6:141351/149‑1 (MQ=255)                                                                                                         ggtggTTTGGGAAGAGTGTGTCGCGGAAGGTCAGGCTCCTGCGCGAGCGTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGTGGCTTAAAAGCGCACCCGGCAGGCGCCTAtg                                           <  6:139508/146‑2 (MQ=255)                                                                                                             gTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGc                                       >  8:222251/1‑149 (MQ=255)                                                                                                             cttggggAAGAGTGGGTCGCGGCAGGTCAGGCGCCTGCGCGATCCGGTGGTGAGGCGGCTATGTTCACGCCAAGCCTGCTCGGTGAGATGTCCGTCGTAGCAGCAATCGCCTGACGCGTAAAAGCGCACCCGGCAGGCGCCTATCCAGc                                       <  4:295889/145‑1 (MQ=255)                                                                                                              tttGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCgt                                                                                             <  8:61270/94‑1 (MQ=255)                                                                                                              tttGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCgt                                                                                             >  7:61270/1‑94 (MQ=255)                                                                                                                 gggAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                                    <  5:29809/148‑1 (MQ=255)                                                                                                                  ggAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                                    <  3:272985/147‑1 (MQ=255)                                                                                                                  ggAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                                    >  8:41752/1‑147 (MQ=255)                                                                                                                  ggAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                                    >  4:272985/1‑147 (MQ=255)                                                                                                                  ggAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                                    <  7:41752/147‑1 (MQ=255)                                                                                                                      gagTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACa                              <  7:222251/149‑1 (MQ=255)                                                                                                                       agTGGGTCTCGGAAGGTCAGTCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTTAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAg                             <  3:37017/149‑1 (MQ=255)                                                                                                                       agTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAg                             <  8:121442/149‑1 (MQ=255)                                                                                                                            gTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGt                           >  6:228682/1‑146 (MQ=255)                                                                                                                             tCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCTGGCAGGCGCCTATCCAGCCTGGTAACAGGTCa                         >  6:171623/1‑147 (MQ=255)                                                                                                                               gcgGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAgg                     <  4:51139/149‑1 (MQ=255)                                                                                                                                  tcaacacAGGCTCGCGCAGGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGgcgc                  >  5:48881/7‑149 (MQ=255)                                                                                                                                  gAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGgcgc                  >  7:37946/1‑149 (MQ=255)                                                                                                                                     acacAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCg               >  2:132799/4‑149 (MQ=255)                                                                                                                                      cacAGGCTCGCGCAGGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCagag                      >  7:271871/3‑141 (MQ=255)                                                                                                                                      cacAGGCTCGCGCAGGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCagag                      <  8:271871/139‑1 (MQ=255)                                                                                                                                       tCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGcc             >  4:601/1‑149 (MQ=255)                                                                                                                                       acAGGCTCGCGCAGGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGcc             <  6:48881/148‑1 (MQ=255)                                                                                                                                             tCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCAc       >  8:329674/1‑149 (MQ=255)                                                                                                                                                   gcgAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCACCGCTa  >  2:33326/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCACCGCTA  >  NC_002947/5032614‑5032907

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