BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,051,163	Δ1 bp	coding (94/231 nt)	sfa →	Sfa

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,051,162	0	A	.	99.6%	1278.4 / ‑2.4	275	coding (93/231 nt)	sfa	Sfa
Reads supporting (aligned to +/- strand): ref base A (1/0); new base . (143/131); total (144/131)
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.72e-01

Only 100 of 291 total aligned reads displayed.

ATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAC  >  NC_000913/1051070‑1051250

aTGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTa                                                                                          >  1:8353154/1‑93 (MQ=255)

aTGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTa                                                                                          <  2:1476223/93‑1 (MQ=255)

aTGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTa                                                                                          <  1:7008966/93‑1 (MQ=255)

    gccgcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTagcga                                                                                      <  2:7508153/93‑5 (MQ=255)

       gcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATaa                                                                                  >  2:5790722/1‑93 (MQ=255)

       gcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATaa                                                                                  >  1:7564868/1‑93 (MQ=255)

       gcTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATaa                                                                                  >  2:5158194/1‑93 (MQ=255)

        cTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAAc                                                                                 >  1:2091794/1‑93 (MQ=255)

         tGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAAcc                                                                                <  2:7033078/93‑1 (MQ=255)

         tGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAAcc                                                                                <  2:3122895/93‑1 (MQ=255)

          ggATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCt                                                                               <  1:5301768/93‑1 (MQ=255)

          ggATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTAAATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCt                                                                               <  1:1575684/93‑1 (MQ=255)

            aTTTCACAAAATAATATCAGGCTCCCTCGTGGAGCGTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGAGATAAACTGt                                                                             >  1:4444242/1‑93 (MQ=255)

            aTTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGt                                                                             <  1:1157109/93‑1 (MQ=255)

             tttCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGtt                                                                            <  2:754894/93‑1 (MQ=255)

             tttCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGtt                                                                            <  2:5479861/93‑1 (MQ=255)

              ttCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTg                                                                           <  1:4857410/93‑1 (MQ=255)

               tCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGCATGTATAGT‑AGCGATAACCTGTTGa                                                                          >  2:7635938/1‑93 (MQ=255)

                cacaAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGat                                                                         >  1:7195588/1‑93 (MQ=255)

                   aaaaTAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGattat                                                                      >  1:4074823/1‑93 (MQ=255)

                   aaaaTAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGattat                                                                      >  1:548676/1‑93 (MQ=255)

                    aaaTAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGattatt                                                                     <  1:6116943/93‑1 (MQ=255)

                    aaaTAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGattatt                                                                     <  2:6312622/93‑1 (MQ=255)

                    aaaTAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGattatt                                                                     <  2:1045350/93‑1 (MQ=255)

                     aataatATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTg                                                                    <  1:6502861/93‑1 (MQ=255)

                     aataatATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGAATGTATAGT‑AGCGATAACCTGTTGATTATTg                                                                    >  1:1081956/1‑93 (MQ=255)

                       taataTCAGGCTCCCTCGAGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGat                                                                  >  2:6161680/1‑92 (MQ=255)

                          tatCAGTCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCt                                                               <  2:2217899/93‑1 (MQ=255)

                           atCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCtt                                                              >  2:3318707/1‑93 (MQ=255)

                            tCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCttt                                                             <  1:4517866/93‑1 (MQ=255)

                            tCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCttt                                                             <  1:3743200/93‑1 (MQ=255)

                             cAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTc                                                            <  2:292742/93‑1 (MQ=255)

                             cAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGGATAGT‑AGCGCTAACCTGTTGATTATTGTATGTTTc                                                            >  1:1456070/1‑93 (MQ=255)

                             cAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTGTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTc                                                            <  2:2424102/93‑1 (MQ=255)

                              aGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCg                                                           <  2:4089636/93‑1 (MQ=255)

                               ggCTCCCTCGTGGAGCCTTTTTTCTATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCgg                                                          <  2:6283997/93‑1 (MQ=255)

                               ggCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCgg                                                          <  2:3955567/93‑1 (MQ=255)

                                gCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCggg                                                         <  1:3546817/93‑1 (MQ=255)

                                  tCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCgggca                                                       >  2:7462402/1‑91 (MQ=255)

                                  tCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGga                                                       <  1:6173538/93‑1 (MQ=255)

                                  tCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGga                                                       <  2:8112733/93‑1 (MQ=255)

                                  tCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGga                                                       <  2:3382695/93‑1 (MQ=255)

                                  tCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGga                                                       <  2:8289333/93‑1 (MQ=255)

                                   cccTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGgag                                                      >  1:2390986/1‑93 (MQ=255)

                                   cccTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGgag                                                      >  2:135506/1‑93 (MQ=255)

                                      tCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATg                                                   <  1:5186065/93‑1 (MQ=255)

                                        gTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGc                                                 >  1:5809657/1‑93 (MQ=255)

                                         tGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCt                                                >  1:1013927/1‑93 (MQ=255)

                                          ggAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCtt                                               >  1:6422174/1‑93 (MQ=255)

                                           gAGCGTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAGTCTTTCGGGGAGATGGCTta                                              <  2:5327104/93‑1 (MQ=255)

                                            aGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTtat                                             >  2:4640482/1‑93 (MQ=255)

                                            aGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTtat                                             <  1:7392407/93‑1 (MQ=255)

                                            aGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTtat                                             >  1:3033555/1‑93 (MQ=255)

                                              ccTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATaa                                           >  2:1868367/1‑93 (MQ=255)

                                              ccTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATaa                                           >  2:2010464/1‑93 (MQ=255)

                                              ccTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATaa                                           >  1:3895270/1‑93 (MQ=255)

                                               cTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAAc                                          <  1:1912551/93‑1 (MQ=255)

                                                ttttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACa                                         <  2:542599/93‑1 (MQ=255)

                                                ttttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACa                                         >  2:5011100/1‑93 (MQ=255)

                                                tttttaATACCTGCCTTATTTTTCTTCAACGCTGTATGGATAGT‑AGCGATAACCTGTTGATTATGGAATCTTACGGGGAGATGGCGAATAACa                                         >  2:6688337/1‑93 (MQ=255)

                                                  ttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAtt                                       >  1:8645180/1‑93 (MQ=255)

                                                  ttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAtt                                       <  2:4198308/93‑1 (MQ=255)

                                                  ttttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACAtt                                       <  1:165191/93‑1 (MQ=255)

                                                    ttATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGTTTATAACCTTTc                                     <  2:5599623/93‑1 (MQ=255)

                                                     tataTCTGCCTTGTTTTTCCTCAAAGCTGAATGTATCGGAAGCGATAACCTGTTGCTTATTGTATCTTGCCGCGAGGTGGGTTATAAAATTTc                                     >  2:4918187/1‑93 (MQ=255)

                                                       tatCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTa                                  >  2:266137/1‑93 (MQ=255)

                                                        atCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTAc                                 >  1:4417526/1‑93 (MQ=255)

                                                         tCTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTAcc                                >  1:338174/1‑93 (MQ=255)

                                                          cTGCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCt                               <  1:8072585/93‑1 (MQ=255)

                                                            gCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGa                             >  1:1040575/1‑93 (MQ=255)

                                                            gCCTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGa                             >  2:6466742/1‑93 (MQ=255)

                                                             ccTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAc                            >  2:8769554/1‑93 (MQ=255)

                                                             ccTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAc                            >  2:3674226/1‑93 (MQ=255)

                                                             ccTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAc                            <  1:800306/93‑1 (MQ=255)

                                                              cTTATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGAcc                           <  2:209924/93‑1 (MQ=255)

                                                               ttATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCa                          >  2:3558037/1‑93 (MQ=255)

                                                                tATTTTTCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAg                         >  2:1610751/1‑93 (MQ=255)

                                                                    tttCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTa                     >  2:3373129/1‑93 (MQ=255)

                                                                    tttCTTCAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTa                     <  2:4804440/93‑1 (MQ=255)

                                                                     ttcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAc                    >  2:8571588/1‑93 (MQ=255)

                                                                     ttcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAc                    >  1:3310448/1‑93 (MQ=255)

                                                                     ttcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAc                    >  1:6856889/1‑93 (MQ=255)

                                                                      tcttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAcc                   <  1:6360645/93‑1 (MQ=255)

                                                                       cttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCg                  <  2:2540116/93‑1 (MQ=255)

                                                                       cttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCg                  >  1:1266041/1‑93 (MQ=255)

                                                                        ttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCgg                 >  2:8384484/1‑93 (MQ=255)

                                                                        ttcAACGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCgg                 <  1:4537019/93‑1 (MQ=255)

                                                                          taaCGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGa               <  1:6817923/92‑1 (MQ=255)

                                                                           aaCGCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGaa              >  1:6635209/1‑93 (MQ=255)

                                                                              gCTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCa           <  2:6072429/93‑1 (MQ=255)

                                                                               cTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCaa          <  1:4218858/93‑1 (MQ=255)

                                                                               cTGTATGTATAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCaa          >  2:4852915/1‑93 (MQ=255)

                                                                                tgtatgtaTAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac         >  2:5952078/1‑93 (MQ=255)

                                                                                tgtatgtaTAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac         >  1:8444175/1‑93 (MQ=255)

                                                                                 gtatgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAaca        >  1:6666232/1‑93 (MQ=255)

                                                                                 gtatgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAaca        <  1:300037/93‑1 (MQ=255)

                                                                                  tatgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAacac       <  1:5600322/93‑1 (MQ=255)

                                                                                    tgtatAGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCt     >  2:6558619/1‑93 (MQ=255)

                                                                                       ataGTAA‑CGATAACCTGGTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCATTGTAACGTGAACAAAGACCTTaa  >  2:4951404/1‑92 (MQ=255)

                                                                                       ataGT‑AGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAc  >  1:3228557/1‑93 (MQ=255)

ATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAC  >  NC_000913/1051070‑1051250

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