Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,605,372	A→C	7.7%	Q263P (CAG→CCG) ‡	PP_4077 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,605,372	0	A	C	7.7%	195.4 / 7.0	65	Q263P (CAG→CCG) ‡	PP_4077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base C (2/3);  total (29/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 = 8.12e-01

AGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGG  >  NC_002947/4605227‑4605518                                                                                                                                                   |                                                                                                                                                   aGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAgg                                                                                                                                                  <  3:326878/148‑1 (MQ=255) aGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAgg                                                                                                                                                  <  2:245460/148‑1 (MQ=255)    ccTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTggg                                                                                                                                              >  1:244931/1‑149 (MQ=255)        aggttgCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTCCCCGAAGCGCTTCCACTTATTCACCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACCTGTGGGCag                                                                                                                                           <  8:133223/143‑2 (MQ=255)        aggaggCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTTCCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACCTGTGGGCag                                                                                                                                           <  1:84818/148‑2 (MQ=255)        aggaggCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATc                                                                                                                                          <  1:41909/149‑1 (MQ=255)        aggaggCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATc                                                                                                                                          >  3:45896/1‑149 (MQ=255)         ggaggCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACCTGTGGGCaggc                                                                                                                                         <  5:24167/149‑4 (MQ=255)              ccTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTc                                                                                                                                    <  7:282661/149‑1 (MQ=255)               cTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTc                                                                                                                                    >  1:185780/1‑148 (MQ=255)                      cccTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGtt                                                                                                                                     >  8:320003/1‑140 (MQ=255)                      cccTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGtt                                                                                                                                     <  7:320003/140‑1 (MQ=255)                      cccTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTc                                                                                                                                    >  6:145589/1‑141 (MQ=255)                      cccTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTc                                                                                                                                    <  5:145589/141‑1 (MQ=255)                      cccTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCAc                                                                                                                            >  5:30944/1‑149 (MQ=255)                        cTACCCGCAGAACGTGTTCGATGGCTACCTTATTCTGCAGGAATACTTCGTTTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTTCCCACAGGTGGGCATCGAGTTCCATTTCAcgc                                                                                                                          <  4:313095/149‑1 (MQ=255)                        cTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCAcgc                                                                                                                          >  1:130144/1‑149 (MQ=255)                        cTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCAcgc                                                                                                                          >  5:191383/1‑149 (MQ=255)                           cccGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGt                                                                                                                        <  2:185780/148‑1 (MQ=255)                            ccGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCa                                                                                                                      <  2:261388/149‑1 (MQ=255)                                 gAACGTGTTCGATGGCTACCTTATTCTTCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCTGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGc                                                                                                                 <  2:304080/149‑1 (MQ=255)                                           gATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACTCGTCAACTGCCAGACACCCt                                                                                                       >  4:105929/1‑149 (MQ=255)                                                 tACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGt                                                                                                 >  8:265511/1‑149 (MQ=255)                                                   ccGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCAcgc                                                                                                                          <  5:13671/122‑1 (MQ=255)                                                   ccGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCAcgc                                                                                                                          >  6:13671/1‑122 (MQ=255)                                                   ccGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTg                                                                                               <  4:293417/149‑1 (MQ=255)                                                    cGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCCCGCGTCAACTGCCAGACACCCTTCCCGTTgg                                                                                              >  8:252095/1‑149 (MQ=255)                                                      tATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGca                                                                                            >  2:128471/1‑149 (MQ=255)                                                           tGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACag                                                                                         >  4:9638/1‑147 (MQ=255)                                                           tGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACag                                                                                         >  3:190977/1‑147 (MQ=255)                                                            gCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACag                                                                                         >  7:297080/1‑146 (MQ=255)                                                             cAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACagaa                                                                                       <  6:191383/147‑1 (MQ=255)                                                             cAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACagaa                                                                                       <  2:202149/147‑1 (MQ=255)                                                              aGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGAt                                                                                    >  2:134965/1‑149 (MQ=255)                                                                  aTACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCa                                                                                <  3:9638/149‑1 (MQ=255)                                                                     cTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGc                                                                              <  8:297080/148‑1 (MQ=255)                                                                      ttCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGca                                                                                            <  6:266290/133‑1 (MQ=255)                                                                            ttCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTg                                                                      <  1:134965/149‑1 (MQ=255)                                                                              cccGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTTCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCa                                                                    >  8:239585/1‑149 (MQ=255)                                                                              cccGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCa                                                                    >  5:177776/1‑149 (MQ=255)                                                                                        ttCCACTTCTTCAGCATTACCGGCCTGCGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTg                                                           >  1:15665/1‑149 (MQ=255)                                                                                        ttCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCAc                                                                   <  4:301500/140‑1 (MQ=255)                                                                                        ttCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCAc                                                                   >  3:301500/1‑140 (MQ=255)                                                                                            aCTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGTCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTACCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAc                                                      <  7:92764/149‑1 (MQ=255)                                                                                             cttcttCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAcc                                                     <  6:30944/149‑1 (MQ=255)                                                                                                 atcAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACCAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGt                                                  <  4:193279/147‑1 (MQ=255)                                                                                                   cAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGtt                                                 <  2:130144/147‑1 (MQ=255)                                                                                                     gCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAg                                             <  2:244931/149‑1 (MQ=255)                                                                                                           ccGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGt                                       <  4:45896/149‑1 (MQ=255)                                                                                                            cGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTg                                      <  6:177776/149‑1 (MQ=255)                                                                                                                cTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCgg                                                              <  6:224714/121‑1 (MQ=255)                                                                                                                cTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCgg                                                              >  5:224714/1‑121 (MQ=255)                                                                                                                       aaGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCAt                            >  4:77381/1‑149 (MQ=255)                                                                                                                               gCCGGCTCAAGCCTGCCCACAGGTGGTCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTg                     <  4:110600/148‑1 (MQ=255)                                                                                                                               gCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTg                     <  4:190977/148‑1 (MQ=255)                                                                                                                                    cacaAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGa                                  <  3:41770/130‑1 (MQ=255)                                                                                                                                    cacaAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGa                                  >  4:41770/1‑130 (MQ=255)                                                                                                                                      actcgatGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGc               >  1:138849/7‑147 (MQ=255)                                                                                                                                       aaGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAg            >  2:111397/1‑149 (MQ=255)                                                                                                                                       aaGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAg            >  6:55766/1‑149 (MQ=255)                                                                                                                                        tcgatGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAgg           >  2:214217/5‑149 (MQ=255)                                                                                                                                        aGCCTGCCCACAGGTGGGCATCGAGTTTCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAgg           <  7:239585/149‑1 (MQ=255)                                                                                                                                        aGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAgg           <  7:265511/149‑1 (MQ=255)                                                                                                                                        aGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAgg           <  3:77381/149‑1 (MQ=255)                                                                                                                                           cTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGACACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTg        >  8:111229/1‑149 (MQ=255)                                                                                                                                                cacaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAg   >  7:289400/1‑149 (MQ=255)                                                                                                                                                 acaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAc                                                      <  7:17266/97‑1 (MQ=255)                                                                                                                                                 acaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAc                                                      >  8:17266/1‑97 (MQ=255)                                                                                                                                                 acaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAgg  <  3:230102/149‑1 (MQ=255)                                                                                                                                                 acaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAgg  <  2:15665/149‑1 (MQ=255)                                                                                                                                                 acaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAgg  <  3:105929/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTG‑GAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGG  >  NC_002947/4605227‑4605518

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