BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,132,550	T→C	G135G (GGT→GGC)	yghT →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,132,550	0	T	C	82.2%	640.7 / 130.7	326	G135G (GGT→GGC)	yghT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (32/25); new base C (143/121); total (176/150)
Fisher's exact test for biased strand distribution p-value = 8.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 334 total aligned reads displayed.

GTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAA  >  NC_000913/3132459‑3132639

gtcgCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑ggtg                                                                                          <  1:2667484/93‑1 (MQ=255)

 tcgCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑ggtgg                                                                                         >  1:2216265/1‑93 (MQ=255)

 tcgCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑ggcgg                                                                                         >  2:3584234/1‑93 (MQ=255)

 tcgCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑ggcgg                                                                                         >  2:5931020/1‑93 (MQ=255)

     gggTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCgc                                                                                     >  2:2990055/1‑92 (MQ=255)

     gggTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCgc                                                                                     >  2:6195600/1‑92 (MQ=255)

     gggTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCgc                                                                                     >  1:4164844/1‑92 (MQ=255)

         tAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTTGTTTTg                                                                                 <  2:6993197/93‑1 (MQ=255)

         tAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTg                                                                                 <  2:6056386/93‑1 (MQ=255)

         tAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTg                                                                                 >  2:1123193/1‑93 (MQ=255)

           aTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGtc                                                                               <  2:2273202/93‑1 (MQ=255)

           aTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGCGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGtc                                                                               >  1:4260485/1‑93 (MQ=255)

              gTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCa                                                                            <  1:7312725/93‑1 (MQ=255)

              gTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCa                                                                            <  2:1354875/93‑1 (MQ=255)

              gTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCa                                                                            <  2:1266412/93‑1 (MQ=255)

               ttCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCAt                                                                           >  1:6986016/1‑93 (MQ=255)

                 cTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCAtgg                                                                         >  2:5816222/1‑93 (MQ=255)

                     ccgaccgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCAtggctgg                                                                     <  1:803569/93‑1 (MQ=255)

                     ccgaccgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCAtggctgg                                                                     <  1:4614667/93‑1 (MQ=255)

                        accgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCaa                                                                  >  1:251860/1‑93 (MQ=255)

                          cgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGc                                                                >  2:2135403/1‑93 (MQ=255)

                           gTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGcc                                                               >  1:3885669/1‑93 (MQ=255)

                           gTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGcc                                                               <  2:4047192/93‑1 (MQ=255)

                           gTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGcc                                                               >  1:5208885/1‑93 (MQ=255)

                           gTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGcc                                                               <  2:2780053/93‑1 (MQ=255)

                            ttATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCa                                                              >  2:3169806/1‑93 (MQ=255)

                            ttATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCa                                                              <  2:7128778/93‑1 (MQ=255)

                             tATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAg                                                             >  2:2400696/1‑93 (MQ=255)

                               tCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGa                                                           >  2:2544420/1‑93 (MQ=255)

                                ccTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGaa                                                          <  2:4608525/93‑1 (MQ=255)

                                 cTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAAc                                                         <  1:1945658/93‑1 (MQ=255)

                                 cTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGGGTTCCCACCTAACGTTGAA‑GGCGGTAGCTTTGTCTCATGGCTGGCAAGCCAGGAAc                                                         >  1:5540920/1‑93 (MQ=255)

                                  tCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACg                                                        >  2:5730048/1‑93 (MQ=255)

                                  tCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACg                                                        >  2:6534219/1‑93 (MQ=255)

                                    aGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGGAAGTa                                                      >  1:320701/1‑93 (MQ=255)

                                    aGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTa                                                      <  2:220982/93‑1 (MQ=255)

                                    aGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTa                                                      <  2:3281487/93‑1 (MQ=255)

                                     ggACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTaa                                                     <  1:5289446/93‑1 (MQ=255)

                                     ggACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTaa                                                     <  1:198572/93‑1 (MQ=255)

                                     ggACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTaa                                                     <  2:1988378/93‑1 (MQ=255)

                                     ggACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTaa                                                     >  1:1618116/1‑93 (MQ=255)

                                      gACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTaaa                                                    <  1:840545/93‑1 (MQ=255)

                                       aCCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAg                                                   >  1:4385897/1‑93 (MQ=255)

                                       aCCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAg                                                   >  2:3896906/1‑93 (MQ=255)

                                        ccAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGc                                                  >  1:1309479/1‑93 (MQ=255)

                                         cAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCg                                                 >  1:2881594/1‑93 (MQ=255)

                                         cAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCg                                                 >  2:2946137/1‑93 (MQ=255)

                                          aaaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGt                                                >  1:6563520/1‑93 (MQ=255)

                                          aaaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGt                                                >  2:2639909/1‑93 (MQ=255)

                                          aaaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGt                                                >  1:5520773/1‑93 (MQ=255)

                                          aaaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGt                                                >  2:2746735/1‑93 (MQ=255)

                                          aaaTTCCTGGCGCTTACGATGGTACGGTGTCCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGt                                                <  1:3375237/93‑1 (MQ=255)

                                           aaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGtt                                               >  2:1758524/1‑93 (MQ=255)

                                           aaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGtt                                               >  2:7464029/1‑93 (MQ=255)

                                           aaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGtt                                               >  2:1815831/1‑93 (MQ=255)

                                           aaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGAAAGCCAGGAACGTAAAGCGtt                                               >  1:6479184/1‑93 (MQ=255)

                                            aTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGttt                                              <  1:6957333/93‑1 (MQ=255)

                                             ttCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTc                                             <  2:5062945/93‑1 (MQ=255)

                                             ttCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTc                                             >  2:425164/1‑93 (MQ=255)

                                             ttCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTc                                             <  2:4296825/93‑1 (MQ=255)

                                              tCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCa                                            <  2:5204771/93‑1 (MQ=255)

                                              tCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCa                                            <  1:1867473/93‑1 (MQ=255)

                                              tCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCa                                            <  2:647239/93‑1 (MQ=255)

                                               ccTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCAc                                           <  1:612132/93‑1 (MQ=255)

                                               ccTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCAc                                           <  1:4801519/93‑1 (MQ=255)

                                               ccTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCAc                                           <  1:4210461/93‑1 (MQ=255)

                                                cTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACt                                          <  1:4736126/93‑1 (MQ=255)

                                                cTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACt                                          >  1:4438897/1‑93 (MQ=255)

                                                cTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACt                                          <  1:7495925/93‑1 (MQ=255)

                                                 tGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTg                                         >  2:6391506/1‑93 (MQ=255)

                                                  ggCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTgg                                        >  2:6878953/1‑93 (MQ=255)

                                                  ggCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTgg                                        <  1:7204910/93‑1 (MQ=255)

                                                   gcgcTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGa                                       <  2:3325840/93‑1 (MQ=255)

                                                     gcTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATg                                     <  1:5048200/93‑1 (MQ=255)

                                                           gATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGc                               >  2:5303263/1‑93 (MQ=255)

                                                           gATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGc                               >  1:1198285/1‑93 (MQ=255)

                                                           gATGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGc                               >  1:4586531/1‑93 (MQ=255)

                                                            aTGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGcc                              >  2:6671715/1‑93 (MQ=255)

                                                            aTGGTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGcc                              <  1:5233771/93‑1 (MQ=255)

                                                               gTACGGTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATa                           <  2:4954845/93‑1 (MQ=255)

                                                                   ggTGTTCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGcc                       <  2:4808830/93‑1 (MQ=255)

                                                                   ggTGTCCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGcc                       <  1:3997713/93‑1 (MQ=255)

                                                                       ttCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGAt                   <  2:4194148/93‑1 (MQ=255)

                                                                        tCCCACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATc                  >  2:5592556/1‑93 (MQ=255)

                                                                          ccACCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTg                >  2:778142/1‑93 (MQ=255)

                                                                           ccccTAACGTTGAG‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTgg               <  1:5897480/91‑1 (MQ=255)

                                                                            aCCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGGGAGCCATAAGCCCGATCTGGt              >  1:3633043/1‑93 (MQ=255)

                                                                            aCCTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGt              <  2:7164904/93‑1 (MQ=255)

                                                                             ccTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtc             >  1:6336163/1‑93 (MQ=255)

                                                                              cTAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtca            >  1:4705290/1‑93 (MQ=255)

                                                                               tAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcat           >  1:6181850/1‑93 (MQ=255)

                                                                               tAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcat           <  2:1482422/93‑1 (MQ=255)

                                                                               tAACGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcat           >  2:3520904/1‑93 (MQ=255)

                                                                               ctaCGTTGAAGGGGGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtca            <  1:2350987/91‑1 (MQ=255)

                                                                                aaCGTTGAA‑GGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCTTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcatc          <  1:5599555/93‑1 (MQ=255)

                                                                                 aCGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcatca         >  2:7287281/1‑93 (MQ=255)

                                                                                 aCGTTGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcatca         >  1:1022906/1‑93 (MQ=255)

                                                                                    ttGAA‑GGCGTCCGCTTTGTCTCAGGGCTGGCAAGCCAGGAACAAAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAAc      <  1:6722029/93‑1 (MQ=255)

                                                                                    ttGAA‑GGCGGTCGCTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAAc      >  1:3757371/1‑93 (MQ=255)

                                                                                        a‑GGCGGTCGCTTTGTCTCATGGCGGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCaa  <  1:1459700/93‑1 (MQ=255)

GTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAA‑GGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAA  >  NC_000913/3132459‑3132639

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 17 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`