BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,269	C→T	5.1%	intergenic (+47/‑100)	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,269	0	C	T	5.1%	1464.6 / 37.4	462	intergenic (+47/‑100)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base C (205/217); new base T (12/11); total (231/231)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 483 total aligned reads displayed.

GTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGC  >  NC_000913/781120‑781418

gTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAAGTGAACCGTGAAGGATAACGATgcgc                                                                                                                                                       >  2:4908761/1‑150 (MQ=255)

 ttCGATCCCGTCATCACCCACCACTCGGGTCGTTATCTCAGTTGGTAGAGCAGTTGACTTTTAAAAAATTGGTCGCAGGTTAGAATCCAGCACGACCCACCAGGTTTAACATAAAAATCAGAAGTTAAGCGTGAAGGATAACGGTGAGcc                                                                                                                                                      >  1:2192090/1‑150 (MQ=14)

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

    gATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGATGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTCGAATAATGCAAGACCCACCAGTTTTAACATAAAAATAAGATGTTAAGAGTGAAGGATAAAGTTGAGACAGc                                                                                                                                                   >  1:5442546/1‑150 (MQ=14)

             atcaCCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATGAAACGCGGATGTTAATAGTGAAGGAAAAACGTGAGACAGCAACGGcccc                                                                                                                                          >  1:5113707/1‑149 (MQ=14)

                    accacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgcg                                                                                                                                   >  2:8201130/1‑150 (MQ=255)

                     ccccccGGGTCCTCAGCCCAGCTGGTAGAGCAGACCACATTCCACCAATTGGTCGCATGTTCGACGCCTGCCCGACCCACCACATATAACAACAACACCATATGTTTAGCGTGCTGTATACCGTTGCGCCAGCAACGGCCCGTAGGgcga                                                                                                                                  <  2:2048458/145‑1 (MQ=2)

                      cacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgcgag                                                                                                                                 >  2:1379309/1‑150 (MQ=255)

                       acTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgcgagc                                                                                                                                <  2:5176045/150‑1 (MQ=255)

                        cTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGGTCGAATCCTGCACGACCCACCAGTGTTAACATCAAACTCAGAAGGTTCGCGGGACGGGTAACGTTGGGCCAGCCACGGCCCGAAGggcgggcg                                                                                                                               >  2:9907173/1‑150 (MQ=14)

                           gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaag                                                                                                                            <  2:8778588/150‑1 (MQ=255)

                           gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaag                                                                                                                            >  2:4532831/1‑150 (MQ=255)

                           gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaag                                                                                                                            >  2:5771714/1‑150 (MQ=255)

                              tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCGGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGTTAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                                                         <  2:2252831/150‑1 (MQ=255)

                                gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGt                                                                                                                       >  2:8971630/1‑150 (MQ=255)

                                    gCTCAGTTGGTAGAGCAGTTGACTTTTCATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACGTCAAACTCAGATTATACGCGGTAAGTAAAACGTTGCTCCAGCAACGGCCCGATGGGCGAGCGAAGCGAGTCATc                                                                                                                   >  2:9479925/1‑150 (MQ=14)

                                     cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATcc                                                                                                                  >  1:9859296/1‑150 (MQ=255)

                                      tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCt                                                                                                                 >  2:7267895/1‑150 (MQ=255)

                                        aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGc                                                                                                               <  2:1816387/150‑1 (MQ=255)

                                         gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                                                              >  1:2616084/1‑150 (MQ=255)

                                            ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGa                                                                                                           >  2:6493765/1‑150 (MQ=255)

                                              tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                                                                                                         <  2:785969/150‑1 (MQ=255)

                                              tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                                                                                                         >  2:7737537/1‑150 (MQ=255)

                                                 cgctgTTGACTTTTAATCATTGGGTCGCTGGTTCGAATCCTGCACGCCCCACCTGTTTTAACATCAAACTCAGATGTTATGGGTGATGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                                                                      <  1:6532906/146‑1 (MQ=14)

                                                    aGTTTATTTTTAATCAATTGGTCGCAGTTTCGAATCCTGCACGCCCCACCAGTTTTAACATCAAACTCAGATTTTATGCGTGAAGGATAAGGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa                                                                                                   <  1:9213596/150‑1 (MQ=14)

                                                     gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                                                                  >  1:2806988/1‑150 (MQ=255)

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

                                                      ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCAGAAGATAAGCGTGAAGGATAACGAAGCGCCACCAACGGCCCGAAGGGCGAGCGACGCGAGTCATCCTGCACGACCCACCAATc                                                                                                 >  2:6270747/1‑150 (MQ=14)

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

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

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

                                                                aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                                                       <  2:7075361/150‑1 (MQ=255)

                                                                 tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGGGAAGCGTGAAGGATAACGTGGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAgggtgg                                                                                      >  1:8121191/1‑145 (MQ=255)

                                                                     ttGGTCGCAGGGTCGAATCCTGCACGACCCACCTTTTTTAACTTCAAACTCCGTTTTTAAGCGTGAAGGATAAGGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                                                  <  1:1721545/150‑1 (MQ=17)

                                                                       ggTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAACGATTGGCCCCGa                                                                                <  1:7515507/150‑1 (MQ=255)

                                                                        gTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACGCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAg                                                                               >  1:4886867/1‑150 (MQ=255)

                                                                         tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                                                              >  1:9106111/1‑150 (MQ=255)

                                                                         tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                                                              >  1:1266636/1‑150 (MQ=255)

                                                                         tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTAATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                                                              >  2:1218406/1‑150 (MQ=255)

                                                                           gCAGGTTCGAATCCTGCACGACCCACCTTTTTTAACATCAAACTCAGATGTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaa                                                                            <  1:602702/150‑1 (MQ=255)

                                                                           gCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACGTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCACCGGCCCGTAGGGCGAGCGAAGCGGGTCAAACGGCACGACCCACCCATCGTAAAGATTGGCCCCGacaaa                                                                            >  2:2569690/1‑146 (MQ=17)

                                                                              ggTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaaa                                                                         >  1:7889578/1‑150 (MQ=255)

                                                                               gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                                                        >  1:8762980/1‑150 (MQ=255)

                                                                                   gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATTTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGACCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                                                    <  1:4512489/150‑1 (MQ=255)

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

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

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

                                                                                     cccccgCACGCCCCACCATTTATAACATCAATCGCCGATGTTATTGGTGATGGTTTGAGTTGCGCCCTGAAGTGTGCGTTTGGCGAGCGCAGGTAGGAATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                                                  <  1:9418210/145‑1 (MQ=255)

                                                                                     aTCCTGCACGACCCACCAGTTTTAACATCTAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                                                  <  2:2958471/150‑1 (MQ=255)

                                                                                     aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                                                  <  1:8971630/150‑1 (MQ=255)

                                                                                        cTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                                                               >  1:1651380/1‑150 (MQ=255)

                                                                                              gACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCATGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAcc                                                         >  1:7830863/1‑150 (MQ=255)

                                                                                                 ccaccaGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGt                                                      >  1:3129631/1‑150 (MQ=17)

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

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

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

                                                                                                    ccaGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                                                   <  2:1020919/150‑1 (MQ=17)

                                                                                                    ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                                                   >  2:3230806/1‑150 (MQ=255)

                                                                                                    ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                                                   <  2:53617/150‑1 (MQ=255)

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

                                                                                                      aGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATggg                                                 <  1:6802274/150‑1 (MQ=255)

                                                                                                       gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                                                <  1:6493765/150‑1 (MQ=255)

                                                                                                       gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                                                <  2:3992736/150‑1 (MQ=255)

                                                                                                       gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                                                <  2:4226545/150‑1 (MQ=255)

                                                                                                        ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTc                                               <  2:9650460/150‑1 (MQ=255)

                                                                                                          ttAACAACAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt                                             <  1:3673907/150‑1 (MQ=255)

                                                                                                            aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa                                           >  1:4773180/1‑150 (MQ=255)

                                                                                                            aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGGATGGGTCGTTa                                           >  2:7971237/1‑150 (MQ=255)

                                                                                                             aCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAg                                          >  2:4434911/1‑150 (MQ=255)

                                                                                                               aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAGAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                                        >  2:7398525/1‑150 (MQ=255)

                                                                                                               aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                                        <  1:1995660/150‑1 (MQ=255)

                                                                                                               aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                                        <  2:4678544/150‑1 (MQ=255)

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

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

                                                                                                                        aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa                               <  1:6350935/150‑1 (MQ=255)

                                                                                                                         gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTag                              >  1:4712678/1‑150 (MQ=255)

                                                                                                                           tGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTagag                            <  2:3129631/150‑1 (MQ=255)

                                                                                                                           tGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTagag                            >  1:8442400/1‑150 (MQ=255)

                                                                                                                           tGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTagag                            <  2:3671901/150‑1 (MQ=255)

                                                                                                                             ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAATATTGGCCCCCAGTAAAAATCTTTCTGGTACCACCCGTATGGGTCGTGAGCACGGTGGGTagagga                          >  2:5260449/1‑148 (MQ=21)

                                                                                                                               aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt                        <  1:7660073/150‑1 (MQ=255)

                                                                                                                               aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt                        <  2:6973989/150‑1 (MQ=255)

                                                                                                                               aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt                        <  2:7478102/150‑1 (MQ=255)

                                                                                                                                   gTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGAc                    >  1:3335494/1‑150 (MQ=255)

                                                                                                                                    tGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACt                   >  1:3475351/1‑150 (MQ=255)

                                                                                                                                     gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtg                  <  1:3995937/150‑2 (MQ=255)

                                                                                                                                     gAAGGAGCACGTTGAGCGTCCTACGGAGCGTAGGTCGTGGGGAGCGAGTCATCCTGGACTTCCCCACCCTCATAAAGATTGGGCCCGAGTTACGATCCCTCACTGAACCACCGTTAGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt                  <  2:2807978/150‑1 (MQ=255)

                                                                                                                                      taGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                                                              <  1:5765379/56‑1 (MQ=18)

                                                                                                                                      aaGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                 <  2:1600923/150‑1 (MQ=255)

                                                                                                                                      aaGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                 <  1:5234163/150‑1 (MQ=255)

                                                                                                                                       aGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtttt                <  1:1379309/150‑1 (MQ=255)

                                                                                                                                        ggATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTa               <  2:3447704/150‑1 (MQ=255)

                                                                                                                                         gATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa              <  1:2830098/150‑1 (MQ=255)

                                                                                                                                         gATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGACTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa              <  1:7273916/150‑1 (MQ=255)

                                                                                                                                           tAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATc            >  2:9364787/1‑150 (MQ=255)

                                                                                                                                             aCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa          <  2:9859296/150‑1 (MQ=255)

                                                                                                                                                ttGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTg       >  1:5911525/1‑150 (MQ=255)

                                                                                                                                                    gcCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCg   >  1:5555557/1‑150 (MQ=255)

                                                                                                                                                    gcCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCg   <  2:5513127/150‑1 (MQ=255)

                                                                                                                                                     ccAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGAATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTTGTCGc  >  1:2460222/1‑150 (MQ=255)

GTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGC  >  NC_000913/781120‑781418

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