Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,616,314	A→G	5.8%	Q74R (CAA→CGA)	bepA →	OM protein maintenance and assembly metalloprotease and chaperone, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,616,314	0	A	G	5.8%	211.0 / 3.1	69	Q74R (CAA→CGA)	bepA	OM protein maintenance and assembly metalloprotease and chaperone, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (30/35);  new base G (2/2);  total (32/37)
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.37e-01

GTTTGCCGACAGCGCAGACACCTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACGTGGTGCTGCACTCTGCCCTGTTCCGTTA  >  NC_000913/2616168‑2616461                                                                                                                                                   |                                                                                                                                                    gTTTGCCGACAGCGCAGACACCTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAat                                                                                                                                                   <  5:87258/149‑1 (MQ=255)          cAGCGCAGACACCTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAAtt                                                                                                                                          <  4:275277/149‑1 (MQ=255)           agggcAGACACGTTGCCGGATATGGGTACCTCCGCAGGAAGCACGCTTTCGATTGGTCATGAAATGCAGATGGGCGACTATTATGTCCGCCATCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTc                                                                                                                                         <  5:33891/146‑1 (MQ=255)              gcAGACACCTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCt                                                                                                                                      >  4:190650/1‑149 (MQ=255)               cAGACACCTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTg                                                                                                                                     <  5:333959/149‑1 (MQ=255)                aGACACCTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCTCTgg                                                                                                                                    >  1:106316/1‑149 (MQ=255)                aGACACCTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTgg                                                                                                                                    <  5:197760/149‑1 (MQ=255)                      cTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGc                                                                                                                              >  3:279295/1‑149 (MQ=255)                      cTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGc                                                                                                                              <  5:90068/149‑1 (MQ=255)                      cTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATACGCTGGGGATGc                                                                                                                              >  6:343805/1‑149 (MQ=255)                          ccGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTc                                                                                                                           <  7:283632/148‑1 (MQ=255)                           cGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTg                                                                                                                         >  1:78088/1‑149 (MQ=255)                              atatGGGAACCTCCGCTGGAATCAAGGTTTCTATTGGTGTGGAATTGCTGATGGGGGACTATTATGTCCGCCAGCTAGGCGGCAGAGCGCAGTTATTTAATGACCCGCTGTTAACGCAATATATTAATTCTCTGGGGATGAGTCTgg                                                                                                                        <  2:106316/147‑1 (MQ=255)                                  gggAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCg                                                                                                                   >  1:49403/1‑148 (MQ=255)                                  gggAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGAACCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCg                                                                                                                   >  7:352595/1‑148 (MQ=255)                                   ggAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCa                                                                                                                 <  1:214461/149‑1 (MQ=255)                                      aCCTCCGCAGGAAGCACGCTTTCCATTTGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGc                                                                                                              <  6:199192/149‑1 (MQ=255)                                                  aGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTaa                                                                                                  <  7:279814/149‑1 (MQ=255)                                                    cACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGa                                                                                                <  1:77445/149‑1 (MQ=255)                                                         tttCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTa                                                                                                   <  7:60896/141‑1 (MQ=255)                                                         tttCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTa                                                                                                   >  8:60896/1‑141 (MQ=255)                                                             cATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCGATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGacac                                                                                             >  5:364408/1‑143 (MQ=255)                                                             cATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCGATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGacac                                                                                             <  6:364408/143‑1 (MQ=255)                                                             cATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAg                                                                                                 >  3:238177/1‑139 (MQ=255)                                                             cATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAg                                                                                                 <  4:238177/139‑1 (MQ=255)                                                              aTTGGTCAGGAAATGCAGATCGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCAGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCa                                                                                      >  4:113251/1‑149 (MQ=255)                                                               ttGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCAt                                                                                     <  2:267827/149‑1 (MQ=255)                                                                   tCAGGAAATGCAGAAGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACAACGTTTCAttttt                                                                                 >  2:344254/1‑149 (MQ=255)                                                                      ggAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTc                                                                              >  1:176532/1‑149 (MQ=255)                                                                      ggAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTc                                                                              <  5:343805/149‑1 (MQ=255)                                                                       gAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGttt                                                                                        >  1:282541/1‑138 (MQ=255)                                                                       gAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGttt                                                                                        <  2:282541/138‑1 (MQ=255)                                                                         aaTGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGa                                                                           <  3:96768/149‑1 (MQ=255)                                                                            gCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATc                                                                         >  7:73086/1‑148 (MQ=255)                                                                                 tGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCGATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCaacaac                                                                   >  4:282823/1‑149 (MQ=255)                                                                                    gCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGcc                                                                                                             <  4:264669/104‑1 (MQ=255)                                                                                    gCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGcc                                                                                                             >  3:264669/1‑104 (MQ=255)                                                                                    gCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAacgac                                                                >  5:208215/1‑149 (MQ=255)                                                                                       aCTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAAcgacga                                                              >  1:315592/1‑148 (MQ=255)                                                                                        cTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGaaa                                                            <  8:73086/149‑1 (MQ=255)                                                                                        cTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGaaa                                                            <  4:252474/149‑1 (MQ=255)                                                                                         tattatGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAAt                                                           >  3:32417/1‑149 (MQ=255)                                                                                           ttatGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAAtt                                                          >  4:119947/1‑148 (MQ=255)                                                                                           ttatGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAAtt                                                          <  1:344254/148‑1 (MQ=255)                                                                                             atGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAAc                                                       >  1:322451/1‑149 (MQ=255)                                                                                              tGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAAc                                                       <  5:297356/148‑1 (MQ=255)                                                                                              tGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAAc                                                       <  2:322451/148‑1 (MQ=255)                                                                                                tCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGc                                                     <  4:79584/148‑1 (MQ=255)                                                                                                  cGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTAAACGCAATATAATAATTCGCTGGCGAAGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCtt                                                  >  6:288621/1‑149 (MQ=255)                                                                                                     cAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGc                                               >  2:15332/1‑149 (MQ=255)                                                                                                        cTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGCTGCGTCTGGTTTCGCATGCGAATTCGGGTAAGCGACCGTTTCATGTTTGTGTGAGCACCACCGACGAAGTTGACGCCTGTGcttt                                            >  8:328688/1‑149 (MQ=255)                                                                                                        cTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGcttt                                            <  2:374622/149‑1 (MQ=255)                                                                                                         tACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCAttttttt                                                                               <  8:63938/113‑1 (MQ=255)                                                                                                           cgcgGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGctttct                                          <  4:279295/148‑1 (MQ=255)                                                                                                            gcgGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGctttcttt                                        <  2:49403/149‑1 (MQ=255)                                                                                                              ggaaGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTgg                                      <  5:329732/146‑1 (MQ=255)                                                                                                              ggCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTgg                                      <  3:190650/149‑1 (MQ=255)                                                                                                               gCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCAGGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTggc                                     <  5:11947/149‑1 (MQ=255)                                                                                                                      gcCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTAATCAACAACGACGAAATTAACGCCTTTGCTTTTTTTGGCGGCAACg                              >  8:13270/1‑149 (MQ=255)                                                                                                                        cGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACgtg                            >  8:60972/1‑149 (MQ=255)                                                                                                                         gTTAATTAATGACCCGCTGTTAACGCAATATATTTATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACgtgg                           <  3:113251/149‑1 (MQ=255)                                                                                                                         gTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACgtgg                           <  1:69038/149‑1 (MQ=255)                                                                                                                            aattaatGACCCGCTGTTAACGCGATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACgtggtg                         <  3:282823/148‑1 (MQ=255)                                                                                                                             attaatGACCCGCTGTTAACGCAATATTTTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAGCGCCTTTGCTTTGTTTGGCTGCAACGTGGtgct                       >  2:128478/1‑149 (MQ=255)                                                                                                                              ttaatGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACGTGGtgctg                      >  3:289988/1‑149 (MQ=255)                                                                                                                                 atGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTGAAGACACCGTTTCATTTTTTTCTTATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACGTGGTGCTGCa                    <  6:301870/148‑1 (MQ=255)                                                                                                                                      ccGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACGTGGTGCTGCACTCTGc              >  8:122479/1‑149 (MQ=255)                                                                                                                                           gTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACgtg                            >  2:34166/1‑130 (MQ=255)                                                                                                                                                  cAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACGTGGTGCTGCACTCTGCCCTGTTCCGTTa  >  2:32549/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GTTTGCCGACAGCGCAGACACCTTGCCGGATATGGGAACCTCCGCAGGAAGCACGCTTTCCATTGGTCAGGAAATGCAGATGGGCGACTATTATGTCCGCCAGCTACGCGGCAGCGCGCCGTTAATTAATGACCCGCTGTTAACGCAATATATTAATTCGCTGGGGATGCGTCTGGTTTCGCATGCCAATTCGGTTAAGACACCGTTTCATTTTTTTCTGATCAACAACGACGAAATTAACGCCTTTGCTTTCTTTGGCGGCAACGTGGTGCTGCACTCTGCCCTGTTCCGTTA  >  NC_000913/2616168‑2616461

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