Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,824,151	A→G	5.7%	F2347L (TTC→CTC)	pvdL ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,824,151	0	A	G	5.7%	210.6 / 2.1	71	F2347L (TTC→CTC)	pvdL	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (28/38);  new base G (2/2);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTC  >  NC_002947/4824003‑4824296                                                                                                                                                     |                                                                                                                                                   tGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGaa                                                                                                                                                    >  4:77228/1‑149 (MQ=255) tGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGaa                                                                                                                                                    >  4:114129/1‑149 (MQ=255)   ggCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAAc                                                                                                                                                   <  5:379472/148‑1 (MQ=255)   ggCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACt                                                                                                                                                  <  2:121718/149‑1 (MQ=255)    gcgcgcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGACGAATGCGGCGATCACCGCCTGGCACTGCATGGCGAACTc                                                                                                                                                 >  4:21539/1‑149 (MQ=255)       cgcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGc                                                                                                                                               >  3:139997/1‑148 (MQ=255)           gTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATgg                                                                                                                                          <  8:168125/149‑1 (MQ=255)           gTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATgg                                                                                                                                          <  5:18251/149‑1 (MQ=255)           gTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATgg                                                                                                                                          <  4:264464/149‑1 (MQ=255)             aCCACGCGCGCACCACACAACAGCGTCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAGTTCGCCATgca                                                                                                                                         <  3:313737/148‑3 (MQ=255)               cACGCGCGCACCACACAACAGCGGAACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAGATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGa                                                                                                                                      >  8:251956/1‑149 (MQ=255)               cACGCCCGCACCACACAACAGCGCCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCTCAATCATCGATGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGa                                                                                                                                      <  4:359843/149‑1 (MQ=255)                 cgcgcgCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGacc                                                                                                                                    <  8:37767/149‑1 (MQ=255)                      gcACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACcaccac                                                                                                                               <  3:93671/149‑1 (MQ=255)                        accacACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGc                                                                                                                             <  5:94419/149‑1 (MQ=255)                                acaGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGt                                                                                                                     <  2:306601/149‑1 (MQ=255)                                  aGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGttt                                                                                                                   <  7:183407/149‑1 (MQ=255)                                  aGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGttt                                                                                                                   <  3:272719/149‑1 (MQ=255)                                      gCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCg                                                                                                               <  6:347335/149‑1 (MQ=255)                                         cccACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTg                                                                                                            <  2:280685/149‑1 (MQ=255)                                          ccACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAGTTCGCCATgcact                                                                                                                                       <  2:288338/121‑5 (MQ=255)                                          ccACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAGTTCGCCATgcact                                                                                                                                       >  1:288338/1‑117 (MQ=255)                                          ccACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTg                                                                                                            >  3:267872/1‑148 (MQ=255)                                          ccACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTg                                                                                                            >  4:283346/1‑148 (MQ=255)                                           cacaTGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGa                                                                                                          <  2:308175/149‑1 (MQ=255)                                           cacaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTg                                                                                                            >  7:220963/1‑147 (MQ=255)                                             caGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAGATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGc                                                                                                        <  7:251956/149‑1 (MQ=255)                                                gcgcTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCt                                                                                                      >  1:207488/1‑148 (MQ=255)                                                  ggtcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAg                                                                                                   <  1:296975/147‑1 (MQ=255)                                                  gctcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAg                                                                                                   <  3:283346/149‑1 (MQ=255)                                                  gctcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAg                                                                                                   <  2:207488/149‑1 (MQ=255)                                                    tcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTgg                                                                                                           >  6:358810/1‑139 (MQ=255)                                                    tcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTgg                                                                                                           <  5:358810/139‑1 (MQ=255)                                                      gctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGgt                                                                                                 <  8:220963/147‑1 (MQ=255)                                                      gctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGACCCTGAGgtgt                                                                                               <  7:245344/149‑1 (MQ=255)                                                           cAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCTATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCGCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATc                                                                                          <  3:21539/149‑1 (MQ=255)                                                             gCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAg                                                                                        <  4:139997/149‑1 (MQ=255)                                                               gTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAg‑g                                                                                      >  3:309155/1‑148 (MQ=255)                                                                   aaaTTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAC‑GTAAga                                                                                 >  2:233520/1‑148 (MQ=255)                                                                        gATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGt                                                                            >  5:98960/1‑149 (MQ=255)                                                                         aTCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGtg                                                                           >  5:245763/1‑149 (MQ=255)                                                                         aTCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGt                                                                            <  3:114129/148‑1 (MQ=255)                                                                            gAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGtgctg                                                                        <  7:149457/149‑1 (MQ=255)                                                                             aaTAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGtgctgc                                                                       >  5:76848/1‑149 (MQ=255)                                                                                 gAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGc                                                                    <  6:98960/148‑1 (MQ=255)                                                                                 gAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCa                                                                   <  7:130654/149‑1 (MQ=255)                                                                                    aTGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGc                                                                >  7:328928/1‑149 (MQ=255)                                                                                         gCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGCCGCAGCACCAGGCttacc                                                           <  3:114549/149‑5 (MQ=37)                                                                                             gCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGtt                                                         <  3:77228/147‑1 (MQ=255)                                                                                               aaTGATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGTGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTc                                                      <  1:233520/149‑1 (MQ=255)                                                                                                    tcGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGACCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCgg                                                >  6:81741/1‑149 (MQ=255)                                                                                                           gcATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGACCTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCg                                         >  3:4019/1‑149 (MQ=255)                                                                                                           gcATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGt                                                                            <  1:47721/114‑1 (MQ=255)                                                                                                           gcATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGt                                                                            >  2:47721/1‑114 (MQ=255)                                                                                                           gcATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAC‑GTAAGCCTGCTGATCATTAAATATGAAGTTGTCCAGCGGTGCATCg                                         >  2:27754/1‑149 (MQ=21)                                                                                                             aTGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGATGCACCGCt                                       >  6:120032/1‑149 (MQ=255)                                                                                                              tGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTg                                      >  1:331012/1‑149 (MQ=255)                                                                                                                 cGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGATGCACCGCTggac                                   <  5:120032/149‑4 (MQ=255)                                                                                                                 cGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGACTGGTTCTGATTCATTATGAGGTTGTCCAGCGGTGCATCGCTGTAg                                   >  2:184223/1‑149 (MQ=255)                                                                                                                        ggcgATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGa                                                             >  6:209972/1‑116 (MQ=255)                                                                                                                        ggcgATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGa                                                             <  5:209972/116‑1 (MQ=255)                                                                                                                              cACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGACTCGTGCTGCGTCATGATGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTgg                      >  8:278548/1‑149 (MQ=255)                                                                                                                               aCCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCa                            <  8:327387/142‑1 (MQ=255)                                                                                                                               aCCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCa                            >  7:327387/1‑142 (MQ=255)                                                                                                                                 cGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggcgg                   >  3:310043/1‑149 (MQ=255)                                                                                                                                       ggtgTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTTGGTTTGCCGGTGGAGCCTGAGGGGTAGATCAG‑GTAAGCCTGGTGCTGCGGCTGGCTGAGGTTGTCCAGCGGTGCTTCGCTGTAGGCCGCCAGGCTGGCGGCCTCGt             >  5:260475/4‑149 (MQ=38)                                                                                                                                       gCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCg              <  6:245763/148‑1 (MQ=255)                                                                                                                                       gCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGt             >  5:260473/1‑149 (MQ=255)                                                                                                                                          agccATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAgg                                                                 >  1:327947/4‑94 (MQ=255)                                                                                                                                          agccATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAgg                                                                 <  2:327947/91‑1 (MQ=255)                                                                                                                                           tGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTc         >  7:52747/1‑149 (MQ=255)                                                                                                                                            ccATGGCGAGTTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTc         >  3:237843/2‑148 (MQ=255)                                                                                                                                                  cGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTc  <  4:257485/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAG‑GTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTC  >  NC_002947/4824003‑4824296

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