Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,661,885	C→G	6.2%	L504V (CTG→GTG) ‡	ligB →	DNA ligase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,661,885	0	C	G	6.2%	241.4 / 4.0	76	L504V (CTG→GTG) ‡	ligB	DNA ligase B
Reads supporting (aligned to +/- strand):  ref base C (42/29);  new base G (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 = 7.73e-01

GCAGCTCGATGCCGAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGTTGGCCGAAGCCG  >  NC_002947/5661737‑5662032                                                                                                                                                     |                                                                                                                                                     gCAGCTCGATGCCGAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCac                                                                                                                                                      >  2:300364/1‑149 (MQ=255)    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                                                                                                                                                <  5:73662/136‑1 (MQ=255)                    cTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCactgactg                                                                                                                                                >  6:73662/1‑136 (MQ=255)                     tGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGgcg                                                                                                                                  >  4:140078/1‑149 (MQ=255)                      gcagcaGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGgcgc                                                                                                                                 >  3:342316/1‑149 (MQ=255)                         gcagGCCCCCGGAATCACCAGACTCACTGCAACACAGATGCTGCGCAGCTTCGACCCAGCCCGCTCACGCCCGTGCGCTCCGCGGCCGCGCGCCCTTGGCGTGCCCACTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTgg                                                                                                                              <  4:71293/149‑1 (MQ=255)                               ccccGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCAc                                                                                                                        >  1:71363/1‑149 (MQ=255)                               ccccGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCAc                                                                                                                        >  2:45655/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                                                         >  2:349953/1‑149 (MQ=255)                                                                                             gCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTtgtg                                                         >  4:58634/1‑149 (MQ=255)                                                                                             gCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTtgtg                                                         >  8:43071/1‑149 (MQ=255)                                                                                             gCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTtgtg                                                         >  6:242799/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)                                                                                                                                                     |                                                                                                                                                     GCAGCTCGATGCCGAGCGCCTGCAGCAGGCCCCCGGAATCAGCAGACTCACTGCAACACAGATGCTGGGCAGTTTCGACCAAGCCCGCTCACGCCCGTTCGATCAGTGGTTGCGCGCCCTTGGCGTGCCCATTGGCAAGCATTTGCCACTGACTGGCAACTGGCAGGCGCTGGCTTCACGCAGTGCCGGGCAATGGCAAACCGTG‑CCCGGCATCGGCGCAAAACGCTCGCGCCAACTTGTGGAATTTTTCGCGGCCTCCGAGGTGCAAGCCATTGCCGCACAGTTGGCCGAAGCCG  >  NC_002947/5661737‑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: ‑