BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,238	T→A	8.1%	intergenic (+16/‑131)	lysY → / → lysZ	tRNA‑Lys/tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	781,238	0	T	A	8.1%	599.4 / 16.8	200	intergenic (+16/‑131)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base T (97/80); new base A (16/4); total (115/85)
Fisher's exact test for biased strand distribution p-value = 3.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

CACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGG  >  NC_000913/781092‑781387

cACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCACATGACTTTTAATCAATTGGTCGCAGGTTCGAATCATGCACGACCAACCAGTATTAACATAAAACTCAg                                                                                                                                                    >  2:3911189/1‑150 (MQ=14)

 aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGa                                                                                                                                                   <  1:288165/150‑1 (MQ=255)

   cTCACTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATg                                                                                                                                                 >  1:3241682/1‑150 (MQ=255)

       cTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTaa                                                                                                                                             <  2:2579083/150‑1 (MQ=255)

         tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGc                                                                                                                                           <  1:2994830/150‑1 (MQ=255)

            aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTg                                                                                                                                        >  2:2728674/1‑150 (MQ=255)

             aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGTGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGa                                                                                                                                       <  1:1754592/150‑1 (MQ=255)

             aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGTTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGa                                                                                                                                       <  1:2601195/150‑1 (MQ=255)

                  ggggTCGGCGGTTCGATCCCGTCAACAACCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGAt                                                                                                                                  <  2:3178767/150‑1 (MQ=255)

                   gggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATa                                                                                                                                 >  2:4204619/1‑150 (MQ=255)

                   gggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATa                                                                                                                                 >  2:3387425/1‑150 (MQ=255)

                   gggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATa                                                                                                                                 >  2:2760963/1‑150 (MQ=255)

                   gggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTTGTCGCAGGTTCGAATCAAGCAACACCCACCACATTTAAAATAAAAATAAGATGTTAAGCGTGAAGGATa                                                                                                                                 >  1:1085908/1‑150 (MQ=14)

                   gggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAAATGGTCGCAGGTTCGAATACAGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATa                                                                                                                                 >  1:1548885/1‑150 (MQ=14)

                   gggTCGGCGGTTCGATCCCGTCATCACCCACCACACGGGTCGTTAGCTCAGTTGGTAGAGCAGTTTACTTTTAATAAAATTGTCGCAGGTTCGAATACTGCACGACCCACCACATTTAAAATCAAACACAGATTTTAATCGTGAATGATa                                                                                                                                 >  1:1782959/1‑150 (MQ=11)

                         gcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTg                                                                                                                           >  2:4066826/1‑150 (MQ=255)

                            gTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCACAGTTGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAATAGTGAAGGATAAAGTTgcgc                                                                                                                        >  1:894474/1‑150 (MQ=14)

                                 aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCa                                                                                                                   >  1:3308113/1‑150 (MQ=255)

                                    ccGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACg                                                                                                                >  1:40788/1‑150 (MQ=255)

                                    ccGTCATCACCCACCACTCGGGTCGTTAGCTCAGGTGGTAGAGCCGTTGACTTTTAATAAATTTGTCGCAGGTTCGAATCCAGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGAGAAGGATAAAGTTGCGCCAGCAAcc                                                                                                                >  1:4203844/1‑149 (MQ=14)

                                    ccGTCATCACCCACCACTCGGGTCGTGAGCTCAGTTGGTAGTGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCACCACTTTTAACATCAAAAACAGATGTTAAGCGTGAAGGTTAAAGTTGGGACAGCaaag                                                                                                                >  1:3960020/1‑148 (MQ=11)

                                     cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAATCGCAGATGTTATGCGTGAGGGATCACGTTGCGCCAGCAACgg                                                                                                               >  2:192912/1‑150 (MQ=14)

                                         atcaCCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGGTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATCAAAAACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCg                                                                                                           >  1:1005765/1‑150 (MQ=255)

                                                     tCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACGTCTAACACAGTAGTGAACCGTAAGGAAAAACGTTCCGCCAGCACAGGCCCGTAGGgcgagcga                                                                                               >  2:1382246/1‑150 (MQ=11)

                                                         gTCGTTAGCTCATTTGGTAGAGCAGTTGACTTTTAATCAATTGGGCGCAGGTTCGACTCCTGCACGACCCACCAGTTCTCACATCCAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcg                                                                                           <  2:2802433/150‑1 (MQ=14)

                                                             ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                                       >  2:211423/1‑150 (MQ=255)

                                                             ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGAGAAGCGCGTc                                                                                       >  1:1926959/1‑150 (MQ=255)

                                                             ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATAAGATGTTAAGCGTGAAGTATAACTTTGCTCAAGCAACGCCCCGAAGGCCGATCGAAGCGAGta                                                                                       >  2:3020326/1‑149 (MQ=14)

                                                               agccccgTTGGTAGAGCAGCTGACTTCGAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCAt                                                                                     <  2:2130396/144‑1 (MQ=14)

                                                                     gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                               <  2:3787756/150‑1 (MQ=255)

                                                                         gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCAGCACGAc                                                                           >  1:742335/1‑150 (MQ=255)

                                                                         gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAAGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGAGAAGCGAGTCATCCTGCACGAc                                                                           >  1:996717/1‑150 (MQ=255)

                                                                                  ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATc                                                                  >  2:2901224/1‑150 (MQ=255)

                                                                                  ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACAAATc                                                                  >  2:1156937/1‑150 (MQ=255)

                                                                                     aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTCTTa                                                               >  1:1461570/1‑150 (MQ=255)

                                                                                       ttttAATCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             <  2:30679/150‑1 (MQ=255)

                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACTTCGAACACTGATGTTAAGCTTGAAGGATAACTTTGCGCCAGCAACGGCCCGAAGGGCGAGCGAAGAGTGTCAACCAGCACCACCCACCAATCATaaa                                                             >  1:4004869/1‑150 (MQ=11)

                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             >  1:2574717/1‑150 (MQ=255)

                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             >  1:1140261/1‑150 (MQ=255)

                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             >  2:58934/1‑150 (MQ=255)

                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             >  2:838050/1‑150 (MQ=255)

                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCAGCACGACCCACCCAACTaaaa                                                             >  1:1151423/1‑150 (MQ=255)

                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCt                                                                                  >  1:1298403/1‑129 (MQ=38)

                                                                                        tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCAGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAg                                                            >  1:645875/1‑150 (MQ=255)

                                                                                           aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATAAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGAta                                                         >  1:2232067/1‑149 (MQ=255)

                                                                                            aTCAATGGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTTAAGATTg                                                        <  2:3643990/150‑1 (MQ=255)

                                                                                              cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGc                                                      >  2:317002/1‑150 (MQ=255)

                                                                                                    gTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAg                                                >  2:4197959/1‑150 (MQ=255)

                                                                                                        cAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                            >  2:1900751/1‑150 (MQ=255)

                                                                                                        cAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                            <  1:2804066/150‑1 (MQ=255)

                                                                                                        cAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                            >  1:3078324/1‑150 (MQ=255)

                                                                                                           gTTCGAAGCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                         <  1:3810644/150‑1 (MQ=255)

                                                                                                            ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                        >  2:2216887/1‑150 (MQ=255)

                                                                                                              cgttctcTGCACGACCAACCAGTTTTAACATCTATGTCAGATTTTAAGCGTGTATTATAACGTTGCGCAAGCAACGGCCCGTAGGGAGAGCGAAGCGAGTCATCCTGCACGACCCACCAATATTAAAGATTGGCCCCGAGTAAAAATCtt                                      <  2:1227871/144‑1 (MQ=17)

                                                                                                              cGAATTCTGCACGACCCACCAGTTTTAACATCTAAGTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                      <  1:3326200/150‑1 (MQ=255)

                                                                                                               gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                     <  1:2103375/150‑1 (MQ=255)

                                                                                                               gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                     <  2:2769949/150‑1 (MQ=255)

                                                                                                               gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                     <  1:3459071/150‑1 (MQ=255)

                                                                                                               gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                     <  1:3184548/150‑1 (MQ=255)

                                                                                                               gAATCCTGCACGACCCACCAGTTTACACAACAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                     <  1:3928408/150‑1 (MQ=255)

                                                                                                                aaTCCTGCACGACCCACCAGTTTTAACATCGAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATc                                                                                    >  1:356564/1‑102 (MQ=255)

                                                                                                                aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATc                                                                                    <  2:356564/102‑1 (MQ=255)

                                                                                                                aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaaaacatac                                    >  2:1291456/1‑144 (MQ=255)

                                                                                                                aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                    <  1:2785156/150‑1 (MQ=255)

                                                                                                                 aTCCTGCACGACCCACCAGTTTTAACATCCAAGTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                   <  2:2448712/150‑1 (MQ=255)

                                                                                                                 aTCCTGCACGACCCACCAGTTTTAACATAAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                   <  1:3716563/150‑1 (MQ=255)

                                                                                                                  tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  >  2:1409571/1‑150 (MQ=255)

                                                                                                                  tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  <  1:3167128/150‑1 (MQ=255)

                                                                                                                  tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  <  1:234002/150‑1 (MQ=255)

                                                                                                                  tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  >  2:441821/1‑150 (MQ=255)

                                                                                                                  tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGGGAAGCGTGACGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGACGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  >  2:1327695/1‑150 (MQ=255)

                                                                                                                   ccTGCACGACCCACCAGTTTTAACATCATACTCAGATGTTAAGCGTGAAGGATATCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAgg                                 <  2:2354331/150‑1 (MQ=255)

                                                                                                                      gCACGACCCACTAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                              <  1:267150/150‑1 (MQ=255)

                                                                                                                      gCACGACCCACCTGTTTTAACATCAAACACAGATGTTTTGGGTTTATGATATGGTTGCGCCATCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                              <  2:1782959/150‑1 (MQ=18)

                                                                                                                      gCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                              <  2:1044796/150‑1 (MQ=17)

                                                                                                                      gCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                              <  2:1548885/150‑1 (MQ=255)

                                                                                                                       cACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAac                             <  2:3591892/150‑1 (MQ=255)

                                                                                                                        aCGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTAAAAGATTGGCCCCGAGTAAAAATCTTGCAGGTAaca                            >  2:4119703/1‑150 (MQ=255)

                                                                                                                             ccaccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTCAAAATCTTTCAGGTAACACCCgg                       >  2:2279238/1‑149 (MQ=255)

                                                                                                                             ccaccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCATAAAGATTGGCCCCGAGTAAAAAGCGGTCAGGTAACAACCga                       >  2:2039272/1‑149 (MQ=18)

                                                                                                                               accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                     <  2:4095524/150‑1 (MQ=255)

                                                                                                                               accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                     >  1:2676514/1‑150 (MQ=255)

                                                                                                                               accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGGAt                     >  1:594210/1‑150 (MQ=255)

                                                                                                                                 caGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                   <  2:4018901/150‑1 (MQ=255)

                                                                                                                                 caGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                   <  2:3856645/150‑1 (MQ=255)

                                                                                                                                   gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                 >  1:4128916/1‑150 (MQ=255)

                                                                                                                                     tttAAAATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCg               <  1:1425453/150‑1 (MQ=255)

                                                                                                                                      ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt              <  2:879425/150‑1 (MQ=255)

                                                                                                                                      ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt              <  1:4132742/150‑1 (MQ=255)

                                                                                                                                      ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt              <  1:1954544/150‑1 (MQ=255)

                                                                                                                                      ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt              <  1:2901224/150‑1 (MQ=255)

                                                                                                                                        aaCATCCAAGTCAGATGTTAAGCGTGCAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa            <  2:2025082/150‑1 (MQ=255)

                                                                                                                                        aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa            <  2:3021893/150‑1 (MQ=255)

                                                                                                                                        aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa            <  1:2216887/150‑1 (MQ=255)

                                                                                                                                             cAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa       <  1:558531/150‑1 (MQ=255)

                                                                                                                                              aaaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAg      >  1:4027606/1‑150 (MQ=255)

                                                                                                                                              aaaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAg      >  1:2731454/1‑150 (MQ=255)

                                                                                                                                                aCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGtt    <  2:1436950/150‑1 (MQ=255)

                                                                                                                                                  tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg  >  1:2436751/1‑150 (MQ=255)

                                                                                                                                                  tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg  >  2:2508073/1‑150 (MQ=255)

CACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGG  >  NC_000913/781092‑781387

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