Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,754,404	T→G	7.2%	V260G (GTG→GGG)	pvdR →	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,754,404	0	T	G	7.2%	176.3 / 3.5	64	V260G (GTG→GGG)	pvdR	ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (26/33);  new base G (2/3);  total (28/36)
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

CGCCGCTGATCCTGCGCATCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGACAACCCCGA  >  NC_002947/4754260‑4754542                                                                                                                                                  |                                                                                                                                           cgccgcTGATCCTGCGCATCGCCAAGCTGTACCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGcc                                                                                                                                       <  2:271307/150‑1 (MQ=255)  gccgcTGATCCTGCGCATCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCa                                                                                                                                      >  1:578968/1‑150 (MQ=255)  gccgcTGATCCTGCGCATCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCa                                                                                                                                      >  1:661314/1‑150 (MQ=255)         aTCCTGCGCATCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCTGGCGCACGGtgctg                                                                                                                               <  2:942015/150‑1 (MQ=255)                   tCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGc                                                                                                                     <  1:17532/150‑1 (MQ=255)                   gtgCCAAGCTGTCCCCCGTGACCGGCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGc                                                                                                                     <  2:453718/148‑1 (MQ=255)                       cAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCTGGCGCACGGtgctg                                                                                                                               <  2:232524/136‑1 (MQ=255)                       cAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCTGGCGCACGGtgctg                                                                                                                               >  1:232524/1‑136 (MQ=255)                        aaGCTGTCCCCTATGGCGGTCTGGGCCGAGGTCTCGGACGCCACTACCGAAAACATCAAGACCAGCATGAACGCCTACTTCACCACTCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAg                                                                                                                <  2:518922/150‑1 (MQ=255)                         aGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGc                                                                                                               <  1:1263788/150‑1 (MQ=255)                              tCCCCCATGGCCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGtgctg                                                                                                                               >  2:676667/1‑129 (MQ=255)                              tCCCCCATGGCCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGtgctg                                                                                                                               <  1:676667/129‑1 (MQ=255)                              tCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTg                                                                                                          >  2:868918/1‑150 (MQ=255)                               cccccATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTgg                                                                                                         >  2:808561/1‑150 (MQ=255)                               cccccATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTgg                                                                                                         >  2:231602/1‑150 (MQ=255)                                 cccATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGTGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgac                                                                                                       >  2:319227/1‑150 (MQ=255)                                  ccATGACCGTCTGGGCCCAGGTCTCGGAGGCGGATAGCTGCAAGATCCATCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgacc                                                                                                      >  2:1019875/1‑150 (MQ=255)                                  ccATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgacc                                                                                                      >  2:29564/1‑150 (MQ=255)                                   cATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgacca                                                                                                     >  1:1409817/1‑150 (MQ=255)                                     tGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgaccaga                                                                                                   <  2:152450/150‑1 (MQ=255)                                        ccGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgaccagacca                                                                                                <  1:667014/150‑1 (MQ=255)                                            cTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTCCTTCCCCACCCTGGCCTTCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagcca                                                                                            <  1:882394/150‑1 (MQ=255)                                             tGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagccag                                                                                           <  1:231602/150‑1 (MQ=255)                                                   cAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGggcggc                                                                                     <  1:1061690/150‑1 (MQ=255)                                                          cGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGcagcc                                                                              <  1:808561/150‑1 (MQ=255)                                                          cGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTACCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGcagcc                                                                              <  2:578968/150‑1 (MQ=255)                                                               gCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGcagcccagcc                                                                         >  2:353626/1‑150 (MQ=255)                                                                  gATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCcagccagc                                                                      >  2:1062911/1‑150 (MQ=255)                                                                  gATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCcagccagc                                                                      <  1:29564/150‑1 (MQ=255)                                                                      tCGGCAAGATCCAGCCTGGTATGACCGCCTACTTCACCCCCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCa                                                                  <  2:939267/150‑1 (MQ=255)                                                                        ggCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCTGGCGCACGGtgctg                                                                                                                               >  1:805581/1‑87 (MQ=255)                                                                        ggCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCTGGCGCACGGtgctg                                                                                                                               <  2:805581/87‑1 (MQ=255)                                                                            aGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCAccgccc                                                            >  2:1417948/1‑149 (MQ=255)                                                                                    ccTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCAcc                                                    >  2:319809/1‑150 (MQ=255)                                                                                        gCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagccag                                                                                           >  1:8867/1‑107 (MQ=255)                                                                                        gCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagccag                                                                                           <  2:8867/107‑1 (MQ=255)                                                                                          atccccGCATACTTCAACACCCCGGCCGGCGGC‑AACCCCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGc                                              <  2:1084069/146‑1 (MQ=255)                                                                                             aCCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGCGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCaa                                           <  2:556936/150‑1 (MQ=255)                                                                                              ccGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCAGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAg                                          >  1:394392/1‑150 (MQ=255)                                                                                               cGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGt                                         >  1:780314/1‑150 (MQ=255)                                                                                                  cTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCAGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGt                                      <  2:394392/150‑1 (MQ=255)                                                                                                   tACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGCGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTg                                     <  1:543389/150‑1 (MQ=255)                                                                                                    aCTTCACCACCCTGGCCGGCGGCAAAACGCC‑CTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGc                                    >  1:1327847/1‑150 (MQ=255)                                                                                                     cTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCa                                   <  1:1417948/150‑1 (MQ=255)                                                                                                      ttCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCaa                                                                                                                 <  1:93161/71‑1 (MQ=255)                                                                                                      ttCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCaa                                                                                                                 >  2:93161/1‑71 (MQ=255)                                                                                                      ttCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGTCCAACAAGGGACTCCTTCGTCTCCTCCCGCAGCCTAGTCGTGCac                                  >  2:156027/1‑149 (MQ=255)                                                                                                      ttCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCaa                                  >  2:345199/1‑150 (MQ=255)                                                                                                           caccCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATacac                             <  2:780314/150‑1 (MQ=255)                                                                                                              ccTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagcca                                                                                            >  1:1060661/1‑84 (MQ=255)                                                                                                              ccTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagcca                                                                                            <  2:1060661/84‑1 (MQ=255)                                                                                                              ccTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagc                                                                                              <  1:1005518/82‑1 (MQ=255)                                                                                                              ccTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagc                                                                                              >  2:1005518/1‑82 (MQ=255)                                                                                                                tGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACGGTGCGCACGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCACCGACACCGCCGGACCCACCGGTAGCCAGGTCGCTCTATACTCCGtgc                        <  2:1156099/150‑1 (MQ=255)                                                                                                                    cggcggC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTgctgct                    <  1:868918/150‑1 (MQ=255)                                                                                                                     ggcggc‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACCCCGAGctgctg                   <  2:354577/150‑1 (MQ=255)                                                                                                                        ggc‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGAc                <  2:1409817/150‑1 (MQ=255)                                                                                                                          c‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGt              <  1:353626/150‑1 (MQ=255)                                                                                                                          c‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGt              <  1:266608/150‑1 (MQ=255)                                                                                                                              gcgcCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGa           <  1:1062911/150‑1 (MQ=255)                                                                                                                              gcgcCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACAAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATCCACCGTGCTGCTGGACGTGGa           >  1:350569/1‑150 (MQ=255)                                                                                                                                       gACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGTCGGCAGCGGCCCAGCCAGCGCAACCGCCGACGCCACCGGCAGCGAAGTCCTCCAATACGCCGGGCTGCTCGACGTGGACGAcccgaa  >  2:140784/1‑147 (MQ=255)                                                                                                                                         ccAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGgcgc                                                        <  2:649410/94‑1 (MQ=255)                                                                                                                                         ccAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGgcgc                                                        >  1:649410/1‑94 (MQ=255)                                                                                                                                                  |                                                                                                                                           CGCCGCTGATCCTGCGCATCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGC‑AAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGACAACCCCGA  >  NC_002947/4754260‑4754542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑