Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,961,628	T→A	9.1%	I175N (ATT→AAT)	lysA‑II →	diaminopimelate decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,961,628	0	T	A	9.1%	185.0 / 7.7	66	I175N (ATT→AAT)	lysA‑II	diaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (3/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

CACCGACGAGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGCCT  >  NC_002947/5961482‑5961775                                                                                                                                                   |                                                                                                                                                     cACCGACGAGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCAtt                                                                                                                                                     <  1:221023/148‑1 (MQ=255)           cTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgccgatgccga                                                                                                                                          <  3:301778/149‑1 (MQ=255)              gAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGGCAAtggcgatgccgaac                                                                                                                                        <  8:86251/148‑3 (MQ=255)              gAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCACTGCGGGTCAACCCGGACGTCGATGCCGGTACACACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGc                                                                                                                                       >  1:90173/1‑149 (MQ=255)                    cTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCta                                                                                                                                 >  6:286956/1‑149 (MQ=255)                             gtCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACgcgc                                                                                                                        >  8:196525/1‑149 (MQ=255)                              tCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgcc                                                                                                                       <  2:136020/149‑1 (MQ=255)                               cgttttcgAAATTGGCAATGTTTCCCCGGTGTCGCTTCGGGTCAACCCGGACGTCGATTCCGGTACCCACCCGTACATCACCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgcc                                                                                                                       <  3:266338/142‑1 (MQ=255)                                gccgccGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgccg                                                                                                                      >  8:150923/1‑148 (MQ=255)                                  cgccgAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTAACCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGATGCCATCTACTTCAGCgcagca                                                                                                                    >  3:266487/1‑148 (MQ=255)                                        aaTGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTAc                                                                                                             >  7:288537/1‑149 (MQ=255)                                           gggCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgccgatgcc                                                                                                                                            <  6:321171/115‑1 (MQ=255)                                           gggCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgccgatgcc                                                                                                                                            >  5:321171/1‑115 (MQ=255)                                             gCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAAc                                                                                                        >  7:91145/1‑149 (MQ=255)                                                 ggTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTg                                                                                                     >  5:7391/1‑148 (MQ=255)                                                   ttGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                                  <  4:77469/149‑1 (MQ=255)                                                   atGCCCCGGTGGCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTAACAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                                  <  8:322879/148‑1 (MQ=255)                                                       cccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTgg                                                                                              >  4:227746/1‑149 (MQ=255)                                                       cccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTgg                                                                                              <  6:7391/149‑1 (MQ=255)                                                       cccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTgg                                                                                              >  6:177251/1‑149 (MQ=255)                                                        ccGGTGTCGCTGCGGGTCAACCCGGACGTCGGTGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGt                                                                                             <  7:239224/149‑1 (MQ=255)                                                        ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCt                                                                                                                                  <  5:187951/112‑1 (MQ=255)                                                        ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCt                                                                                                                                  >  6:187951/1‑112 (MQ=255)                                                        ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGt                                                                                             <  8:91145/149‑1 (MQ=255)                                                        ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGt                                                                                             <  3:227746/149‑1 (MQ=255)                                                          ggTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCg                                                                                           <  7:150923/149‑1 (MQ=255)                                                              tCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGgtg                                                                                        <  8:46568/148‑1 (MQ=255)                                                              tCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACATGGAAGTGGTCGgtgt                                                                                       >  4:56280/1‑149 (MQ=255)                                                                             ccGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCgg                                                                        <  4:16556/149‑1 (MQ=255)                                                                              cGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGc                                                                       <  3:56280/149‑1 (MQ=255)                                                                                        tGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑Ga                                                            <  7:196525/149‑1 (MQ=255)                                                                                          ccGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑Gac                                                           <  5:177251/148‑1 (MQ=255)                                                                                                cccacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTc                                                    >  6:164939/1‑149 (MQ=255)                                                                                                 ccacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTc                                                    >  4:180065/1‑148 (MQ=255)                                                                                                  cacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTc                                                    >  3:241627/1‑147 (MQ=255)                                                                                                  cacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTc                                                    >  5:128874/1‑147 (MQ=255)                                                                                                   acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGcc                                                                               >  8:306981/1‑120 (MQ=255)                                                                                                   acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGcc                                                                               <  7:306981/120‑1 (MQ=255)                                                                                                   acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAg                                                 <  3:180065/149‑1 (MQ=255)                                                                                                   acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAg                                                 >  2:199558/1‑149 (MQ=255)                                                                                                    cccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGc                                                >  6:270281/1‑149 (MQ=255)                                                                                                    cccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGc                                                <  5:286956/149‑1 (MQ=255)                                                                                                     ccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATc                                                                          <  8:132340/123‑1 (MQ=255)                                                                                                     ccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATc                                                                          >  7:132340/1‑123 (MQ=255)                                                                                                              tcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGa                                      >  5:342464/1‑149 (MQ=255)                                                                                                              tcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGa                                      >  5:18880/1‑149 (MQ=255)                                                                                                                    ggCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCt                                <  4:241627/149‑1 (MQ=255)                                                                                                                     gCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTcc                                          >  5:3448/1‑138 (MQ=255)                                                                                                                     gCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTcc                                          <  6:3448/138‑1 (MQ=255)                                                                                                                     gCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTc                               >  3:110131/1‑149 (MQ=255)                                                                                                                       cTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGa                             <  2:90173/149‑1 (MQ=255)                                                                                                                        ttAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGAc                                     <  5:345511/140‑1 (MQ=255)                                                                                                                        ttAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGACCCATCGGCTAGCAGCT‑GACCACCGTCGAGCCCTTCCTCGAc                                     >  6:345511/1‑140 (MQ=255)                                                                                                                          aaaGAAAACAAGTTCGGCATCGCCATTGCCGCTGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCg                          >  1:19299/1‑149 (MQ=255)                                                                                                                            aGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTCGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGc                         <  1:243970/149‑1 (MQ=255)                                                                                                                            aGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGcc                        <  6:128874/149‑1 (MQ=255)                                                                                                                                aaCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCctgct                    >  3:256870/1‑149 (MQ=255)                                                                                                                                   aaGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggt                 <  5:270281/149‑1 (MQ=255)                                                                                                                                     gTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGAAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgc               >  1:15650/1‑149 (MQ=255)                                                                                                                                     cctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgc               >  2:265513/3‑149 (MQ=255)                                                                                                                                      ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTc                                        >  2:39221/2‑123 (MQ=255)                                                                                                                                      ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTc                                        <  1:39221/122‑1 (MQ=255)                                                                                                                                       tcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctg             >  8:334842/1‑149 (MQ=255)                                                                                                                                            aTCGGCAATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGa        <  1:265513/149‑1 (MQ=255)                                                                                                                                               gCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCg     >  6:174590/1‑149 (MQ=255)                                                                                                                                                ccATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGc    >  5:190892/1‑149 (MQ=255)                                                                                                                                                  aTTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGCCt  <  6:342464/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CACCGACGAGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGCCT  >  NC_002947/5961482‑5961775

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