BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,269	C→T	7.2%	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	7.2%	585.0 / 22.0	198	intergenic (+47/‑100)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base C (93/86); new base T (8/6); total (104/94)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

GTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCG  >  NC_000913/781120‑781417

gTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGGTTTAACATCAAAATAAGATGTTAAGCGGGAAGGAAAACGTTgcgc                                                                                                                                                        >  1:3390779/1‑150 (MQ=14)

  tCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCCGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTAAGATGTTAAGCGTGAAGGATAACGTTGCGCCa                                                                                                                                                      >  1:224614/1‑150 (MQ=255)

   cGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAg                                                                                                                                                     <  1:1573491/150‑1 (MQ=255)

   cGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGAGTTAAGGATAACGTTGCGCCAg                                                                                                                                                     >  1:3521502/1‑150 (MQ=255)

    gATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGc                                                                                                                                                    >  1:5179051/1‑150 (MQ=255)

     aTCCCGTCATCACCCACCACTCGGGTCGTTAGCACAGTTGGTAGATCAGTTTACTTTTAATAAATTTGTCGCAGGTTAGAATCCTGCACGAACCACCAATTTTAACAACAAAATAAGATGTTAATCGGGAAGTATAACGTTGCGACAgaa                                                                                                                                                   >  1:1170411/1‑148 (MQ=14)

      tCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCaa                                                                                                                                                  <  1:1868761/150‑1 (MQ=255)

      tCCCGTCATCACCCACCACTCGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCAGCACGACCCACCACGTTTAACATCCAACACGGTTGTTTAGCGTGAGGGTTAACGTTGCGACCGCaa                                                                                                                                                  >  1:2098465/1‑150 (MQ=11)

      tCCCGTCAGCA‑CCACCCCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAAc                                                                                                                                                 <  2:3715072/150‑1 (MQ=255)

        ccGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGAAGTTGAATTTTAATCAATTTGTCGCAGGTTCGCATCCTGCACGACCCACCAGTTTTAACATAAAACTAAGTAGATAAGAGTGATGGAAAACGTGTCGCcaacaaac                                                                                                                                                >  1:3349580/1‑148 (MQ=14)

         cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGTCTTTTAATCAATTTGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCAGATGTTAAGCGGGAAGGATAACGTTGCGCCAGCAACgg                                                                                                                                               >  1:1770098/1‑150 (MQ=255)

         cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCAGTTGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACgg                                                                                                                                               >  1:5488619/1‑150 (MQ=255)

          gTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGc                                                                                                                                              <  2:3115467/150‑1 (MQ=255)

           tcatcaCCCCCCCCGCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGcc                                                                                                                                             <  2:1478304/150‑1 (MQ=255)

               cccccaccACTCGGGTCGTTAGGTCAGTTGGTAGAGCTGTTGACTTTTAATCACTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTa                                                                                                                                         <  2:4694444/148‑1 (MQ=255)

                acccaccACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTTGTCGCAGGTTCGAATCCTGCACGACCCACCACTTTTAACATCAAAAAAAGGTGTTAAGCGTGAAGCATAACGTAGCCCCACCAACGGccccaag                                                                                                                                        >  1:4966501/1‑146 (MQ=14)

                acccaccACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTACCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAg                                                                                                                                        <  2:2448430/150‑1 (MQ=255)

                  ccaccaCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTgggg                                                                                                                                      >  2:611764/1‑150 (MQ=255)

                  ccaccaCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAggg                                                                                                                                      >  1:2839943/1‑150 (MQ=255)

                         tCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATCAAACTAAGAAGTTAAGCGTCAAGCAAAACGGTGCGCCAGCAACGGCCCCGAGCgcgag‑cga                                                                                                                              >  2:2694042/1‑150 (MQ=14)

                              tCGTTAGCACAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTATTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGag‑cgaag‑cga                                                                                                                        <  2:3620958/150‑1 (MQ=255)

                                      tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGag‑cgaag                                                                                                                            <  2:1916964/139‑1 (MQ=255)

                                      tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGag‑cgaag                                                                                                                            >  1:1916964/1‑139 (MQ=37)

                                         gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCATTTTAAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTCCGCCAGCAACGCCCCGTAGGCCTAG‑CGAAG‑AGAGTCaaacagcc                                                                                                             >  2:399511/1‑143 (MQ=14)

                                         gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCa                                                                                                             >  1:2696488/1‑150 (MQ=255)

                                          ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTAAGATGTTAAGCGTGAAGGATAACTTTGCGCCAGCAACGGCCCGTAGGGCGAG‑AGAAG‑CGAGTCATCCAGCAc                                                                                                            >  2:2884463/1‑150 (MQ=255)

                                            ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGa                                                                                                          >  2:512609/1‑150 (MQ=255)

                                             gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGAc                                                                                                         >  1:3719560/1‑150 (MQ=255)

                                                 agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTa                                                                                                                                         <  2:2745890/116‑1 (MQ=255)

                                                 agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCAGCCCGAccccc                                                                                                     >  2:3909004/1‑148 (MQ=255)

                                                  gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACccacc                                                                                                    >  1:1367853/1‑150 (MQ=255)

                                                  gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACccacc                                                                                                    >  1:1367834/1‑150 (MQ=255)

                                                  gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACcca                                                                                                      >  1:172044/1‑150 (MQ=14)

                                                     gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAAt                                                                                                 <  2:255097/150‑1 (MQ=255)

                                                      tttaattttAATCATTGGGTCGCTGTTTCGAATCCGGCACGACCCACCATTTTTAACAACAAACTCTTATTTTAAGCGTGATGGATAACGTTGCGCCACCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATc                                                                                                <  1:5082784/145‑1 (MQ=14)

                                                           ttttAATCATTTGGTCGCTGGTGCGAATCCTGCACGACCCACAATTTTTAACAACAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                                           <  1:2043209/150‑1 (MQ=17)

                                                                 tCAATTGGTCGCGGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTACAGATTgg                                                                                     >  1:5072270/1‑150 (MQ=255)

                                                                 tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                                                     >  1:3623800/1‑150 (MQ=255)

                                                                 tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                                                     >  1:2558687/1‑150 (MQ=255)

                                                                  cAATTGGGCGCTGTTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCATATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTTAAGATTggt                                                                                    <  1:1287053/150‑2 (MQ=14)

                                                                     ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACTTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                                                 <  2:4774373/150‑1 (MQ=255)

                                                                     ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                                                 <  2:1341443/150‑1 (MQ=255)

                                                                     ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                                                 <  2:3552828/150‑1 (MQ=255)

                                                                     ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                                                 >  1:1196253/1‑150 (MQ=255)

                                                                         tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGTAG‑CGAGTCAGACTGCACGACCCACCAAGCTTAGAGTGTGGCCCCGAGt                                                                             >  2:2631744/1‑150 (MQ=14)

                                                                         tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                                                             >  1:146132/1‑150 (MQ=255)

                                                                         tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGAGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                                                             >  2:2964245/1‑150 (MQ=255)

                                                                               gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCATAAAGATTGGCCCCGAGTaaaaaa                                                                       >  1:3285419/1‑149 (MQ=255)

                                                                                ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                                                      >  1:4527223/1‑150 (MQ=255)

                                                                                ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGGAGGGCGAG‑CGGAG‑CGAGTCATCCTGCACGACCCACCAAGCTTAGAGTTTGGCCCCGAGTAAAAATc                                                                      >  2:4909583/1‑150 (MQ=18)

                                                                                 tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                                                     >  1:804214/1‑150 (MQ=255)

                                                                                 tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                                                     >  1:2751659/1‑150 (MQ=255)

                                                                                 tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAACCTGAATGATTGGCCCCGCGTaaacaacc                                                                     >  1:3188224/1‑145 (MQ=17)

                                                                                   gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGCG‑CGAAG‑CGAGTCCTCCTGCACGACCCACCAATCTTAAAGAGTGGCCCCGAGTAAAAATCtgt                                                                   >  2:543912/1‑148 (MQ=255)

                                                                                   gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                                                   >  1:4883443/1‑150 (MQ=255)

                                                                                   gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                                                   <  2:4736470/150‑1 (MQ=255)

                                                                                    aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                                  <  2:5421502/150‑1 (MQ=255)

                                                                                     aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                                                 <  2:4196931/150‑1 (MQ=255)

                                                                                     aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                                                 <  2:3301178/150‑1 (MQ=255)

                                                                                       ccTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAgg                                                               <  1:1962558/150‑1 (MQ=255)

                                                                                         tGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                                                             >  2:960476/1‑150 (MQ=17)

                                                                                         tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTATCAGGTa                                                             >  1:4277054/1‑150 (MQ=255)

                                                                                          gCACGACCCACCAGTTTTAACATCAACCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                                            <  2:3798866/150‑1 (MQ=255)

                                                                                           cACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAac                                                           >  2:5216048/1‑150 (MQ=255)

                                                                                                 ccaccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGt                                                     >  1:473240/1‑150 (MQ=255)

                                                                                                  caccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGGAAAAATCTTTCAGGTAACACCCGTa                                                    >  2:4093442/1‑150 (MQ=255)

                                                                                                    ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                                                  >  1:614694/1‑150 (MQ=255)

                                                                                                    ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                                                  <  1:4650165/150‑1 (MQ=255)

                                                                                                          ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt                                            <  2:1099267/150‑1 (MQ=255)

                                                                                                           tAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                                                 <  2:232577/128‑1 (MQ=14)

                                                                                                           tAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                                                 >  1:232577/1‑128 (MQ=14)

                                                                                                            aaCATCTAACTCATATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGTCCGGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa                                          <  2:4816922/150‑1 (MQ=255)

                                                                                                             aCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAg                                         >  2:2017558/1‑150 (MQ=255)

                                                                                                               aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCCATCTAAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGcg                                       >  1:4986541/1‑149 (MQ=255)

                                                                                                               aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                                       <  2:5224252/150‑1 (MQ=255)

                                                                                                                tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAACACTGACGGGTAACACCCGAATGGGTCGTTAGCTc                                      >  2:4707525/1‑150 (MQ=21)

                                                                                                                tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTc                                      >  1:3037766/1‑150 (MQ=255)

                                                                                                                tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCAACCTGAACGACCCACAAATCAAACAGAATGGCCCCGACTAAAAATCTATCAGCAAACACACGAATGGGTCGTTAGCTc                                      >  2:4379346/1‑150 (MQ=18)

                                                                                                                 cacacccAGAAGTTAAGCGTGACGGATACCGTAGCGCCAGCTACGGCGCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa                                     <  2:2998259/144‑1 (MQ=25)

                                                                                                                   aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATATTTCAGGAAACACCCGTATGGGTCGTTAGCTCAGt                                   >  1:1303249/1‑150 (MQ=255)

                                                                                                                   aaCTCAGATGATAAGCGTGAAGGATAACGTTGCGCAAGCAACGGCCCGTAGGGCGAG‑AGAAG‑CGAGTCATCCAGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTAAGGTAACACACGTATGGGTCGTTATCTCAGt                                   >  1:3820475/1‑150 (MQ=25)

                                                                                                                    aatCATATGTAAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CTAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGtt                                  <  1:1168985/148‑1 (MQ=255)

                                                                                                                       cAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGtt                                  <  1:1969185/147‑1 (MQ=255)

                                                                                                                        aGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                                                 <  2:4348332/131‑1 (MQ=255)

                                                                                                                         gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTag                             <  2:412340/150‑1 (MQ=255)

                                                                                                                         gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTag                             <  2:4672990/150‑1 (MQ=255)

                                                                                                                            gTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGc                          >  1:4496051/1‑150 (MQ=255)

                                                                                                                            gTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGc                          >  1:4496286/1‑150 (MQ=255)

                                                                                                                             ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa                         >  1:1794592/1‑150 (MQ=255)

                                                                                                                              tAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAg                        >  1:1224423/1‑150 (MQ=255)

                                                                                                                               aaGCGTGAAGTATAACTTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACTACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt                       <  1:2431640/150‑1 (MQ=255)

                                                                                                                                 gCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGTGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTg                     <  2:4580299/150‑1 (MQ=255)

                                                                                                                                      aaGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                >  1:5368383/1‑150 (MQ=255)

                                                                                                                                       aGGATAACGTTGCGCCAGCAACGGCCCGTAGGgcgag                                                                                                                                  <  2:4073305/37‑1 (MQ=40)

                                                                                                                                        ggATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTa              <  1:960476/150‑1 (MQ=255)

                                                                                                                                              cGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAt        <  1:3063007/150‑1 (MQ=255)

                                                                                                                                               gTTGCGCCAGCTTCGGCCCGTAGGGCGCG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAtt       <  2:4564872/150‑1 (MQ=255)

                                                                                                                                               gTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGGAAGGGTCGTTAGCTCAGTTGGTAGAGCAGTTTGAATTTAATCaaat       >  2:120387/1‑148 (MQ=37)

                                                                                                                                                  gcgcCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGt    <  1:518896/150‑1 (MQ=255)

                                                                                                                                                    gcCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCg  <  1:4365498/150‑1 (MQ=255)

GTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAG‑CGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCG  >  NC_000913/781120‑781417

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