Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,857,061	C→T	6.1%	R36W (CGG→TGG)	PP_5671 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,857,061	0	C	T	6.1%	192.7 / 2.1	65	R36W (CGG→TGG)	PP_5671	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base T (2/2);  total (33/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.58e-01

GCCCGCCCCCCAGCACCTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCACCGCGCCGCCCAGGATCACA  >  NC_002947/4856913‑4857207                                                                                                                                                     |                                                                                                                                                   gcccgcccCCCAGCACCTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTcc                                                                                                                                                    <  8:125502/149‑1 (MQ=255)  cccgcccCCCAGCACCTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCg                                                                                                                                                   <  8:548/149‑1 (MQ=255)      ccccccAGCACCTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCag                                                                                                                                               >  8:105839/1‑149 (MQ=255)        ccccAGCACCTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCagga                                                                                                                                             <  3:313597/149‑1 (MQ=255)         cccAGCACCTTCGCCGCAACAGAATAGACAGCCCGATTCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCaggag                                                                                                                                            <  2:92248/149‑1 (MQ=255)         cccAGCACCTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAAGTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCaggag                                                                                                                                            <  4:24665/149‑1 (MQ=255)         cccAGCACCTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCaggag                                                                                                                                            <  7:248664/149‑1 (MQ=255)           cAGCACCTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCTGGCAGGAGGc                                                                                                                                          >  1:221810/1‑149 (MQ=255)                 cTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCagga                                                                                                                                             <  3:62447/140‑1 (MQ=255)                 cTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCagga                                                                                                                                             >  4:62447/1‑140 (MQ=255)                 cTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTg                                                                                                                                    >  5:137275/1‑149 (MQ=255)                        gcAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCa                                                                                                                             <  6:131280/149‑1 (MQ=255)                               aaTAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCaa                                                                                                                      <  7:100829/149‑1 (MQ=255)                                    aCAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGt                                                                                                                 >  6:347444/1‑149 (MQ=255)                                      aGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGtt                                                                                                                >  1:51905/1‑148 (MQ=255)                                        cccGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTcag                                                                                                             <  1:342615/149‑1 (MQ=255)                                           gATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTcagcag                                                                                                          <  7:105839/149‑1 (MQ=255)                                                  cAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCAc                                                                                                   >  8:48210/1‑149 (MQ=255)                                                   aGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTg                                                                                                                                    >  2:56380/1‑115 (MQ=255)                                                   aGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTg                                                                                                                                    <  1:56380/115‑1 (MQ=255)                                                    ggATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGa                                                                                                 <  2:51905/149‑1 (MQ=255)                                                     gatgatCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGAc                                                                                                >  5:43738/1‑149 (MQ=255)                                                         atCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCTGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCAt                                                                                             <  2:221810/148‑1 (MQ=255)                                                          tCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCa                                                                                                                             <  8:179640/115‑1 (MQ=255)                                                          tCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCa                                                                                                                             >  7:179640/1‑115 (MQ=255)                                                                    ctcGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCTGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAAt                                                                                 >  4:42918/1‑149 (MQ=255)                                                                       gAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAAt                                                                                 <  3:185912/146‑1 (MQ=255)                                                                           ttGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTcagc                                                                                                            >  8:285598/1‑115 (MQ=255)                                                                           ttGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTg                                                                          >  8:14451/1‑149 (MQ=255)                                                                             gagaAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGc                                                                        >  4:327678/1‑149 (MQ=255)                                                                              agaAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGc                                                                        >  3:286320/1‑148 (MQ=255)                                                                               gaAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGc                                                                        >  3:250201/1‑147 (MQ=255)                                                                                     aCAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAg                                                                 <  6:43738/148‑1 (MQ=255)                                                                                       aaaaTGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGctc                                                              >  7:252014/1‑149 (MQ=255)                                                                                       aaaaTGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGctc                                                              >  5:41164/1‑149 (MQ=255)                                                                                       aaaaTGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGctc                                                              >  5:262969/1‑149 (MQ=255)                                                                                         aaTGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTca                                                                                                              >  5:211015/1‑99 (MQ=255)                                                                                         aaTGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTca                                                                                                              <  6:211015/99‑1 (MQ=255)                                                                                         aaTGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGctct                                                             >  4:108072/1‑148 (MQ=255)                                                                                          aTGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGctct                                                             >  1:172809/1‑147 (MQ=255)                                                                                               ggCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTc                                                      <  6:41164/149‑1 (MQ=255)                                                                                                 ccGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCTGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGa                                                    <  3:42918/149‑1 (MQ=255)                                                                                                 ccGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGa                                                    >  6:39752/1‑149 (MQ=255)                                                                                                 ccGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGa                                                    >  6:45247/1‑149 (MQ=255)                                                                                                      aGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAAt                                               <  4:286320/149‑1 (MQ=255)                                                                                                      aGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAAt                                               <  6:137275/149‑1 (MQ=255)                                                                                                       gTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATc                                              <  3:108072/149‑1 (MQ=255)                                                                                                          aaGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCa                                            <  2:172809/148‑1 (MQ=255)                                                                                                            gggCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGc                                          >  6:150099/1‑148 (MQ=255)                                                                                                                  aGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAgg                                   <  6:262969/149‑1 (MQ=255)                                                                                                                   ggCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGa                                  <  4:250201/149‑1 (MQ=255)                                                                                                                     cacaGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGt                                <  7:48210/149‑1 (MQ=255)                                                                                                                      acaGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGt                                >  8:59393/1‑148 (MQ=255)                                                                                                                      acaGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACCACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGt                                >  4:348350/1‑148 (MQ=255)                                                                                                                           aaTTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACaaa                          <  7:14451/149‑1 (MQ=255)                                                                                                                                gTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCa                      <  5:150099/148‑1 (MQ=255)                                                                                                                                gTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCa                      <  7:59393/148‑1 (MQ=255)                                                                                                                                  aGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGc                                          >  8:276217/1‑126 (MQ=255)                                                                                                                                   gCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCACcgc                  <  5:347444/149‑1 (MQ=255)                                                                                                                                   gCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCACcg                   <  5:39752/148‑1 (MQ=255)                                                                                                                                   gCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCACcg                   <  5:45247/148‑1 (MQ=255)                                                                                                                                    cGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCTACGGGTTCAGCAGTTGTCACGACCATAGATGGCGTAATGGCACTGGCCAGCAAGTTGTGAGCTTCGATTAAGCCAGCCACAAGGAGTAACAAATGCACagcg                 >  6:208709/1‑149 (MQ=255)                                                                                                                                    cGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCACcgcg                 >  5:203015/1‑149 (MQ=255)                                                                                                                                           gcgTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCACCGCgccgcc            >  5:237859/1‑147 (MQ=255)                                                                                                                                                   tCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCACCGCGCCGCCCAGGATcaca  >  7:262742/1‑149 (MQ=255)                                                                                                                                                   tCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCACCGCGCCGCCCAGGATcaca  >  6:15548/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCCCGCCCCCCAGCACCTTCGCCGCAACAGAATAGACAGCCCGATGCATCAGGATGATCGTGAACATCTCGAACTTGAGAAGTTACAAAATGAAGGCCGCAAGTTAAGGGCCGAGGCACAGAAATTTGTAGCGGAAGCGCGTAAGCTCCGGCAGGAGGCGTTCTGGTACCCAATCGCAACGGGTTCAGCAGTTGTCACGACCATAGTTGGCGTAATGGCACTGGCCAGCAAGCTCTGAGCTTCGATTAATCCAGCCACAAGGAGTAACAAATGCACCGCGCCGCCCAGGATCACA  >  NC_002947/4856913‑4857207

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