BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,013,177	+C	coding (543/999 nt)	ysgA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,013,177	1	.	C	99.5%	769.0 / ‑3.8	214	coding (543/999 nt)	ysgA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (112/101); total (113/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 220 total aligned reads displayed.

TTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCCG‑GCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAC  >  NC_000913/4013092‑4013260

ttAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGGCGAGAATCGGCGCGTTa                                                                               >  1:5731797/1‑93 (MQ=255)

  agcagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAg                                                                             <  2:5151982/93‑1 (MQ=255)

    cagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGAt                                                                           <  1:5008949/93‑1 (MQ=255)

     agcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATc                                                                          <  1:1184297/93‑1 (MQ=255)

      gcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCg                                                                         >  2:1641810/1‑93 (MQ=255)

      gcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCg                                                                         <  2:5864933/93‑1 (MQ=255)

       cccGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCgg                                                                        >  1:1115129/1‑93 (MQ=255)

       cccGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCgg                                                                        >  2:228179/1‑93 (MQ=255)

       cccGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCgg                                                                        <  1:3100473/93‑1 (MQ=255)

        ccGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGt                                                                       >  2:4585038/1‑93 (MQ=255)

         cGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGtt                                                                      >  2:5316778/1‑93 (MQ=255)

         cGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGtt                                                                      <  1:2979991/93‑1 (MQ=255)

          gCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTg                                                                     >  2:5955879/1‑93 (MQ=255)

           cAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGc                                                                    >  1:5878407/1‑93 (MQ=255)

           cAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGc                                                                    <  2:539669/93‑1 (MQ=255)

             gcgcCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGa                                                                  <  2:3575168/93‑1 (MQ=255)

               gcCTGGCGCATGTTTTCACCCTTCTCCTGCGGAATGCTGTTATCCTGCCCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGata                                                                <  1:2410790/93‑1 (MQ=255)

               gcCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGata                                                                <  1:4872544/93‑1 (MQ=255)

                ccTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGatat                                                               <  2:4979792/93‑1 (MQ=255)

                ccTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGatat                                                               >  2:586273/1‑93 (MQ=255)

                 cTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATc                                                              >  1:1822659/1‑93 (MQ=255)

                  tGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCa                                                             >  1:4585146/1‑93 (MQ=255)

                    gcgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAAc                                                           <  2:1941545/93‑1 (MQ=255)

                    gcgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAAc                                                           <  2:95904/93‑1 (MQ=255)

                      gcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAg                                                         >  2:6034990/1‑93 (MQ=255)

                      gcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAg                                                         >  1:5960750/1‑93 (MQ=255)

                      gcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAg                                                         >  2:262397/1‑93 (MQ=255)

                      gcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAg                                                         >  2:3630158/1‑93 (MQ=255)

                       cATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAgg                                                        <  1:1452262/93‑1 (MQ=255)

                        aTGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgt                                                       <  1:4304564/93‑1 (MQ=255)

                         tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                      >  1:3920224/1‑93 (MQ=255)

                         tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                      >  2:5544064/1‑93 (MQ=255)

                         tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                      >  2:3720462/1‑93 (MQ=255)

                            tttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCACTAGgttggt                                                   >  1:1624778/1‑91 (MQ=255)

                             ttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGttt                                                  <  2:4233804/93‑1 (MQ=255)

                             ttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGttt                                                  <  1:2940241/93‑1 (MQ=255)

                              tCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTc                                                 <  1:6087827/93‑1 (MQ=255)

                                 aCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                                              >  2:3412803/1‑93 (MQ=255)

                                 aCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                                              <  1:4169830/93‑1 (MQ=255)

                                    ctctCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGaa                                           >  1:2164466/1‑93 (MQ=255)

                                     cctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAt                                          <  1:2448043/92‑1 (MQ=255)

                                      ctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                                         <  2:3326606/93‑1 (MQ=255)

                                      ctcCTGCGGAACGCTGTTATCCTGGCCGCCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                                         <  1:1095943/93‑1 (MQ=255)

                                       tcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTc                                        >  1:1538954/1‑93 (MQ=255)

                                       tcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTc                                        >  1:5352689/1‑93 (MQ=255)

                                       tcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTc                                        <  2:3427817/93‑1 (MQ=255)

                                        ccTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCa                                       >  2:668920/1‑93 (MQ=255)

                                        ccTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCa                                       >  2:1424699/1‑93 (MQ=255)

                                        ccTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCa                                       <  2:3555944/93‑1 (MQ=255)

                                         cTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAg                                      >  2:4308571/1‑93 (MQ=255)

                                          tGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGc                                     >  2:3033345/1‑93 (MQ=255)

                                           gCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCg                                    <  2:2368425/93‑1 (MQ=255)

                                            cGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGa                                   <  1:5653005/93‑1 (MQ=255)

                                             ggAATGCTGTTATCCTGACCATCATATAAGCCGAGAATCGGCGCGTTAAGATCTGTTGCTATATCAACAGGTTGTTTCGGTGAAGTCAGCGAc                                  >  1:5480763/1‑93 (MQ=255)

                                             ggAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGAc                                  >  1:850640/1‑93 (MQ=255)

                                             ggAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGAc                                  >  1:180282/1‑93 (MQ=255)

                                              gAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACt                                 <  2:2645868/93‑1 (MQ=255)

                                               aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACtt                                <  2:1992835/93‑1 (MQ=255)

                                                aTGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTg                               >  2:1070006/1‑93 (MQ=255)

                                                 tGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGt                              <  1:5267999/93‑1 (MQ=255)

                                                 tGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGt                              <  1:3067483/93‑1 (MQ=255)

                                                 tGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGt                              >  2:5270280/1‑93 (MQ=255)

                                                 tGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGt                              <  2:5107974/93‑1 (MQ=255)

                                                 tGCTGTTATCCTGACCACCATATAAGACGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGt                              >  1:2048365/1‑93 (MQ=255)

                                                  gCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGATGTTTCGGTGAATTCAGCGACTTGTc                             >  1:3353080/1‑93 (MQ=255)

                                                   cTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCg                            <  1:3616102/93‑1 (MQ=255)

                                                    tGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGc                           >  2:3119586/1‑93 (MQ=255)

                                                    tGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGc                           >  2:1728629/1‑93 (MQ=255)

                                                    tGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGc                           >  2:2783334/1‑93 (MQ=255)

                                                    tGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGc                           >  1:3752426/1‑93 (MQ=255)

                                                    tGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGc                           >  2:4757786/1‑93 (MQ=255)

                                                     gTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGcc                          >  2:3876111/1‑93 (MQ=255)

                                                      ttATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGccc                         >  2:2197267/1‑93 (MQ=255)

                                                       tATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCg                        >  1:5990185/1‑93 (MQ=255)

                                                       tATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCg                        >  1:4667985/1‑93 (MQ=255)

                                                        aTCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGt                       >  1:5031874/1‑93 (MQ=255)

                                                         tCCTGACCACCATATAAGGCGAGAAGCGGCGCGTTAAGATCGGGTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTc                      <  2:4733300/93‑1 (MQ=255)

                                                          ccTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCa                     >  2:5989354/1‑93 (MQ=255)

                                                            tGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGt                   >  1:3835483/1‑93 (MQ=255)

                                                            tGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGt                   >  1:785597/1‑93 (MQ=255)

                                                            tGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGt                   >  1:3199945/1‑93 (MQ=255)

                                                             gACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGtt                  <  1:1689820/93‑1 (MQ=255)

                                                              gccaccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGttt                 <  2:3914905/92‑1 (MQ=255)

                                                              accaccATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGttt                 <  2:2392404/93‑1 (MQ=255)

                                                              accaccATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGttt                 >  1:5194826/1‑93 (MQ=255)

                                                               ccaccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTg                >  2:1367537/1‑93 (MQ=255)

                                                               ccaccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTg                >  2:1793516/1‑93 (MQ=255)

                                                                 accaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGcc              >  2:2889935/1‑93 (MQ=255)

                                                                   caTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGt            >  1:1958441/1‑93 (MQ=255)

                                                                   caTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGt            >  2:32202/1‑93 (MQ=255)

                                                                    atatAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTa           >  1:2673101/1‑93 (MQ=255)

                                                                      ataAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAcc         <  2:4853678/93‑1 (MQ=255)

                                                                      ataAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAcc         <  2:2441847/93‑1 (MQ=255)

                                                                       taAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCa        <  2:41708/93‑1 (MQ=255)

                                                                       taAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCa        <  2:4343957/93‑1 (MQ=255)

                                                                       aaaGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCa        <  1:3409851/92‑1 (MQ=255)

                                                                          gCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGc     >  2:1764694/1‑93 (MQ=255)

                                                                          gCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGc     <  2:1243922/93‑1 (MQ=255)

                                                                             gagaATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAc  <  2:5281839/93‑1 (MQ=255)

                                                                             gagaATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAc  <  2:5599431/93‑1 (MQ=255)

TTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCCG‑GCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAC  >  NC_000913/4013092‑4013260

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