Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,032,757	T→C	5.9%	P92P (CCT→CCC)	PP_4435 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,032,757	0	T	C	5.9%	197.1 / 2.1	67	P92P (CCT→CCC)	PP_4435	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base C (2/2);  total (31/36)
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.34e-01

CCCTGGCGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCAC  >  NC_002947/5032609‑5032902                                                                                                                                                     |                                                                                                                                                  cccTGGCGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCt                                                                                                                                                   <  2:129283/149‑1 (MQ=255)  ccTGGCGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATATATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTg                                                                                                                                                  >  4:37017/1‑149 (MQ=255)    tGGCGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTgcg                                                                                                                                                >  4:192183/1‑149 (MQ=255)      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                                                                                                  <  4:342346/123‑1 (MQ=255)                                                                            tGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTg                                                                                                  >  3:342346/1‑123 (MQ=255)                                                                              aaCGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCt                                                                       >  3:212139/1‑148 (MQ=255)                                                                              aaCGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCt                                                                       <  5:409965/148‑1 (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                                                                                        >  7:61270/1‑94 (MQ=255)                                                                                                                   tttGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCgt                                                                                        <  8:61270/94‑1 (MQ=255)                                                                                                                      gggAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                               <  5:29809/148‑1 (MQ=255)                                                                                                                       ggAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                               <  3:272985/147‑1 (MQ=255)                                                                                                                       ggAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                               <  7:41752/147‑1 (MQ=255)                                                                                                                       ggAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                               >  8:41752/1‑147 (MQ=255)                                                                                                                       ggAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                               >  4:272985/1‑147 (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                 <  8:271871/139‑1 (MQ=255)                                                                                                                                           cacAGGCTCGCGCAGGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCagag                 >  7:271871/3‑141 (MQ=255)                                                                                                                                            tCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGcc        >  4:601/1‑149 (MQ=255)                                                                                                                                            acAGGCTCGCGCAGGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGcc        <  6:48881/148‑1 (MQ=255)                                                                                                                                                  tCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCAc  >  8:329674/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCCTGGCGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCAC  >  NC_002947/5032609‑5032902

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