Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,890,792	A→T	11.9%	intergenic (+80/+40)	PP_5164 → / ← plpB	penicillin amidase family protein/NLPA lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,890,792	0	A	T	11.9%	180.5 / 13.0	67	intergenic (+80/+40)	PP_5164/plpB	penicillin amidase family protein/NLPA lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base T (5/3);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTT  >  NC_002947/5890644‑5890936                                                                                                                                                      |                                                                                                                                                 cAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACaa                                                                                                                                                  <  5:59826/149‑1 (MQ=255) cAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACaa                                                                                                                                                  <  4:125195/149‑1 (MQ=255)      tCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGATGCTTTGCGGCCGCTTCCCACCAAAAGGCTTCTACGGTCCTGTGCGGGAAGGGCT‑CATGAAAAACGCCCACACAAGttgt                                                                                                                                             >  5:109764/1‑149 (MQ=255)        aGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTc                                                                                                                                            <  2:133711/148‑1 (MQ=255)             gcgcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGtt                                                                                                                                      >  7:318024/1‑149 (MQ=255)               gcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACCTTGTGTCGGCGt                                                                                                                                       >  4:310132/1‑146 (MQ=255)               gcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACCTTGTGTCGGCGt                                                                                                                                       <  3:310132/146‑1 (MQ=255)               gcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTgg                                                                                                                                    <  8:318024/149‑1 (MQ=255)                 aGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGca                                                                                                                                  >  5:2736/1‑149 (MQ=255)                       gggTGATTCCGGCGCAAAGTACCTTACGGCTGGTAGCCCCGAAGTAAGAGGCTGGTATGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACatat                                                                                                                             >  5:215265/1‑148 (MQ=255)                       gggTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACatat                                                                                                                             >  6:376578/1‑148 (MQ=255)                       gggTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATc                                                                                                                            <  2:17150/149‑1 (MQ=255)                             ttCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGt                                                                                                                      <  6:109744/149‑1 (MQ=255)                                cGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGc                                                                                                                   <  4:304385/149‑1 (MQ=255)                                    gcAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCACGCGTTGCTAGc                                                                                                               >  1:139900/1‑149 (MQ=255)                                              ttACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGaa                                                                                                     >  4:37226/1‑149 (MQ=255)                                              ttACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGaa                                                                                                     >  4:2698/1‑149 (MQ=255)                                                  ggCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGgtgt                                                                                                                                             <  5:176507/105‑1 (MQ=255)                                                  ggCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGgtgt                                                                                                                                             >  6:176507/1‑105 (MQ=255)                                                    cTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGGTTGGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAATAACGCCGAGACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                               <  7:248175/149‑1 (MQ=255)                                                    cTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                               <  3:2698/149‑1 (MQ=255)                                                      ggTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGCCACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTa                                                                                                                 >  2:387432/1‑129 (MQ=255)                                                      ggTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGCCACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTa                                                                                                                 <  1:387432/129‑1 (MQ=255)                                                       gTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAg                                                                                            >  2:188881/1‑149 (MQ=255)                                                        tACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGt                                                                                           <  5:292766/149‑1 (MQ=255)                                                           ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                     >  1:153216/1‑104 (MQ=255)                                                           ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                     >  7:117252/1‑104 (MQ=255)                                                           ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                     <  8:117252/104‑1 (MQ=255)                                                           ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                     <  2:153216/104‑1 (MQ=255)                                                           ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGATGCTAGCGAAATCAGAACGCCGGCACTACg                                                                                        >  2:295735/1‑149 (MQ=255)                                                                     agagGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGc                                                                                                               <  5:37423/116‑1 (MQ=255)                                                                     agagGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGc                                                                                                               >  6:37423/1‑116 (MQ=255)                                                                     agagGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTgg                                                                               >  3:185129/1‑148 (MQ=255)                                                                       agGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTTCGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACTTATCAGGCGTTGCTAGCGAAATCAGATCGCCGGCAGTACGGCACCTTGGTAc                                                                            <  6:215265/149‑1 (MQ=255)                                                                          ctgctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACAAAAGGTGTCGGCGTTGGCACATATCATGCGTTGCAAGCGAAATCAGAACGCCGGCAGTACGGCCCCTTGGTAcgtc                                                                         >  2:402206/1‑149 (MQ=255)                                                                             ctTTGCGGCCGCTTCGCAGCAAAAGGATTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGTTGTCGTCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTActtctt                                                                       <  1:402206/148‑1 (MQ=255)                                                                                      cGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCaatgaat                                                             >  8:24922/1‑149 (MQ=255)                                                                                         ttCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAg                                                                                                                                                 <  8:116870/62‑1 (MQ=255)                                                                                            gctgcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCtt                                                        <  6:53082/145‑1 (MQ=255)                                                                                            gcagcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCtt                                                        <  4:185129/148‑1 (MQ=255)                                                                                            gcagcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCtt                                                        <  3:37226/148‑1 (MQ=255)                                                                                              agcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGAc                                                     <  2:37698/149‑1 (MQ=255)                                                                                                caAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGACGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCCCCTTGGTACTTCTTGGCAATGAATGCCTTGACt                                                    >  1:377847/1‑148 (MQ=255)                                                                                                 aaaaGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACCTTGTGTCGGCGttt                                                                                                                                     >  7:296308/1‑65 (MQ=37)                                                                                                 aaaaGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACCTTGTGTCGGCGttt                                                                                                                                     <  8:296308/66‑2 (MQ=37)                                                                                                 aaaaGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTc                                                  <  1:188881/149‑1 (MQ=255)                                                                                                        ttCTACGGTGCTGTGCGGGCAGGGCTGGAAGAAAAACCCCGACACAAGTTGTCGGCGTTGGCACATATAAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACGTTGGTACTTCTTGGAAATGAATGCATTGACTTCCGGCCt                                            >  1:213169/1‑149 (MQ=255)                                                                                                         tCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCaa                                                                  <  1:261903/125‑1 (MQ=255)                                                                                                         tCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCaa                                                                  >  2:261903/1‑125 (MQ=255)                                                                                                             cGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAg                                      <  5:376578/149‑1 (MQ=255)                                                                                                                   tgtgCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGCCTTCCGGGCTGGTCAGGGCCtt                                >  3:333916/1‑149 (MQ=255)                                                                                                                     tgCGGGCAGGGCT‑GAAGAAAACCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTg                               <  2:318929/148‑1 (MQ=255)                                                                                                                               gCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTg                    >  5:419234/1‑149 (MQ=255)                                                                                                                               gCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTg                    <  6:2736/149‑1 (MQ=255)                                                                                                                                   gaagaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTTCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGAt                 <  6:109764/149‑1 (MQ=255)                                                                                                                                    acctaaaaCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATg                <  3:62468/145‑1 (MQ=255)                                                                                                                                     agaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATgg               >  2:159359/1‑149 (MQ=255)                                                                                                                                         caaCGCCGACACCTTGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCt            >  5:158223/2‑148 (MQ=255)                                                                                                                                         aaaCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGt                                         >  5:71988/1‑119 (MQ=255)                                                                                                                                         aaaCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGt                                         <  6:71988/119‑1 (MQ=255)                                                                                                                                         aaaCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATg                >  4:62468/1‑144 (MQ=255)                                                                                                                                          aaCGCCGACACCTTGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCa          >  6:160300/1‑149 (MQ=255)                                                                                                                                           aCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCAc         <  3:325583/149‑1 (MQ=255)                                                                                                                                           aCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCAc         >  1:162065/1‑149 (MQ=255)                                                                                                                                              ccGACACCTTGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTg                                                                      >  3:233359/1‑85 (MQ=255)                                                                                                                                              ccGACACCTTGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTg                                                                      <  4:233359/85‑1 (MQ=255)                                                                                                                                                gACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTcc                                                 <  5:129727/104‑1 (MQ=255)                                                                                                                                                gACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTcc                                                 >  6:129727/1‑104 (MQ=255)                                                                                                                                                  cacaAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCtt  >  2:130481/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                 CAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTT  >  NC_002947/5890644‑5890936

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑