BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,583,050	G→A	Q446* (CAA→TAA)	hsdR ←	endonuclease R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,583,050	0	G	A	99.7%	1046.6 / ‑6.4	284	Q446* (CAA→TAA)	hsdR	endonuclease R
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (148/135); minor base C (0/1); total (148/136)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 290 total aligned reads displayed.

TTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGA  >  NC_000913/4582960‑4583141

tttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACg                                                                                           >  1:1103714/1‑93 (MQ=255)

tttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACg                                                                                           <  2:167375/93‑1 (MQ=255)

 ttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACgg                                                                                          <  1:2771902/93‑1 (MQ=255)

 ttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACgg                                                                                          >  2:2788787/1‑93 (MQ=255)

   cAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGtt                                                                                        <  2:2689527/93‑1 (MQ=255)

     aCCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCg                                                                                      <  1:1629755/93‑1 (MQ=255)

      ccGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGc                                                                                     <  1:4198741/93‑1 (MQ=255)

      ccGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGc                                                                                     <  2:2615323/93‑1 (MQ=255)

       cgcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCt                                                                                    <  1:5188191/93‑1 (MQ=255)

        gcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTg                                                                                   <  2:5699321/93‑1 (MQ=255)

         cgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGc                                                                                  <  2:4882374/93‑1 (MQ=255)

         cgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGc                                                                                  <  2:5920227/93‑1 (MQ=255)

         cgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGc                                                                                  <  1:226653/93‑1 (MQ=255)

         cgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGc                                                                                  <  2:162619/93‑1 (MQ=255)

           cgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCg                                                                                >  1:4826980/1‑93 (MQ=255)

           cgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCg                                                                                <  2:1031886/93‑1 (MQ=255)

           cgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCg                                                                                <  1:5523471/93‑1 (MQ=255)

            ggACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGa                                                                               >  2:596184/1‑93 (MQ=255)

                 gCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAtt                                                                          >  1:3790181/1‑93 (MQ=255)

                 gCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAtt                                                                          <  1:6155125/93‑1 (MQ=255)

                  cATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttt                                                                         >  1:1236950/1‑93 (MQ=255)

                  cATCTTCCTGATAATAACGCAGGACCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttt                                                                         >  1:475061/1‑93 (MQ=255)

                  cATCTTCCTGATAATAACGCAGCCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttt                                                                         <  2:2183584/93‑1 (MQ=255)

                    tCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTcc                                                                       <  2:4075165/93‑1 (MQ=255)

                    tCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTcc                                                                       >  1:1489564/1‑93 (MQ=255)

                      ttCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGTTGCCGAGCATTTcacg                                                                     >  2:1203106/1‑90 (MQ=255)

                      ttCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcag                                                                     >  1:3271388/1‑93 (MQ=255)

                       tCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCccgc                                                                    >  2:5035214/1‑93 (MQ=255)

                       tCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagc                                                                    <  1:230266/93‑1 (MQ=255)

                       tCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagc                                                                    <  1:852146/93‑1 (MQ=255)

                           gATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                                                <  2:3165325/93‑1 (MQ=255)

                            ataataACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                                               >  2:890360/1‑93 (MQ=255)

                             taataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctct                                                              <  1:2438649/93‑1 (MQ=255)

                               ataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTc                                                            >  1:1405795/1‑93 (MQ=255)

                               ataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTc                                                            <  2:5847503/93‑1 (MQ=255)

                                   cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCggg                                                        <  1:3969855/93‑1 (MQ=255)

                                   cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCggg                                                        <  2:1342636/93‑1 (MQ=255)

                                      aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCgg                                                     <  2:4814208/93‑1 (MQ=255)

                                      aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCgg                                                     >  1:309457/1‑93 (MQ=255)

                                       ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGt                                                    <  1:3334877/93‑1 (MQ=255)

                                        gCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTg                                                   <  2:46086/93‑1 (MQ=255)

                                        gCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTg                                                   >  2:1481467/1‑93 (MQ=255)

                                          ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGcc                                                 >  1:1818187/1‑93 (MQ=255)

                                          ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGcc                                                 >  1:4879727/1‑93 (MQ=255)

                                             gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACt                                              <  2:3905700/93‑1 (MQ=255)

                                             gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACt                                              >  1:5891720/1‑93 (MQ=255)

                                                 gctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggg                                          <  1:5648078/93‑1 (MQ=255)

                                                 gctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggg                                          <  1:787880/93‑1 (MQ=255)

                                                 gctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggg                                          >  1:5445982/1‑93 (MQ=255)

                                                  ctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGt                                         >  2:248782/1‑93 (MQ=255)

                                                   tcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTa                                        >  1:3224740/1‑93 (MQ=255)

                                                   tcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTa                                        >  2:2036271/1‑93 (MQ=255)

                                                    cATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTaa                                       <  2:1760101/93‑1 (MQ=255)

                                                       gCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGc                                    >  2:685157/1‑93 (MQ=255)

                                                        ccAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCCGGCGGGGCCACTCGGGGTAGGCt                                   >  2:4786977/1‑93 (MQ=255)

                                                         cAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                                  >  2:4272047/1‑93 (MQ=255)

                                                         cAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                                  >  1:5839901/1‑93 (MQ=255)

                                                         cAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                                  >  2:3203942/1‑93 (MQ=255)

                                                         cAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGCAAGGCtt                                  >  2:2374451/1‑93 (MQ=255)

                                                          aGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCttt                                 <  1:4319323/93‑1 (MQ=255)

                                                          aGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCttt                                 <  2:5369379/93‑1 (MQ=255)

                                                          aGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCttt                                 <  1:709554/93‑1 (MQ=255)

                                                           gggTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTc                                <  1:5770784/93‑1 (MQ=255)

                                                             gTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGa                              <  2:4690628/93‑1 (MQ=255)

                                                             gTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGa                              >  2:437661/1‑93 (MQ=255)

                                                             gTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGa                              >  1:5980/1‑93 (MQ=255)

                                                              ttATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGAc                             >  1:1528972/1‑93 (MQ=255)

                                                               tATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACa                            >  2:623759/1‑93 (MQ=255)

                                                               tATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACa                            >  1:1983498/1‑93 (MQ=255)

                                                                aTCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACat                           <  1:2099639/93‑1 (MQ=255)

                                                                 tCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACata                          >  2:4498283/1‑93 (MQ=255)

                                                                  cGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACatat                         <  1:1749917/93‑1 (MQ=255)

                                                                  cGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACatat                         <  2:4191405/93‑1 (MQ=255)

                                                                   ggCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATAtt                        <  2:754186/93‑1 (MQ=255)

                                                                   ggCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATAtt                        >  2:3416289/1‑93 (MQ=255)

                                                                   ggCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGGGCCACTCGGGTAAGGCTTTCGACATAtt                        >  1:3535140/1‑93 (MQ=255)

                                                                      aaaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCg                     >  2:30097/1‑93 (MQ=255)

                                                                      aaaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCg                     <  1:3025668/93‑1 (MQ=255)

                                                                      aaaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCg                     >  2:1702882/1‑93 (MQ=255)

                                                                      aaaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCg                     <  1:2560069/93‑1 (MQ=255)

                                                                       aaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCGACTCGGGGAAGGCTTTCGACATATTGCgg                    >  2:335540/1‑93 (MQ=255)

                                                                        aCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCggg                   >  2:2868074/1‑93 (MQ=255)

                                                                          cACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTa                 >  2:4956528/1‑93 (MQ=255)

                                                                          cACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTa                 >  1:5622848/1‑93 (MQ=255)

                                                                          cACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTa                 <  1:6047640/93‑1 (MQ=255)

                                                                            cTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATc               <  1:2183253/93‑1 (MQ=255)

                                                                            cTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATc               >  2:5141640/1‑93 (MQ=255)

                                                                              gATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcac             <  1:20896/93‑1 (MQ=255)

                                                                              gATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcac             >  2:1958284/1‑93 (MQ=255)

                                                                              gATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcac             >  1:2953030/1‑93 (MQ=255)

                                                                               aTTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacg            >  2:1593337/1‑93 (MQ=255)

                                                                                ttCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgc           <  1:4685822/93‑1 (MQ=255)

                                                                                ttCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgc           >  2:1521776/1‑93 (MQ=255)

                                                                                ttCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgc           >  1:4095095/1‑93 (MQ=255)

                                                                                 tCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgca          >  1:4495996/1‑93 (MQ=255)

                                                                                 tCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgca          >  1:5250873/1‑93 (MQ=255)

                                                                                  cTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgcac         <  2:1777575/93‑1 (MQ=255)

                                                                                   tGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgcacg        >  2:734039/1‑93 (MQ=255)

                                                                                     tttttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTc      <  2:5571497/93‑1 (MQ=255)

                                                                                     tttttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTc      >  2:5692843/1‑93 (MQ=255)

TTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGA  >  NC_000913/4582960‑4583141

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