BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	683,225	A→G	H137H (CAT→CAC)	ybeK ←	putative tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	683,225	0	A	G	84.4%	602.2 / 105.5	263	H137H (CAT→CAC)	ybeK	putative tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (21/20); new base G (113/109); total (134/129)
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.87e-01

Only 100 of 264 total aligned reads displayed.

TAAATATTAAATTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTCGCCATC  >  NC_000913/683133‑683315

    tatTAAATTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGCGCAATTTTGCTATGCAGTTCCGGGTGGc                                                                                        <  2:6290640/93‑1 (MQ=255)

     atTAAATTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCt                                                                                       >  2:1804632/1‑93 (MQ=255)

      ttAAATTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTa                                                                                      >  2:4969985/1‑93 (MQ=255)

      ttAAATTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTa                                                                                      <  2:5752441/93‑1 (MQ=255)

       tAAATTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTAt                                                                                     >  2:5075043/1‑93 (MQ=255)

         aaTTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTg                                                                                   >  1:4910026/1‑93 (MQ=255)

         aaTTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTg                                                                                   >  2:5255702/1‑93 (MQ=255)

          aTTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGa                                                                                  <  1:5408643/93‑1 (MQ=255)

          aTTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGa                                                                                  <  2:1323619/93‑1 (MQ=255)

          aTTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGa                                                                                  <  2:6868124/93‑1 (MQ=255)

            tCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGagc                                                                                <  1:1584333/93‑1 (MQ=255)

               gCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGagcagc                                                                             <  1:6428222/93‑1 (MQ=255)

                 cGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCaa                                                                           >  2:3076047/1‑93 (MQ=255)

                 cGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCaa                                                                           >  2:7140020/1‑93 (MQ=255)

                  gCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAg                                                                          <  2:2725356/93‑1 (MQ=255)

                   cAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAgg                                                                         <  1:875932/93‑1 (MQ=255)

                   cAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAgg                                                                         <  2:1909188/93‑1 (MQ=255)

                    aGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGc                                                                        >  1:1109652/1‑93 (MQ=255)

                     ggCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCa                                                                       >  1:3233888/1‑93 (MQ=255)

                     ggCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCa                                                                       >  1:2886306/1‑93 (MQ=255)

                      gCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCaa                                                                      >  2:4282587/1‑93 (MQ=255)

                      gCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCaa                                                                      <  2:5203250/93‑1 (MQ=255)

                       cGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAAc                                                                     <  2:4452285/93‑1 (MQ=255)

                       cGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAAc                                                                     <  1:1268708/93‑1 (MQ=255)

                         tCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGt                                                                   <  1:3326231/93‑1 (MQ=255)

                         tCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGt                                                                   >  1:728492/1‑93 (MQ=255)

                         tCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGt                                                                   >  2:7340577/1‑93 (MQ=255)

                          ccAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGtt                                                                  >  1:5367235/1‑93 (MQ=255)

                          ccAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGtt                                                                  >  1:1877547/1‑93 (MQ=255)

                           cAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTa                                                                 <  2:1989207/93‑1 (MQ=255)

                            aGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAg                                                                >  1:2360210/1‑93 (MQ=255)

                             gTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGt                                                               >  1:6538734/1‑93 (MQ=255)

                             gTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGt                                                               >  2:6825490/1‑93 (MQ=255)

                             gTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGt                                                               >  2:3321812/1‑93 (MQ=255)

                              ttACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGtt                                                              <  1:7344892/93‑1 (MQ=255)

                                 cccAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGc                                                           >  2:5624/1‑93 (MQ=255)

                                 cccAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGc                                                           >  2:4739349/1‑93 (MQ=255)

                                 cccAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGc                                                           >  2:4087615/1‑93 (MQ=255)

                                  ccAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCg                                                          >  1:3155303/1‑93 (MQ=255)

                                  ccAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCg                                                          <  2:3261071/93‑1 (MQ=255)

                                    aGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGt                                                        >  2:3953952/1‑93 (MQ=255)

                                     gCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTc                                                       >  1:268286/1‑93 (MQ=255)

                                      ccccATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTcc                                                      <  2:3989755/93‑1 (MQ=255)

                                      ccccATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTcc                                                      <  2:1544912/93‑1 (MQ=255)

                                      ccccATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTcc                                                      <  1:4020853/93‑1 (MQ=255)

                                       cccATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCg                                                     >  1:100698/1‑93 (MQ=255)

                                          aTAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTa                                                  >  1:252971/1‑93 (MQ=255)

                                              cgcCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGaca                                              <  2:6146676/93‑1 (MQ=255)

                                              cgcCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGaca                                              <  2:4224953/93‑1 (MQ=255)

                                               gcCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGacac                                             >  2:2451436/1‑93 (MQ=255)

                                                ccACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGAcaca                                            >  2:5440735/1‑93 (MQ=255)

                                                ccACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCAGGTAGAcaca                                            >  1:2827791/1‑93 (MQ=255)

                                                 cACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACaa                                           >  1:5162132/1‑93 (MQ=255)

                                                 cACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACaa                                           <  1:5071931/93‑1 (MQ=255)

                                                 cACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACaa                                           <  1:4361394/93‑1 (MQ=255)

                                                 cACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACaa                                           >  1:3710323/1‑93 (MQ=255)

                                                  aCCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAAt                                          <  1:5903402/93‑1 (MQ=255)

                                                  aCCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAAt                                          <  1:2482288/93‑1 (MQ=255)

                                                  aCCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAAt                                          <  1:2149806/93‑1 (MQ=255)

                                                  aCCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAAt                                          >  1:6139983/1‑93 (MQ=255)

                                                   cccATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATg                                         <  2:437230/93‑1 (MQ=255)

                                                   cccATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATg                                         <  2:6698983/93‑1 (MQ=255)

                                                   cccATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATg                                         <  2:773871/93‑1 (MQ=255)

                                                    ccATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATgg                                        >  1:5567454/1‑93 (MQ=255)

                                                     cATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGt                                       >  2:3748663/1‑93 (MQ=255)

                                                     cATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGt                                       >  2:4600854/1‑93 (MQ=255)

                                                      aTAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTg                                      >  1:7073551/1‑93 (MQ=255)

                                                       tAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGa                                     <  2:221375/93‑1 (MQ=255)

                                                       tAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGa                                     >  1:3437170/1‑93 (MQ=255)

                                                       tAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGa                                     >  1:6641485/1‑93 (MQ=255)

                                                         aTCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACa                                   >  2:5359056/1‑93 (MQ=255)

                                                         aTCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACa                                   <  1:1951447/93‑1 (MQ=255)

                                                          tCACGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAg                                  >  1:7260005/1‑93 (MQ=255)

                                                            aCGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGt                                >  2:1963717/1‑93 (MQ=255)

                                                             cGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGtt                               >  1:1835413/1‑93 (MQ=255)

                                                             cGATACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGtt                               >  1:7482135/1‑93 (MQ=255)

                                                               aTACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTcc                             >  1:4275436/1‑93 (MQ=255)

                                                                tACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCg                            <  1:1423901/93‑1 (MQ=255)

                                                                tACGGGCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCg                            >  1:6800555/1‑93 (MQ=255)

                                                                   gggCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCAc                         <  2:3572696/93‑1 (MQ=255)

                                                                   gggCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCAc                         >  1:4698902/1‑93 (MQ=255)

                                                                   gggCAATTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCAc                         <  2:1877850/93‑1 (MQ=255)

                                                                     gCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACtt                       <  2:5471724/93‑1 (MQ=255)

                                                                       aaTTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTc                     >  1:2029468/1‑93 (MQ=255)

                                                                        aTTTTGCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCa                    >  2:1170825/1‑93 (MQ=255)

                                                                             gCTATGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAg               >  1:3867136/1‑93 (MQ=255)

                                                                                 tGCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGtt           >  1:6716000/1‑93 (MQ=255)

                                                                                  gCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGttt          >  2:2340275/1‑93 (MQ=255)

                                                                                  gCAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGttt          >  2:3108918/1‑93 (MQ=255)

                                                                                   cAGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGtttt         <  1:944119/93‑1 (MQ=255)

                                                                                    aGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTc        <  2:6458923/93‑1 (MQ=255)

                                                                                    aGTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTc        <  1:4425976/93‑1 (MQ=255)

                                                                                     gTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTCg       >  1:4345105/1‑93 (MQ=255)

                                                                                     gTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTCg       >  1:4812432/1‑93 (MQ=255)

                                                                                     gTTCCGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTCg       >  2:7298706/1‑93 (MQ=255)

                                                                                        ccGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTCGCCa    <  1:4123476/93‑1 (MQ=255)

                                                                                        ccGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTCGCCa    <  2:1277212/93‑1 (MQ=255)

                                                                                        ccGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTCGCCa    <  2:2821876/93‑1 (MQ=255)

                                                                                         cGGGTGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTCGCCAt   >  2:169631/1‑93 (MQ=255)

                                                                                          ggATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTCGCCATc  >  1:789013/1‑93 (MQ=255)

TAAATATTAAATTCAGCCGCAGGCGTCCAGTTACCCAGCCCCATAGCGCCACCCATAATCACGATACGGGCAATTTTGCTATGCAGTTCCGGATGGCTATTGAGCAGCAAGGCAACGTTAGTTTGCGGTCCGGTAGACACAATGGTGACAGGTTCCGCACTTTCACGCAGCGTTTTCGCCATC  >  NC_000913/683133‑683315

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