Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,687,701	Δ1 bp	6.4%	coding (677/1074 nt)	PP_4148 →	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,687,701	0	T	.	6.4%	197.8 / 6.5	62	coding (677/1074 nt)	PP_4148	ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base . (2/2);  total (30/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 = 9.65e-01

CATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCATGC  >  NC_002947/4687555‑4687847                                                                                                                                                    |                                                                                                                                                    cATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTgg                                                                                                                                                   >  5:144668/1‑149 (MQ=255)  aTCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGt                                                                                                                                                  <  7:355541/149‑1 (MQ=255)               tGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTc                                                                                                                                     <  5:168317/149‑1 (MQ=255)               tGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTc                                                                                                                                     >  3:256800/1‑149 (MQ=255)                gCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCg                                                                                                                                    >  8:7321/1‑149 (MQ=255)                 cTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGc                                                                                                                                   <  1:299883/149‑1 (MQ=255)                   gTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGc                                                                                                                                           >  1:276147/1‑139 (MQ=39)                   gTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGc                                                                                                                                           <  2:276147/139‑1 (MQ=39)                      cGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGccacc                                                                                                                               <  5:18758/148‑1 (MQ=255)                        ccGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAGCAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGCa                                                                                                                                          >  1:292963/1‑135 (MQ=255)                        ccGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAGCAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGCa                                                                                                                                          <  2:292963/135‑1 (MQ=255)                                     cTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGAcc                                                                                                                <  1:6223/148‑1 (MQ=255)                                           cTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACa                                                                                                         <  5:54927/149‑1 (MQ=255)                                            tGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTc                                                                                                                                                 <  8:308861/108‑1 (MQ=255)                                            tGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTc                                                                                                                                                 >  7:308861/1‑108 (MQ=255)                                                 cccGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTc                                                                                                    <  5:169877/148‑1 (MQ=255)                                                 cccGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTc                                                                                                    <  6:371925/148‑1 (MQ=255)                                                   cGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTc                                                                                                 >  2:156141/1‑149 (MQ=255)                                                    gCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGc                                                                                                                                   <  2:48539/114‑1 (MQ=255)                                                    gCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGc                                                                                                                                   >  1:48539/1‑114 (MQ=255)                                                    gCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTc                                                                                                 >  1:100744/1‑148 (MQ=255)                                                              atCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGACGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGAACCTGCTGACAAAGAACAGCTTCCTCGACGAAATcaa                                                                                      >  1:242205/2‑149 (MQ=255)                                                               tCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATcaa                                                                                      >  3:177477/1‑148 (MQ=255)                                                                  ccTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCaacaaa                                                                                  <  2:100744/149‑1 (MQ=255)                                                                   cTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCaacaaa                                                                                  >  7:108345/1‑148 (MQ=255)                                                                     ccAACAACTTCGACGAACTCAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAg                                                                               <  4:307969/149‑1 (MQ=255)                                                                      caacaaCTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTAATCGGCAACTTATCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGt                                                                              <  5:228996/149‑1 (MQ=255)                                                                       aacaacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTa                                                                             <  8:108345/149‑1 (MQ=255)                                                                       aacaacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTa                                                                             >  6:37469/1‑149 (MQ=255)                                                                         caacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACg                                                                           >  3:345372/1‑149 (MQ=255)                                                                          aacTTCGACGCACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTAGACACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAAGCAACAACCAGTACGt                                                                          >  6:279898/1‑149 (MQ=39)                                                                          aacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCAGCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGt                                                                          <  5:228487/149‑1 (MQ=255)                                                                          aacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGt                                                                          >  8:301187/1‑149 (MQ=255)                                                                           acTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTg                                                                         <  3:155424/149‑1 (MQ=255)                                                                             ttCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGt                                                                       >  6:193/1‑149 (MQ=255)                                                                             ttCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGt                                                                       >  1:146550/1‑149 (MQ=255)                                                                                gacgaACTTAGCACCACCTGCATGGTGCTGGACTACTTCTGGCACGTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAATAACAGCTTCGTCGACGAAATCAACAAACAGTACGTGGTCAc                                                                    <  5:279898/149‑1 (MQ=255)                                                                                    aaCTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTTACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTTGTCACCgcc                                                                <  4:177477/149‑1 (MQ=255)                                                                                    aaCTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgc                                                                 >  4:371648/1‑148 (MQ=255)                                                                                         aGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggc                                                           >  3:367183/1‑149 (MQ=255)                                                                                            accaccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggcc                                                          >  8:205464/1‑147 (MQ=255)                                                                                              caccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGccaaggccaagg                                                      >  5:203673/1‑149 (MQ=255)                                                                                              caccGGCAAGGTGCTGGACTACTTCTGGCACATGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGATATCAACAAACAGTACGTGGTCACCCccaaggccaagg                                                      >  6:283744/1‑149 (MQ=255)                                                                                                  ggCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCt                                                  <  1:66033/149‑1 (MQ=255)                                                                                                   gCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTTGTCACCGCCAAGGCCAAGGGCCTg                                                 <  6:203673/149‑1 (MQ=255)                                                                                                      aGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGc                                              <  4:367183/149‑1 (MQ=255)                                                                                                       ggTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGcc                                             <  7:163670/149‑1 (MQ=255)                                                                                                        gTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGcc                                             >  8:86352/1‑148 (MQ=255)                                                                                                            tGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCa                                        >  3:291719/1‑149 (MQ=255)                                                                                                              gACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACg                                      >  7:286037/1‑149 (MQ=255)                                                                                                               actactTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACgg                                     <  5:193/149‑1 (MQ=255)                                                                                                               actactTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACgg                                     <  2:23528/149‑1 (MQ=255)                                                                                                                    ttCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACAAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCt                                >  7:317586/1‑149 (MQ=255)                                                                                                                        ggCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTAc                           >  5:181231/1‑149 (MQ=255)                                                                                                                         gCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTAc                           <  3:371648/148‑1 (MQ=255)                                                                                                                                gTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCAc                     <  6:181231/147‑1 (MQ=255)                                                                                                                                gTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCAACATGACCCTGCTGACCTAGAACAGCTTCGTCGACGAAATCAACATTCAGTACGGTGTCACCGCCAAGGCCAAGGTCCTTAGCCTATCACGGGTTCT‑GTACTGCCACgt                   <  2:242205/149‑1 (MQ=255)                                                                                                                                       cGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACgt                   <  4:126207/142‑1 (MQ=255)                                                                                                                                       cGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACgt                   >  3:126207/1‑142 (MQ=255)                                                                                                                                             cc‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGATGCT‑GTACGGCCACGTGTTCCGCAACGcc      <  1:156141/149‑1 (MQ=255)                                                                                                                                              c‑gcgCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCa     >  3:232101/1‑149 (MQ=255)                                                                                                                                                  gcTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCAtgc  <  4:256800/149‑1 (MQ=255)                                                                                                                                                    tGGTCATCGGAAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCa     >  2:3213/1‑144 (MQ=255)                                                                                                                                                    tGGTCATCGGAAACATCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAACGCCAAGGGCCTGAGCCGATCACGGGTGCTCGTACGGCCACGTGTTCCGCAACGCCa     <  1:3213/145‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCATGC  >  NC_002947/4687555‑4687847

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