Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,890,793	G→T	12.0%	intergenic (+81/+39)	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,793	0	G	T	12.0%	182.5 / 13.8	66	intergenic (+81/+39)	PP_5164/plpB	penicillin amidase family protein/NLPA lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base T (5/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

TCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTT  >  NC_002947/5890649‑5890936                                                                                                                                                  |                                                                                                                                                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                                                                                                                                             >  6:176507/1‑105 (MQ=255)                                             ggCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGgtgt                                                                                                                                             <  5:176507/105‑1 (MQ=255)                                               cTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGGTTGGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAATAACGCCGAGACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                               <  7:248175/149‑1 (MQ=255)                                               cTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                               <  3:2698/149‑1 (MQ=255)                                                 ggTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGCCACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTa                                                                                                                 <  1:387432/129‑1 (MQ=255)                                                 ggTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGCCACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTa                                                                                                                 >  2:387432/1‑129 (MQ=255)                                                  gTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAg                                                                                            >  2:188881/1‑149 (MQ=255)                                                   tACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGt                                                                                           <  5:292766/149‑1 (MQ=255)                                                      ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                     <  8:117252/104‑1 (MQ=255)                                                      ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                     >  1:153216/1‑104 (MQ=255)                                                      ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                     >  7:117252/1‑104 (MQ=255)                                                      ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                     <  2:153216/104‑1 (MQ=255)                                                      ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGATGCTAGCGAAATCAGAACGCCGGCACTACg                                                                                        >  2:295735/1‑149 (MQ=255)                                                                agagGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGc                                                                                                               >  6:37423/1‑116 (MQ=255)                                                                agagGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGc                                                                                                               <  5:37423/116‑1 (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                                                        <  3:37226/148‑1 (MQ=255)                                                                                       gcagcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCtt                                                        <  4:185129/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                                                                  >  2:261903/1‑125 (MQ=255)                                                                                                    tCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCaa                                                                  <  1:261903/125‑1 (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                    <  6:2736/149‑1 (MQ=255)                                                                                                                          gCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTg                    >  5:419234/1‑149 (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                                                                      <  4:233359/85‑1 (MQ=255)                                                                                                                                         ccGACACCTTGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTg                                                                      >  3:233359/1‑85 (MQ=255)                                                                                                                                           gACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTcc                                                 <  5:129727/104‑1 (MQ=255)                                                                                                                                           gACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTcc                                                 >  6:129727/1‑104 (MQ=255)                                                                                                                                             cacaAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCtt  >  2:130481/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                TCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCT‑GAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTT  >  NC_002947/5890649‑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: ‑