Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,687,693	C→G	7.2%	I223M (ATC→ATG)	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,693	0	C	G	7.2%	212.4 / 5.1	68	I223M (ATC→ATG)	PP_4148	ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (36/27);  new base G (3/2);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACGTGTTCCG  >  NC_002947/4687548‑4687836                                                                                                                                                  |                                                                                                                                                  tGTTCGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGAc                                                                                                                                                 >  8:41060/1‑148 (MQ=255) tGTTCGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCAcc                                                                                                                                                >  1:252317/1‑149 (MQ=255)     cGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑gcgc                                                                                                                                           >  7:54788/1‑149 (MQ=255)       ccATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTCCCGATCCCC‑GCGCTg                                                                                                                                         >  5:279419/1‑149 (MQ=255)                tgATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGc                                                                                                                                <  8:54788/149‑1 (MQ=255)                  aTCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGCag                                                                                                                              <  7:41060/149‑2 (MQ=39)                  aTCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGCa                                                                                                                               <  3:185818/148‑1 (MQ=39)                   tCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCAcc                                                                                                                                                >  3:270368/1‑131 (MQ=255)                   tCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCAcc                                                                                                                                                <  4:270368/131‑1 (MQ=255)                           ttCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGccac                                                                                                                    >  7:18537/1‑149 (MQ=255)                             cGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGccacca                                                                                                                  >  1:144152/1‑149 (MQ=255)                                     gCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATTGGCAAC‑TTCGCCACCATGACCctg                                                                                                          >  2:153354/1‑149 (MQ=255)                                           aCTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGc                                                                                                                       <  7:250750/130‑1 (MQ=255)                                           aCTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGc                                                                                                                       >  8:250750/1‑130 (MQ=255)                                           aCTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGAc                                                                                                     >  1:367565/1‑148 (MQ=255)                                               cGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAg                                                                                                 >  8:162155/1‑148 (MQ=255)                                               cGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAg                                                                                                 <  6:35768/148‑1 (MQ=255)                                               cGACTGGTTCCCGCTGCGTGCCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCaa                                                                                                  <  4:113251/147‑1 (MQ=255)                                                     gTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGc                                                                                           <  8:18537/148‑1 (MQ=255)                                                      ttCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAc                                                                                                                             <  2:59199/114‑1 (MQ=255)                                                      ttCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAc                                                                                                                             >  1:59199/1‑114 (MQ=255)                                                      ttCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCtt                                                                                         >  4:295843/1‑149 (MQ=255)                                                             tGCGTGGCCTCACCTACAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTcgac                                                                                  >  4:169996/1‑149 (MQ=255)                                                               cGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTcgacga                                                                                <  1:159882/149‑1 (MQ=255)                                                               cGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTcgac                                                                                  <  7:162155/147‑1 (MQ=255)                                                                 tGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATTGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGaaa                                                                              <  1:153354/149‑1 (MQ=255)                                                                   gCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATc                                                                            <  2:252317/149‑1 (MQ=255)                                                                     cTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTc                                                                                        >  2:8097/1‑135 (MQ=255)                                                                     cTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTc                                                                                        <  1:8097/135‑1 (MQ=255)                                                                        aCCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATcaaca                                                                        >  7:75053/1‑148 (MQ=255)                                                                        aCCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATcaac                                                                         <  3:295843/147‑1 (MQ=255)                                                                          cTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCaa                                                                                                                              >  7:19500/1‑93 (MQ=255)                                                                          cTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCaa                                                                                                                              <  8:19500/93‑1 (MQ=255)                                                                              aacaacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCATC‑GCACTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCaacaaa                                                                      <  1:371063/144‑1 (MQ=255)                                                                              aacaacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCAAC‑GCACTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCaacaaa                                                                      >  2:371063/1‑144 (MQ=255)                                                                                  acTTCGACGAACTGAGCACCACCGGCTAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTg                                                             <  1:346255/149‑1 (MQ=255)                                                                                  acTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTg                                                             <  3:169996/149‑1 (MQ=255)                                                                                    ttCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGAAC‑GCGCTGGTCATCGGCAAC‑TTCGCCCCCATGACCCTGCTTACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTgg                                                            >  4:195030/1‑148 (MQ=255)                                                                                    ttCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGt                                                           >  5:296388/1‑149 (MQ=255)                                                                                       gacgaACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCa                                                         <  8:75053/148‑1 (MQ=255)                                                                                                aGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaag                                                 <  2:367565/147‑1 (MQ=255)                                                                                                   accaccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACGCTGCTGACCAAGAAc                                                                                              <  4:248998/99‑1 (MQ=255)                                                                                                   accaccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCtttc                                                                                       <  3:126397/106‑3 (MQ=255)                                                                                                   accaccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCtttc                                                                                       >  4:126397/1‑104 (MQ=255)                                                                                                   accaccGGCAAGGTCCTGGACAACTTCTGGCACCTGGTCCTGCCCATCACA‑CCGCTGGTCATCGGCAAC‑TTCGCCACCATGACACTGCTGACCAATAAc                                                                                              >  3:248998/1‑99 (MQ=255)                                                                                                    ccaccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCaa                                                                                                  >  8:263947/1‑94 (MQ=255)                                                                                                    ccaccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggccaag                                           <  5:341086/149‑1 (MQ=255)                                                                                                    ccaccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggccaag                                           >  6:341086/1‑149 (MQ=255)                                                                                                     caccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCctgct                                                                                                        <  6:48826/87‑1 (MQ=255)                                                                                                             aGGTGCTGGACTACTTCTGGCACGTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAAGCATCAAACAGTACGTTGTCACCGCCAAGGCCAAGGGGCTGAGc                                  <  6:279419/149‑1 (MQ=255)                                                                                                                    ggACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCAc                           >  2:173614/1‑149 (MQ=255)                                                                                                                     gACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACg                          >  4:88339/1‑149 (MQ=255)                                                                                                                     gACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACg                          >  2:97209/1‑149 (MQ=255)                                                                                                                        tactTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTAc                                                                <  4:109630/108‑1 (MQ=255)                                                                                                                        tactTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTAc                                                                >  3:109630/1‑108 (MQ=255)                                                                                                                          ctTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACGTTCGCCACCATGACCCTGCTGACAAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTg                      >  1:268831/1‑149 (MQ=255)                                                                                                                           ttCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCt                    >  8:207270/1‑149 (MQ=255)                                                                                                                           ttCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCt                    >  8:6282/1‑149 (MQ=255)                                                                                                                                gCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACg               >  5:274528/1‑149 (MQ=255)                                                                                                                                 cACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACgg              >  8:277683/1‑149 (MQ=255)                                                                                                                                      ggTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACg         >  6:326614/1‑149 (MQ=255)                                                                                                                                      gaagttgcCGATGACCAGCGC‑GGTGATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACg         >  5:303073/6‑149 (MQ=255)                                                                                                                                       gTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAGATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACgt        >  5:219479/1‑149 (MQ=255)                                                                                                                                       gTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCAc          <  1:173614/147‑1 (MQ=255)                                                                                                                                           tgcCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACGTGTTc    >  8:69406/1‑149 (MQ=255)                                                                                                                                            gcCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACGTGTTcc   >  8:69508/1‑149 (MQ=255)                                                                                                                                             ccGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACGTGTTCCg  <  2:144152/149‑1 (MQ=255)                                                                                                                                             ccGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACGTGTTCCg  >  4:186390/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TGTTCGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAAC‑TTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACGTGTTCCG  >  NC_002947/4687548‑4687836

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