BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	332,530	T→C	P312P (CCT→CCC)	yahA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	332,530	0	T	C	99.8%	1378.6 / ‑5.6	427	P312P (CCT→CCC)	yahA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (189/237); total (190/237)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 441 total aligned reads displayed.

TAAGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCTGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTATTT  >  NC_000913/332440‑332619

taAGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCgg                                                                                         >  2:7970506/1‑93 (MQ=255)

 aaGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCggc                                                                                        >  1:5546092/1‑92 (MQ=255)

 aaGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGt                                                                                        <  1:1610389/93‑1 (MQ=255)

 aaGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAACCCGGt                                                                                        <  2:6001601/93‑1 (MQ=255)

     tCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGa                                                                                    <  2:2849018/93‑1 (MQ=255)

       gATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGt                                                                                  <  1:1529852/93‑1 (MQ=255)

       gATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGt                                                                                  >  1:2920118/1‑93 (MQ=255)

        aTAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTACGCCCGGTCTGAGTa                                                                                 >  1:8019631/1‑93 (MQ=255)

        aTAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTa                                                                                 >  2:8594543/1‑93 (MQ=255)

        aTAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGCAAGCCCGGTCTGAGTa                                                                                 <  2:2239971/93‑1 (MQ=255)

         tAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTAt                                                                                >  1:9508060/1‑93 (MQ=255)

         tAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTAt                                                                                <  2:3841331/93‑1 (MQ=255)

          aaGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATc                                                                               >  2:4288717/1‑93 (MQ=255)

           aGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCg                                                                              <  1:1698050/93‑1 (MQ=255)

           aGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCg                                                                              >  1:6101832/1‑93 (MQ=255)

            gTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAATCCCGGTCTGAGTATCGt                                                                             <  1:366027/93‑1 (MQ=255)

            gTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGt                                                                             >  2:7025116/1‑93 (MQ=255)

             tCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTg                                                                            <  2:9254146/93‑1 (MQ=255)

              cATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTgg                                                                           <  1:7496419/93‑1 (MQ=255)

              cATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTgg                                                                           <  1:3829979/93‑1 (MQ=255)

               aTTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGc                                                                          >  2:2697649/1‑93 (MQ=255)

               aTTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGc                                                                          >  2:3006201/1‑93 (MQ=255)

               aTTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGc                                                                          >  2:9613635/1‑93 (MQ=255)

               aTTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGc                                                                          <  2:7152767/93‑1 (MQ=255)

                tttGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCg                                                                         >  2:4782047/1‑93 (MQ=255)

                tttGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCg                                                                         >  1:3230373/1‑93 (MQ=255)

                 ttGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCgg                                                                        >  1:9181553/1‑93 (MQ=255)

                   gtgCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGaa                                                                      <  2:9681601/93‑1 (MQ=255)

                     gCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAgg                                                                    >  1:5163515/1‑93 (MQ=255)

                     gCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAgg                                                                    <  2:3138589/93‑1 (MQ=255)

                      cAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAggg                                                                   >  1:7663035/1‑93 (MQ=255)

                        aaTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGt                                                                 <  2:2043097/93‑1 (MQ=255)

                         aTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCTGGTCTGAGTATCGTGGCGGAAGGCGTa                                                                >  2:4159976/1‑93 (MQ=255)

                         aTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTa                                                                <  1:7810481/93‑1 (MQ=255)

                         aTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTa                                                                <  2:2275607/93‑1 (MQ=255)

                           ggCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGa                                                              <  1:9658542/93‑1 (MQ=255)

                           ggCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGa                                                              >  1:5012804/1‑93 (MQ=255)

                             cGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGaaa                                                            <  2:6300392/93‑1 (MQ=255)

                             cGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGaaa                                                            <  1:6010444/93‑1 (MQ=255)

                              gAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAc                                                           <  2:4368212/93‑1 (MQ=255)

                                gtgtCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAccc                                                         <  1:1992309/93‑1 (MQ=255)

                                gtgtCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAccc                                                         >  1:3487218/1‑93 (MQ=255)

                                gtgtCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAccc                                                         <  1:8281666/93‑1 (MQ=255)

                                  gtCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAg                                                       <  1:4612113/93‑1 (MQ=255)

                                   tCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAgg                                                      >  1:4444965/1‑93 (MQ=255)

                                    agacgaAAGCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGa                                                     <  2:5536298/92‑1 (MQ=255)

                                     gacgaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAg                                                    <  1:9907507/93‑1 (MQ=255)

                                      acgaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGc                                                   >  2:5684539/1‑93 (MQ=255)

                                       cgaAATCCTCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCa                                                  <  1:2705677/93‑1 (MQ=255)

                                        gaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAg                                                 >  2:1065156/1‑93 (MQ=255)

                                        gaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAg                                                 >  2:1655372/1‑93 (MQ=255)

                                         aaaTCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGGAGACACCCAGGAGCAgg                                                >  2:89554/1‑93 (MQ=255)

                                            tctcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCgg                                             <  1:4127231/93‑1 (MQ=255)

                                            tctcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCgg                                             >  2:4973292/1‑93 (MQ=255)

                                            tctcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCgg                                             >  1:6381710/1‑93 (MQ=255)

                                              tcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGAt                                           >  1:4735457/1‑93 (MQ=255)

                                              tcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGAt                                           <  2:8182316/93‑1 (MQ=255)

                                                  gtaatatTGTGGAAAATATTGTCGAACTAGCGCGTAAGCGGGGTCTGAGTATGGTGGGGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTaa                                       <  2:7425600/90‑1 (MQ=255)

                                                  gTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGTGTAGCGCGGCGGAAGGGGGAGAAACCCAGGAGCAGGCGGATTTaa                                       >  1:8976695/1‑93 (MQ=255)

                                                   tCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAAt                                      >  2:7736170/1‑93 (MQ=255)

                                                   tCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAAt                                      >  1:657471/1‑93 (MQ=255)

                                                    cATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATg                                     >  2:1253782/1‑93 (MQ=255)

                                                     atagTGTGGACAATAGTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGa                                    <  1:4798290/93‑1 (MQ=255)

                                                      tatTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGAt                                   >  1:2560290/1‑93 (MQ=255)

                                                      tatTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGAt                                   <  2:6333879/93‑1 (MQ=255)

                                                        ttGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTAGCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCg                                 <  2:9506514/93‑1 (MQ=255)

                                                           tgGACAATATTTTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTa                              <  1:6454081/93‑1 (MQ=255)

                                                            ggACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTaa                             <  2:6556501/93‑1 (MQ=255)

                                                            ggACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTaa                             >  1:5405750/1‑93 (MQ=255)

                                                            ggACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTaa                             <  1:7412682/93‑1 (MQ=255)

                                                             gACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGTTTTAATGATCGGTaaa                            <  2:3716879/93‑1 (MQ=255)

                                                             gACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTaag                            >  1:8762245/1‑92 (MQ=255)

                                                               cAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAgg                          >  1:6849511/1‑93 (MQ=255)

                                                                 atatTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCg                        >  2:5837505/1‑93 (MQ=255)

                                                                 atatTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCg                        >  1:9088292/1‑93 (MQ=255)

                                                                 atatTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCg                        >  1:2279199/1‑93 (MQ=255)

                                                                  tatTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGt                       >  2:4452248/1‑93 (MQ=255)

                                                                    ttGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGGGCAGGCGGATTTAATGATCGGTAAAGGCGTTc                     <  2:8416367/93‑1 (MQ=255)

                                                                    ttGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTc                     <  1:6218278/93‑1 (MQ=255)

                                                                         gAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACttt                >  1:9332409/1‑93 (MQ=255)

                                                                          caCTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGGTTTAATGATTGGTAAATCCTTTCAGtttt               >  1:1286260/2‑93 (MQ=255)

                                                                          aaCTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACtttt               >  2:335801/1‑93 (MQ=255)

                                                                          aaCTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGCAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCAAtttt               >  2:3579974/1‑93 (MQ=255)

                                                                           aCTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACttttt              <  2:8060023/93‑1 (MQ=255)

                                                                           aCTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACttttt              >  1:8412710/1‑93 (MQ=255)

                                                                            cTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTg             <  2:7231517/93‑1 (MQ=255)

                                                                             tAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGGAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGc            >  2:5300824/1‑93 (MQ=255)

                                                                              aGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCa           <  1:10001064/93‑1 (MQ=255)

                                                                              aGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCa           <  1:5275546/93‑1 (MQ=255)

                                                                               gcgcGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAg          <  2:6693703/93‑1 (MQ=255)

                                                                                cgcgTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAgg         >  2:4722762/1‑93 (MQ=255)

                                                                                 gcgTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGAAAAGGCGTTCACTTTTTGCAggg        >  2:6932709/1‑93 (MQ=255)

                                                                                   gTAAGCCCGGCCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGATCAGGCGGGTTTAATGATCGGTAAAGGCGTTTACTTTTTGCTGGGCt      >  2:6791522/1‑93 (MQ=255)

                                                                                     aaGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAt    <  2:9950022/93‑1 (MQ=255)

                                                                                      aGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAtt   <  2:2460047/93‑1 (MQ=255)

                                                                                      aGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAtt   >  1:7143284/1‑93 (MQ=255)

                                                                                      aGCCCGGTCTGAGTATCGTGGCGGAAGGGGGAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAtt   >  1:5889152/1‑93 (MQ=255)

                                                                                       gCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAttt  <  2:8756309/93‑1 (MQ=255)

                                                                                       gCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAttt  <  1:315384/93‑1 (MQ=255)

                                                                                       gCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAttt  <  1:1160870/93‑1 (MQ=255)

TAAGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCTGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTATTT  >  NC_000913/332440‑332619

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