Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,661,887	G→C	6.0%	L504L (CTG→CTC) ‡	ligB →	DNA ligase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,661,887	0	G	C	6.0%	245.2 / 2.2	76	L504L (CTG→CTC) ‡	ligB	DNA ligase B
Reads supporting (aligned to +/- strand):  ref base G (42/29);  new base C (2/3);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCGATGCCGAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGTTGGCCGAAGCCG  >  NC_002947/5661740‑5662032                                                                                                                                                    |                                                                                                                                                   gCTTGATGCCGAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCactga                                                                                                                                                   <  1:16895/149‑1 (MQ=255) gCTCGATGCCGAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCactga                                                                                                                                                   <  1:308450/149‑1 (MQ=255)    cGATGCCGAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCactgactg                                                                                                                                                >  2:104158/1‑149 (MQ=255)      aTGCCGAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGactgg                                                                                                                                               >  1:119838/1‑148 (MQ=255)          cgcgcgcCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGactggcaact                                                                                                                                          >  5:238899/4‑149 (MQ=255)          cGAGCGCCTGCAGCAGGCCGCCGGAATCAGGAGACTCTCTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGACACTGactggcaact                                                                                                                                          <  7:61979/149‑1 (MQ=255)           gAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGactggcaactg                                                                                                                                         >  5:100141/1‑149 (MQ=255)           gAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGactggcaactg                                                                                                                                         >  2:68516/1‑149 (MQ=255)             gcgcCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCAGTCagtggcaaatgct                                                                                                                                       <  6:234756/149‑3 (MQ=255)              cgcCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCAGTCagtggcaaatgctt                                                                                                                                      <  4:165524/149‑4 (MQ=255)                 cTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCactgactg                                                                                                                                                >  6:73662/1‑136 (MQ=255)                 cTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCactgactg                                                                                                                                                <  5:73662/136‑1 (MQ=255)                  tGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGgcg                                                                                                                                  >  4:140078/1‑149 (MQ=255)                   gcagcaGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGgcgc                                                                                                                                 >  3:342316/1‑149 (MQ=255)                      gcagGCCCCCGGAATCACCAGACTCACTGCAACACAGATGCTGCGCAGCTTCGACCCAGCCCGCTCACGCCCGTGCGCTCCGCGGCCGCGCGCCCTTGGCGTGCCCACTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTgg                                                                                                                              <  4:71293/149‑1 (MQ=255)                            ccccGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCAc                                                                                                                        >  2:45655/1‑149 (MQ=255)                            ccccGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCAc                                                                                                                        >  1:71363/1‑149 (MQ=255)                              ccGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACg                                                                                                                       >  2:68284/1‑148 (MQ=255)                              ccGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTAGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACg                                                                                                                       >  7:178781/1‑148 (MQ=255)                                     cagcagACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCg                                                                                                               >  3:68477/1‑149 (MQ=255)                                      agcagACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCgg                                                                                                              >  4:39664/1‑149 (MQ=255)                                           aCTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAAt                                                                                                         <  6:100141/149‑1 (MQ=255)                                             taaatGCAACACAGATGATGGGCAGTTTCGACCAAGCCAGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGACCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCTCTGGCTTCACGCAGTGCCGGGCAATgg                                                                                                       >  4:24971/5‑149 (MQ=255)                                             tCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATgg                                                                                                       >  6:160190/1‑149 (MQ=255)                                                 tGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCaaa                                                                                                   >  7:67954/1‑149 (MQ=255)                                                  gCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAAt                                                                                                         >  2:311036/1‑142 (MQ=255)                                                  gCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTACGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAAt                                                                                                         <  1:311036/142‑1 (MQ=255)                                                   cAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCaaa                                                                                                   <  5:219568/147‑1 (MQ=255)                                                   cAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAAcc                                                                                                 <  3:39664/149‑1 (MQ=255)                                                   cAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAAcc                                                                                                 >  5:13545/1‑149 (MQ=255)                                                     acacAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCGTGGCGTGCCCATTGGCAAGCATTTGCCCCGGACGGGCAACGGGCAGGCGCTGGCTGCACGCAGTGCCGGGCAAGGGCAAACCgg                                                                                               >  2:341300/1‑148 (MQ=255)                                                      cacaGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTg                                                                                              <  3:167623/149‑1 (MQ=255)                                                         aGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑ccc                                                                                          >  5:131782/1‑149 (MQ=255)                                                             gCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCa                                                                                      >  7:195736/1‑149 (MQ=255)                                                                  gCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGgc                                                                                 <  3:140078/149‑1 (MQ=255)                                                                         cGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCaaaa                                                                           <  2:71363/148‑1 (MQ=255)                                                                             cAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTc                                                                      <  1:68284/149‑1 (MQ=255)                                                                                 cccGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCgcgc                                                                  <  7:3619/149‑1 (MQ=255)                                                                                  ccGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCgcgc                                                                  >  4:102974/1‑148 (MQ=255)                                                                                     cTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAAc                                                              >  2:99779/1‑149 (MQ=255)                                                                                         cGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTtgt                                                          <  8:195736/149‑1 (MQ=255)                                                                                          gCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTGCCCCGGCATCGGCGCAAAACGCTCGCGCCAACTtgt                                                          <  1:45655/149‑1 (MQ=255)                                                                                          gCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTtgtg                                                         >  6:242799/1‑149 (MQ=255)                                                                                          gCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTtgtg                                                         >  2:349953/1‑149 (MQ=255)                                                                                          gCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTtgtg                                                         >  4:58634/1‑149 (MQ=255)                                                                                          gCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTtgtg                                                         >  8:43071/1‑149 (MQ=255)                                                                                          gCCCGTTCGATCAGTGGTTGCGCGCCCTTAGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTtgtg                                                         <  8:178781/149‑1 (MQ=255)                                                                                            ccGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTgg                                                        >  6:194933/1‑148 (MQ=255)                                                                                               ttCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAAtt                                                    <  8:67954/149‑1 (MQ=255)                                                                                               ttCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAAtt                                                    >  5:164370/1‑149 (MQ=255)                                                                                                  gATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGt                                                                                                                   >  6:121101/1‑84 (MQ=255)                                                                                                  gATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGt                                                                                                                   <  5:121101/84‑1 (MQ=255)                                                                                                  gATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAAttttt                                                 >  2:12338/1‑149 (MQ=255)                                                                                                      aGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCTGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTcgc                                              <  7:255069/148‑1 (MQ=255)                                                                                                      aGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTcgcg                                             <  1:99779/149‑1 (MQ=255)                                                                                                      aGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTcgc                                              <  4:342316/148‑1 (MQ=255)                                                                                                      aGTGGGTGCGCGCCCTTGGCGTGCCCATTGTCAAGCTTTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGGAAACGGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGTAATTTTTcgc                                              <  1:205464/148‑1 (MQ=255)                                                                                                       gTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCgg                                            >  5:98521/1‑149 (MQ=255)                                                                                                          gTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCt                                         >  6:257657/1‑149 (MQ=255)                                                                                                             gcgcgcCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTcc                                       <  1:372110/148‑1 (MQ=255)                                                                                                               gcgcCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAg                                    >  5:191753/1‑149 (MQ=255)                                                                                                               gcgcCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAg                                    <  1:349953/149‑1 (MQ=255)                                                                                                               gcgcCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAg                                    <  3:24971/149‑1 (MQ=255)                                                                                                               gcgcCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAg                                    <  3:58634/149‑1 (MQ=255)                                                                                                                   ccTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTg                                 >  6:235130/1‑148 (MQ=255)                                                                                                                   ccTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCACTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTg                                 >  8:34103/1‑148 (MQ=255)                                                                                                                    cTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTACCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGc                                >  3:337089/1‑148 (MQ=255)                                                                                                                       ggCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCaa                              <  6:13545/147‑1 (MQ=255)                                                                                                                       ggCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAg                             <  3:358864/148‑1 (MQ=255)                                                                                                                       ggCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCCCCGAGGTGCaa                              <  1:341300/147‑1 (MQ=255)                                                                                                                                aTTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGc                   <  6:164370/149‑1 (MQ=255)                                                                                                                                 ttGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGca                  <  5:235130/149‑1 (MQ=255)                                                                                                                                     ctgccagttGCCAGTCAGTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAACCGCTCGCGCCACCTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGt              >  6:364070/8‑149 (MQ=255)                                                                                                                                     cAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGt              >  6:140044/1‑149 (MQ=255)                                                                                                                                      aaGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCTGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGATCGCGCCAACTTGTGGATTTTTTAGCGGACTCCTAGGTGCAATACATTGCAGCAAAGtt             >  8:352287/1‑149 (MQ=255)                                                                                                                                      aaGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGtt             >  7:217076/1‑149 (MQ=255)                                                                                                                                      aaGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGtt             >  6:126251/1‑149 (MQ=255)                                                                                                                                        gCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGTTgg           >  2:204792/1‑149 (MQ=255)                                                                                                                                         cagttGCCAGTCAGTGGCAAATGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAg                                    <  6:350055/120‑1 (MQ=255)                                                                                                                                         cagttGCCAGTCAGTGGCAAATGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAg                                    >  5:350055/4‑123 (MQ=255)                                                                                                                                          aTTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGTTGGcc         >  5:140021/1‑149 (MQ=255)                                                                                                                                                 actgactgGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGTTGGCCGAAGCCg  >  4:30500/1‑149 (MQ=255)                                                                                                                                                 actgactgGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGTTGGCCGAAGCCg  <  6:131782/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCTCGATGCCGAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGTTGGCCGAAGCCG  >  NC_002947/5661740‑5662032

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