BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,234	A→G	5.2%	intergenic (+12/‑135)	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,234	0	A	G	5.2%	593.2 / 11.4	194	intergenic (+12/‑135)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base A (104/71); new base G (7/3); total (112/82)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

GGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAG  >  NC_000913/781085‑781383

gggAGTTCAACTCCCTTACAAGGTTGGGGTGGGGGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCa                                                                                                                                                       <  1:753709/150‑1 (MQ=11)

gggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCa                                                                                                                                                       >  1:1740476/1‑150 (MQ=35)

gggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCa                                                                                                                                                       >  1:2306735/1‑150 (MQ=35)

  gagaGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCaaa                                                                                                                                                     >  2:3897838/1‑150 (MQ=37)

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

        aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGa                                                                                                                                               >  2:3550633/1‑150 (MQ=255)

          cTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATCAAACTAAGATg                                                                                                                                             >  1:3950159/1‑150 (MQ=255)

             ccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCCCGACCCACCAGTTTTAACATCAAACTCAGATGTTa                                                                                                                                          >  1:985023/1‑150 (MQ=255)

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

                   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCggg                                                                                                                                    >  2:2133832/1‑148 (MQ=255)

                    aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGa                                                                                                                                   <  1:2601121/150‑1 (MQ=255)

                          gggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACTAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATa                                                                                                                             >  2:267150/1‑150 (MQ=255)

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

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

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

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

                          gggTCGGCGGTTCGATCCCGTCATCACCCACCACACGGGTCGTTAGCGCAGTTTGTAGAGCAGTTGAATTTTTATAAAATGGTCGCAGGTTCGAATACAGCACCACCCACCACATTTAAAAACAAACACAGATGTTAAGAGTGAAGGATa                                                                                                                             >  1:1782952/1‑150 (MQ=11)

                           ggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATaa                                                                                                                            >  2:3716563/1‑150 (MQ=255)

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

                                  ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgcg                                                                                                                     >  1:3504912/1‑150 (MQ=255)

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

                                     tCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCa                                                                                                                  <  1:2056515/150‑1 (MQ=255)

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

                                           ccGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACGTCAAAATCAGAAGTTAAGCGTGAAGGAAAACTTTGAGACAGCacct                                                                                                            >  2:3788204/1‑147 (MQ=14)

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

                                            cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAGGTTAAGAGTGAAGGATAACGTGGCGCCAGCAACgg                                                                                                           >  2:3176427/1‑150 (MQ=255)

                                            cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCAGATGTTAAGCGTGAAGGATAACGTTGCGCCCGCAACgg                                                                                                           >  1:2887717/1‑150 (MQ=255)

                                                atcaCCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAAACCTGCACGACCCACCCGTTTTAAAATCAAACTCAGATGTTAAGCGTGAAGGATAACGTGGCGCCAGCAAAGGCCCg                                                                                                       >  1:3312235/1‑150 (MQ=14)

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

                                                         cacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATACTGCACGACCCACCAGTTTTAACATCAAAATAAGATGTTAAGCGTGAAGGATAACCTTTCGCCAGCAACGGCCCGGAGGgcgag                                                                                              >  1:2632427/1‑150 (MQ=14)

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

                                                                ctCGCTAGCTCAGTTGGTAGAGGAGTTGACTTTTCATCAATTGGTCGCTGGTCCGAATCCCGCACGACCAACCACTTATTACATCTAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcg                                                                                       <  2:2893165/149‑1 (MQ=14)

                                                                             ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCAc                                                                          >  2:2785156/1‑150 (MQ=255)

                                                                                       aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa                                                                >  2:3061119/1‑150 (MQ=255)

                                                                                        gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACGCAGATGTTAAGCGGCAAGGATAACGGTGCGCCAGCAACGGCCCGAAGGGCGAGCGAAGCGAGTCAACATGCCCGACCCACCAAt                                                               >  2:3228821/1‑150 (MQ=14)

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

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

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

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

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

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

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

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

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

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

                                                                                               tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAg                                                        >  1:641613/1‑150 (MQ=255)

                                                                                                  aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGAtt                                                     <  2:1928054/150‑1 (MQ=255)

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

                                                                                                        ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCAt                                                                                 >  1:2565394/1‑116 (MQ=255)

                                                                                                        ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCCAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCAt                                                                                 <  2:2565394/116‑1 (MQ=255)

                                                                                                        ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAAGGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                               <  1:3867452/150‑1 (MQ=255)

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

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

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

                                                                                                                aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCCAACTTAAAGATTGGCCCCGAGAaaaa                                       >  1:3971118/1‑150 (MQ=255)

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

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

                                                                                                                      gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTCTTAAAGATTGGCCCCGAGTAAAAATCttt                                 <  2:2765345/150‑1 (MQ=255)

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

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

                                                                                                                       aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                >  1:4018901/1‑150 (MQ=255)

                                                                                                                        aTCCTGCACGCCCCACCAGTTTTAACATCCAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                               <  1:2869011/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                              >  2:441821/1‑150 (MQ=255)

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

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

                                                                                                                             gCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                          <  2:1880607/150‑1 (MQ=14)

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

                                                                                                                             gCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                          <  2:1782952/150‑1 (MQ=255)

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

                                                                                                                             gCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCCTTCAGGTaa                          <  1:2760963/150‑1 (MQ=255)

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

                                                                                                                                cGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAacac                       <  1:2790655/150‑1 (MQ=255)

                                                                                                                                 gACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAATGATTGGCCCCGAGTAAATATCTTTCAGGTAACAcc                      >  1:658326/1‑150 (MQ=255)

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

                                                                                                                                    ccaccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCga                   >  2:2039287/1‑149 (MQ=255)

                                                                                                                                      accaGTTTTAACATCGAAGACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGAtt                                                     <  1:2373277/114‑1 (MQ=255)

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

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

                                                                                                                                      accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAACTTAAAGATTGGCCCCGCGTAAAAAACATTCAGGTACCACCCGGAt                 >  1:2676519/1‑150 (MQ=18)

                                                                                                                                      accaGATTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                 <  2:2907611/150‑1 (MQ=255)

                                                                                                                                       ccaGTTGTAACACCACCCGCTGTTGTTTAGCGTGAAGGATTACGTTGGGCCGGCATCGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                <  2:3589703/150‑1 (MQ=25)

                                                                                                                                        caGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg               >  2:2120872/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                                                                   tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTc    >  1:3241258/1‑150 (MQ=255)

                                                                                                                                                   tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTc    >  2:3073243/1‑150 (MQ=255)

                                                                                                                                                   tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTc    >  2:1455732/1‑150 (MQ=255)

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

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

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

GGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAG  >  NC_000913/781085‑781383

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