BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,392,582	A→C	7.8%	intergenic (+35/‑1)	glyX → / → glyY	tRNA‑Gly/tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,582	0	A	C	7.8%	460.9 / 3.3	165	intergenic (+35/‑1)	glyX/glyY	tRNA‑Gly/tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base A (90/57); new base C (9/8); total (100/65)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

CAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAC  >  NC_000913/4392434‑4392729

cAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGTAGGGGGACGCGGTTACGGTGAACGTTTCCCGCTCAAAAATTTGTAAGTGCTGTGAGTCAGAGACCACCCAAg                                                                                                                                                    >  2:2729482/1‑150 (MQ=17)

 aaaaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGc                                                                                                                                                   <  1:1099709/150‑1 (MQ=255)

   aaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCGCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAGAGTTGAAAGTGCTGTACGGCACAGACCACCCAAGCgg                                                                                                                                                 >  1:1316809/1‑150 (MQ=255)

   aaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCgg                                                                                                                                                 <  1:595709/150‑1 (MQ=255)

      ttGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCAAATATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGaa                                                                                                                                              <  1:1386553/150‑1 (MQ=255)

        gAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATa                                                                                                                                            <  1:3091386/150‑1 (MQ=255)

         aaaaGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAg                                                                                                                                           >  2:133308/1‑150 (MQ=255)

         aaaaGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGAGCTGTAAGGCACAGACCACCCACGCGGGAATAg                                                                                                                                           >  2:3689900/1‑150 (MQ=255)

            aGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGAGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTc                                                                                                                                        >  2:1885406/1‑150 (MQ=255)

             gTGCTGCAACGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCa                                                                                                                                       <  2:3361757/150‑1 (MQ=255)

              tgctgcAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAGTTGAAAGTGCTGTAAGGCACAGACCACCCACGCGGGAAGAGCTCAg                                                                                                                                      >  1:2822053/1‑150 (MQ=255)

              tgctgcAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGGTCGAGTCTCGTTTCCCGCTCCCAAATGTGCAAGTGCTGTAAGGCCCAGACCACCCCAGCGGGAATAGCGCAg                                                                                                                                      >  1:1448144/1‑150 (MQ=18)

                 tgcAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGCAAGTGCTGTAAGGCACAGACCACCCAAGCGGGCATAGCGCTGTTg                                                                                                                                   >  2:635480/1‑150 (MQ=255)

                  gcAAATCACAGACCACCCAAGCGGTAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGATCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTgg                                                                                                                                  <  2:446551/150‑1 (MQ=255)

                    acaGCACAGACCACCCAAGCGGGAATAGCTCAGGTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTGTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTa                                                                                                                                <  1:3884439/148‑1 (MQ=255)

                    aaaTCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTa                                                                                                                                <  1:3312903/150‑1 (MQ=255)

                    aaaGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTa                                                                                                                                <  1:3315281/150‑1 (MQ=255)

                      aGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTaga                                                                                                                              <  2:718009/150‑1 (MQ=255)

                        cacaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGc                                                                                                                            >  1:3470081/1‑150 (MQ=255)

                         acaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCa                                                                                                                           <  1:1917734/150‑1 (MQ=255)

                           aGACCACCACCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCATGGTCGGGGTCGCTAGTTCGAGTCTCGTTTCCCGCACCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACg                                                                                                                         <  2:3914533/150‑1 (MQ=14)

                            gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGa                                                                                                                        >  1:2660215/1‑150 (MQ=255)

                               caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCACTCGGGAATAGCTCAGTTGGTAGAGCACGAcca                                                                                                                     >  2:1644247/1‑149 (MQ=255)

                               caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATCTGAAAGGGCTGTAAGGCACAGACCGCCCACGCGCGAATAGCGCAGTTGGTAGAGCCCGAccg                                                                                                                     >  2:3218333/1‑149 (MQ=14)

                                accCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCCAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGAGCTGTAAGGCACAGACCACCCACGCGGGAATAGCGCAGTTGGTCGAGCACGACCtt                                                                                                                    >  2:3020279/1‑150 (MQ=14)

                                accCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCtt                                                                                                                    >  1:3854132/1‑150 (MQ=255)

                                 cccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTGAAAGTGCTGTAAGGCCCAGACCACCCCCGCGGGAATAGCTCAGTGGGTAGTGCCCGACCGTg                                                                                                                   >  1:881489/1‑150 (MQ=14)

                                         ggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTc                                                                                                           >  2:2254284/1‑150 (MQ=38)

                                           aaTAGCTCTTTTGGTTGAGCACGACCGTGCCATGGTTGGGGTGGCGAGTTCGAGTCTCGTTTCCCGCTCAAAAATTGTAATGTGCTGTATGGCACAGACCACCCACGGGGGATTAGCTATTTTGGTAGCGCACTACCTTGCCACGGTCgg                                                                                                         <  2:3681606/150‑1 (MQ=2)

                                           aaTAGCTCCGTTTGTAGAGCACTACCTTGCCAAGGTCGGGGTCGCGAGTTCGATACTCGTTTCCCCCTCCAAAATTTTAAAGTGCTGAAAGGCACAGACCACCCAAGCGGTAATAGCTCAGTTGGTAGAGCACTACCTTGGCAAGGTCgt                                                                                                         <  1:1092933/150‑2 (MQ=11)

                                             tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATTAAAGACATGAAATAAACACACCACCCAACAGCGAATATCAGAGTTGGTAGAGCTCGAACGTGCCAAGGTCgggg                                                                                                       >  2:3662757/1‑150 (MQ=11)

                                                cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGAGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCATGCCAAGGTCGGGGTcg                                                                                                    >  1:2225129/1‑150 (MQ=255)

                                                 cctgTTGGTAGCGCCCGACCTTGCCAAGGTGGGGTTGGCTAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTGAAAGTGCTGTAAGGCACAGACCAACCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                   <  2:1370167/147‑1 (MQ=11)

                                                   aGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                 <  1:2069488/150‑1 (MQ=38)

                                                   aGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                 >  2:3252123/1‑150 (MQ=38)

                                                   aGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATGAAAGGGCTGAAAGGCACACACCCCCCAAGCGCGAATAGCGCAGTTGGTAGAGCACGACCAGGCCACGGGCGCGGTCGCGa                                                                                                 >  2:3041445/1‑150 (MQ=11)

                                                    gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTagag                                                                                                                             <  1:3126284/121‑1 (MQ=40)

                                                     ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAATGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGt                                                                                               >  1:2469080/1‑150 (MQ=25)

                                                       ggTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                                             <  2:3585435/150‑1 (MQ=38)

                                                          agagCACGACCTTGCCACGGTCGGGGTCGCGAGTTCGCGGGTCGTTTCCACCTCCAAATTATGAAAGTGCTGTAAGTCACAGACCACCAACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                                          <  2:2725982/150‑1 (MQ=11)

                                                           gagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAAGTGCTGTAATTCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                         >  2:1373927/1‑150 (MQ=255)

                                                             gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtct                                                                                       >  2:1584223/1‑150 (MQ=38)

                                                              cccGACCTTGGCAAGGTCGTGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAATATTTGAAAGTGCTGGAAGGCCCAGGCCACCCACGCGGGAATAGCGCAGTGGGTAGCGCACGTCCGTGCCACGGGCGGGGTCGCGCGTTCGGGtctg                                                                                      >  2:475083/3‑149 (MQ=11)

                                                              cACGACCTTGCCAAGGTCGTGGTCGCGAGTTCGAGTCACGATTCCCGCGCCAAAATTTGAAAGAGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCACAGTTGGTAGAGCACCACCATGCCAAGGACGGGGTCGCGAGTTCGAGtctc                                                                                      >  1:1107771/1‑150 (MQ=11)

                                                              cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                      <  1:1834542/150‑1 (MQ=38)

                                                              cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                      >  2:3923672/1‑150 (MQ=38)

                                                              cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTACCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                      <  1:213029/150‑1 (MQ=38)

                                                              cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGATTCCCGCGCCAAAAGTTGAAAGTGCTGTAAGGCACACACCACCCAAGCGGGAATCGCGCCGGGGGTAGAGCACGACCTTGCCCAGGTCGGGGTCGCGCGTTCGAGtctc                                                                                      >  1:794022/1‑150 (MQ=11)

                                                               aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCg                                                                                     <  2:3576060/150‑1 (MQ=38)

                                                                   ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                 <  1:3349143/150‑1 (MQ=38)

                                                                       gCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGc                                                                             >  1:1276371/1‑150 (MQ=38)

                                                                         cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATATCTCAGTTGGTAGAGAAAGACCTTGACAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                           >  2:346991/1‑150 (MQ=255)

                                                                         cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                           >  1:1191787/1‑150 (MQ=38)

                                                                         cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                           <  1:3421151/150‑1 (MQ=38)

                                                                         cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACATTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCACGCTc                                                                           >  2:2153781/1‑150 (MQ=255)

                                                                            ggTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaa                                                                        >  2:2141745/1‑150 (MQ=38)

                                                                             gTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCACCAAAATTTTAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAACGTCGGGGTCGCGAGTTCGAGTCTCGTTTCACGCTCCaaa                                                                       >  1:1349700/1‑150 (MQ=255)

                                                                                 gggTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttct                                                                   >  2:1521513/1‑150 (MQ=255)

                                                                                      gctaGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctct                                                              <  2:2374709/147‑1 (MQ=255)

                                                                                            tCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaa                                                        >  2:1483095/1‑150 (MQ=255)

                                                                                            tCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaa                                                        >  1:1160672/1‑150 (MQ=255)

                                                                                              gAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaaaa                                                      >  1:3113042/1‑150 (MQ=255)

                                                                                              gAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaaaa                                                      >  2:1085318/1‑150 (MQ=255)

                                                                                                gTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaaaaaa                                                    >  1:699091/1‑148 (MQ=255)

                                                                                                 tctcGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAatat                                                   >  2:1488197/1‑150 (MQ=255)

                                                                                                 tctcGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAatat                                                   >  2:1998428/1‑150 (MQ=255)

                                                                                                      tttCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCcaca                                              >  2:3936232/1‑150 (MQ=255)

                                                                                                       ttCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAg                                             >  2:1793261/1‑150 (MQ=255)

                                                                                                        tCCCTCTCCAAAATTTTAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTATGTTGGTTGAGCACGACCTTGCCATGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTATAAATAAATTATCCACAGc                                            <  1:1885406/150‑1 (MQ=17)

                                                                                                         cccGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCg                                           >  2:694452/1‑150 (MQ=255)

                                                                                                           cGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAc                                         >  2:2785147/1‑150 (MQ=255)

                                                                                                           cGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAc                                         >  1:621046/1‑150 (MQ=255)

                                                                                                            gCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAcg                                        >  1:1658047/1‑150 (MQ=255)

                                                                                                                 aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATg                                   >  2:2406967/1‑150 (MQ=255)

                                                                                                                  aaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGc                                  >  1:1774573/1‑150 (MQ=255)

                                                                                                                  aaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGc                                  >  1:1761551/1‑150 (MQ=255)

                                                                                                                  aaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGc                                  >  1:3731387/1‑150 (MQ=255)

                                                                                                                        gAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCaa                                                           >  2:2190444/1‑119 (MQ=255)

                                                                                                                        gAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCaa                                                           <  1:2190444/119‑1 (MQ=255)

                                                                                                                        gAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTAtt                            >  2:3315588/1‑150 (MQ=255)

                                                                                                                           aGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCt                         <  1:1483095/150‑1 (MQ=255)

                                                                                                                             tGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTgg                       >  2:574257/1‑150 (MQ=255)

                                                                                                                             tGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTgg                       >  2:1123021/1‑150 (MQ=255)

                                                                                                                               cTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtt                     >  2:3801496/1‑150 (MQ=255)

                                                                                                                                  tAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                  >  2:1768895/1‑150 (MQ=255)

                                                                                                                                  tAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                  >  1:2418334/1‑150 (MQ=255)

                                                                                                                                  caaGGCACAGACCACCCCAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                  <  1:3824788/149‑1 (MQ=255)

                                                                                                                                   aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttg                 <  1:3662757/150‑1 (MQ=255)

                                                                                                                                   aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttg                 <  1:3606331/150‑1 (MQ=255)

                                                                                                                                      gCACAGACCACCCAAGCGGGAATAGCTCCGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAACGTCTTCTCTCAATAAAATATCCACAGCGCCGCGATGCGGTCTTGCTGGTTTttgtgg              >  1:548331/1‑148 (MQ=18)

                                                                                                                                      gCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttgttg              >  1:1377925/1‑150 (MQ=255)

                                                                                                                                        acaGACCTCCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtc            <  2:3368066/150‑1 (MQ=255)

                                                                                                                                        acaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtc            <  2:1891140/150‑1 (MQ=255)

                                                                                                                                         caGACCCCCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTGTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtct           <  2:1816520/150‑1 (MQ=255)

                                                                                                                                         caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtct           <  1:2963696/150‑1 (MQ=255)

                                                                                                                                           gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTctct         <  1:916519/150‑1 (MQ=255)

                                                                                                                                               accCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACa     >  1:1304155/1‑150 (MQ=255)

                                                                                                                                                cccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaa    >  1:2475265/1‑150 (MQ=255)

                                                                                                                                                 ccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaaa   >  1:1812346/1‑150 (MQ=255)

                                                                                                                                                  cAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAc  >  2:1228346/1‑150 (MQ=255)

CAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAC  >  NC_000913/4392434‑4392729

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