BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,005,601	+C	coding (237/282 nt)	b4568 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,005,600	1	.	C	99.8%	1704.7 / ‑4.9	498	coding (236/282 nt)	b4568	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (266/231); total (267/231)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 531 total aligned reads displayed.

CGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GG‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTG  >  NC_000913/4005508‑4005687

 gTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GG‑c                                                                                        >  2:1141579/1‑93 (MQ=255)

  tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GG‑cc                                                                                       >  2:3617158/1‑92 (MQ=255)

   aaTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GG‑cct                                                                                      >  1:6529004/1‑91 (MQ=255)

   aaTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GG‑cct                                                                                      >  2:5274699/1‑91 (MQ=255)

        tttGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTAtc                                                                                  >  1:4433895/1‑93 (MQ=255)

        tttGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTAtc                                                                                  >  1:1908135/1‑93 (MQ=255)

          tGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTAtctc                                                                                <  1:2256390/93‑1 (MQ=255)

             gCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTctgc                                                                             <  2:4573028/93‑1 (MQ=255)

             gCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTctgc                                                                             <  1:343344/93‑1 (MQ=255)

               tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTctgctg                                                                           >  2:6691410/1‑93 (MQ=255)

                  gATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGAt                                                                        <  2:2959309/93‑1 (MQ=255)

                  gATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGAt                                                                        >  2:1663594/1‑93 (MQ=255)

                   aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATg                                                                       >  1:5261608/1‑93 (MQ=255)

                    ttGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATgg                                                                      <  1:1470552/93‑1 (MQ=255)

                     tGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATggg                                                                     <  2:5876566/93‑1 (MQ=255)

                      gcgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGt                                                                    <  1:599475/93‑1 (MQ=255)

                      gcgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGt                                                                    <  2:3926288/93‑1 (MQ=255)

                       cgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGtt                                                                   >  1:2633819/1‑93 (MQ=255)

                         ctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTAc                                                                 <  1:1799857/93‑1 (MQ=255)

                          tgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTAcc                                                                <  2:792790/93‑1 (MQ=255)

                           gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCa                                                               <  2:1312660/93‑1 (MQ=255)

                           gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCa                                                               >  2:6401874/1‑93 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                             <  2:4024526/93‑1 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                             <  2:6975360/93‑1 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                             >  2:313812/1‑93 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                             >  1:5577464/1‑93 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                             <  2:1949597/93‑1 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                             <  2:2235570/93‑1 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTATGCTGGATGGGTTACCAGa                                                             <  1:7376803/93‑1 (MQ=255)

                              tGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGAt                                                            <  2:3762320/93‑1 (MQ=255)

                               ggCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGG‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATg                                                           >  1:1971348/1‑93 (MQ=255)

                                gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGc                                                          >  1:3024176/1‑93 (MQ=255)

                                 ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCt                                                         <  1:613543/93‑1 (MQ=255)

                                 ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCt                                                         >  2:688781/1‑93 (MQ=255)

                                 ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCt                                                         >  2:2843898/1‑93 (MQ=255)

                                  cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTa                                                        <  2:2156200/93‑1 (MQ=255)

                                  cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTa                                                        >  2:2492894/1‑93 (MQ=255)

                                  cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTa                                                        <  2:3880101/93‑1 (MQ=255)

                                   tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                       >  1:5775899/1‑93 (MQ=255)

                                   tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                       >  1:2454184/1‑93 (MQ=255)

                                   tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                       >  1:6916393/1‑93 (MQ=255)

                                     cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgt                                                     <  2:6208052/93‑1 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                    >  2:2334949/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                    <  2:1818693/93‑1 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                    >  1:1036178/1‑93 (MQ=255)

                                       tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtgg                                                   <  2:1006822/93‑1 (MQ=255)

                                        ttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggt                                                  >  1:5001712/1‑93 (MQ=255)

                                         tGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggtg                                                 >  2:4830232/1‑93 (MQ=255)

                                          gATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGc                                                >  2:5556421/1‑93 (MQ=255)

                                           aTTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCa                                               >  2:7300930/1‑93 (MQ=255)

                                           aTTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCa                                               >  2:3371485/1‑93 (MQ=255)

                                            ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                              >  2:7103609/1‑93 (MQ=255)

                                            ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                              <  2:5359185/93‑1 (MQ=255)

                                             tATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACt                                             <  1:5386395/93‑1 (MQ=255)

                                             tATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACt                                             >  1:1621422/1‑93 (MQ=255)

                                                cGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaa                                          >  2:1515509/1‑93 (MQ=255)

                                                 gAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaa                                         >  2:2074300/1‑93 (MQ=255)

                                                  aaaaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaa                                        <  2:2622451/93‑1 (MQ=255)

                                                      aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAaga                                    <  2:1428358/93‑1 (MQ=255)

                                                       aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGTTGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                   >  1:6720714/1‑93 (MQ=255)

                                                       aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                   <  1:3457081/93‑1 (MQ=255)

                                                       aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                   >  1:136585/1‑93 (MQ=255)

                                                       aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                   <  2:6489559/93‑1 (MQ=255)

                                                       aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                   >  1:5323008/1‑93 (MQ=255)

                                                       aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                   >  1:6268393/1‑93 (MQ=255)

                                                        tGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAgg                                  >  2:5684524/1‑93 (MQ=255)

                                                        tGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAgg                                  <  2:4447553/93‑1 (MQ=255)

                                                        tGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAgg                                  <  1:6461854/93‑1 (MQ=255)

                                                         ggCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGc                                 <  1:6985042/93‑1 (MQ=255)

                                                         ggCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGc                                 >  2:5169061/1‑93 (MQ=255)

                                                         ggCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGc                                 <  2:7193206/93‑1 (MQ=255)

                                                            cTGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt                              >  1:1120542/1‑93 (MQ=255)

                                                             tGGCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGtt                             >  1:2376486/1‑93 (MQ=255)

                                                               gCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTc                           <  1:784100/93‑1 (MQ=255)

                                                               gCTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTc                           <  1:7109484/93‑1 (MQ=255)

                                                                cTGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCt                          >  2:7483598/1‑93 (MQ=255)

                                                                 tGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTg                         >  1:7460089/1‑93 (MQ=255)

                                                                 tGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTg                         <  1:852444/93‑1 (MQ=255)

                                                                 tGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTg                         >  1:126277/1‑93 (MQ=255)

                                                                 tGCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGCGGTGCACTGAAAAAAGAGGCGGTTTCTg                         >  1:7201181/1‑93 (MQ=255)

                                                                  gCATACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGc                        >  1:2713409/1‑93 (MQ=255)

                                                                    aTACGCTGATTATGGTGGGCGGTGGG‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAc                      >  1:2766729/1‑93 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                     <  1:3666145/93‑1 (MQ=255)

                                                                        gCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgc                  <  1:5466515/93‑1 (MQ=255)

                                                                        gCTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgc                  >  2:4958848/1‑93 (MQ=255)

                                                                         cTGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcg                 >  2:2742636/1‑93 (MQ=255)

                                                                          tGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcgc                >  2:5067407/1‑93 (MQ=255)

                                                                          tGATTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcgc                <  1:4282791/93‑1 (MQ=255)

                                                                            aTTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCca              >  1:472999/1‑93 (MQ=255)

                                                                            aTTATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGGTTCTGCACCTGCGCca              >  1:4070212/1‑93 (MQ=255)

                                                                             ttATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcac             <  1:1679237/93‑1 (MQ=255)

                                                                              tATGGTGGGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcaca            >  2:4230492/1‑93 (MQ=255)

                                                                                tggtggGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg          >  2:6020269/1‑93 (MQ=255)

                                                                                tggtggGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg          <  2:6757314/93‑1 (MQ=255)

                                                                                 ggtggGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGt         <  1:3751881/93‑1 (MQ=255)

                                                                                   tggGCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCg       >  1:288847/1‑93 (MQ=255)

                                                                                    gggCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGa      <  1:4562712/93‑1 (MQ=255)

                                                                                      gCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGc    <  2:4702304/93‑1 (MQ=255)

                                                                                      gCGGT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGc    <  1:4814001/93‑1 (MQ=255)

                                                                                        ggT‑GGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTg  >  1:5080189/1‑93 (MQ=255)

CGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑GG‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTG  >  NC_000913/4005508‑4005687

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