BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	526,302	(G)_6→7	coding (3818/4281 nt)	rhsD →	RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	526,296	1	.	G	99.2%	823.1 / ‑2.3	260	coding (3812/4281 nt)	rhsD	RhsD
Reads supporting (aligned to +/- strand): ref base . (2/0); new base G (139/119); total (141/119)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 305 total aligned reads displayed.

GAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGA  >  NC_000913/526204‑526383

 aaTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑g                                                                                         >  1:2118525/1‑93 (MQ=255)

  aTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑gg                                                                                        >  2:404769/1‑93 (MQ=255)

  aTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑gg                                                                                        <  1:7870596/93‑1 (MQ=255)

   tttATATCGGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑ggg                                                                                       <  2:2593955/93‑1 (MQ=255)

     tataTCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTTGGATGATGCCAGATCTGGAGCATGTAC‑‑ggggg                                                                                     >  1:5748484/1‑93 (MQ=255)

     tataTCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑ggggg                                                                                     <  2:6686628/93‑1 (MQ=255)

      atatCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑gggggg                                                                                    <  1:1827287/93‑1 (MQ=255)

        atCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑ggggggga                                                                                  >  2:150434/1‑91 (MQ=255)

        atCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑ggggggga                                                                                  <  2:5530427/93‑3 (MQ=255)

         tCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑gggggggag                                                                                 >  2:6599411/1‑90 (MQ=255)

         tCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑gggggggag                                                                                 <  1:2860770/93‑4 (MQ=255)

            gTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑Gttt                                                                              >  2:2900879/1‑93 (MQ=255)

            gTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑Gttt                                                                              >  1:3515911/1‑93 (MQ=255)

             tATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTtg                                                                             <  2:2674430/93‑1 (MQ=255)

               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGG‑GTTtgtg                                                                           >  2:6944148/1‑93 (MQ=255)

               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTtgtg                                                                           >  1:5040451/1‑93 (MQ=255)

               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTtgtg                                                                           >  2:3845203/1‑93 (MQ=255)

                ccTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTgtgg                                                                          <  1:3086108/93‑1 (MQ=255)

                 cTTTAAATCCACTACAACACATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTgtggt                                                                         >  1:5128641/1‑93 (MQ=255)

                  tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTgtggtg                                                                        >  1:2329012/1‑93 (MQ=255)

                  tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTgtggtg                                                                        >  1:4637970/1‑93 (MQ=255)

                    tAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGtt                                                                      <  1:1955372/93‑1 (MQ=255)

                    tAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGtt                                                                      >  2:1083498/1‑93 (MQ=255)

                      aaTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCt                                                                    <  2:6147544/93‑1 (MQ=255)

                       aTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCtt                                                                   >  1:241429/1‑93 (MQ=255)

                       aTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCtt                                                                   >  1:5813625/1‑93 (MQ=255)

                         ccACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTc                                                                 <  2:3580837/93‑1 (MQ=255)

                          cACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCa                                                                <  1:7105495/93‑1 (MQ=255)

                            cTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACg                                                              >  2:1646223/1‑93 (MQ=255)

                                aacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCAGGtata                                                          >  1:1584659/1‑93 (MQ=255)

                                 acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTataa                                                         <  2:6288103/93‑1 (MQ=255)

                                 acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTataa                                                         <  1:6497022/93‑1 (MQ=255)

                                  caaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTataat                                                        <  2:4354969/93‑1 (MQ=255)

                                  caaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTataat                                                        <  2:466908/93‑1 (MQ=255)

                                  caaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTataat                                                        <  2:6667749/93‑1 (MQ=255)

                                   aaaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTataata                                                       >  2:569608/1‑93 (MQ=255)

                                   aaaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTataata                                                       >  1:3375740/1‑93 (MQ=255)

                                    aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAg                                                      <  1:7237958/93‑1 (MQ=255)

                                    aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAg                                                      >  2:4027758/1‑93 (MQ=255)

                                    aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAg                                                      >  1:3902250/1‑93 (MQ=255)

                                     aTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAgg                                                     >  2:2844071/1‑93 (MQ=255)

                                      ttGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGa                                                    <  1:4101060/93‑1 (MQ=255)

                                      ttGACCCTATGGGATTATTGCAGACTTGGGATGAAACCAGATCTGGAGCATGTAC‑‑GGGGGGAGGTTTGTGATGTTATTTCACGTATACTAGGa                                                    >  1:2501207/1‑93 (MQ=255)

                                      ttGACCCTATGGGATTATTGCAGACGTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGa                                                    <  1:5345989/93‑1 (MQ=255)

                                          cccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCaa                                                >  1:3671328/1‑93 (MQ=255)

                                            cTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGt                                              <  2:7571753/93‑1 (MQ=255)

                                            cTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGt                                              >  1:2680493/1‑93 (MQ=255)

                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                             >  1:5606092/1‑93 (MQ=255)

                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                             <  1:4704098/93‑1 (MQ=255)

                                               tGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                           <  1:6303783/93‑1 (MQ=255)

                                                 ggATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAtt                                         <  2:3247130/93‑1 (MQ=255)

                                                 ggATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAtt                                         <  1:424712/93‑1 (MQ=255)

                                                  gATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAttt                                        >  1:6078938/1‑93 (MQ=255)

                                                    ttattGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTCTCACGTATAATAGGACCAAGTAAATTTgt                                      >  2:2246825/1‑92 (MQ=255)

                                                      attGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATa                                    <  2:2729249/93‑1 (MQ=255)

                                                       ttGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                   >  2:317077/1‑93 (MQ=255)

                                                       ttGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                   <  2:5199316/93‑1 (MQ=255)

                                                       ttGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                   >  1:1284585/1‑93 (MQ=255)

                                                        tGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGGGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGt                                  >  1:2196822/1‑93 (MQ=255)

                                                           aGACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACt                               >  1:6969587/1‑93 (MQ=255)

                                                            gACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTg                              >  1:2812065/1‑93 (MQ=255)

                                                            gACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTg                              >  2:3452026/1‑93 (MQ=255)

                                                            gACTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTg                              <  1:5453988/93‑1 (MQ=255)

                                                              cTTGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGGCCAAGTAAATTTGATAGTACTGCa                            <  1:4868300/93‑1 (MQ=255)

                                                               ttGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAg                           >  2:7168135/1‑93 (MQ=255)

                                                                tGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGa                          <  1:145292/93‑1 (MQ=255)

                                                                tGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGa                          <  1:3253993/93‑1 (MQ=255)

                                                                tGGGATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGa                          <  2:369765/93‑1 (MQ=255)

                                                                 gggATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAt                         >  2:6400742/1‑93 (MQ=255)

                                                                 gggATGATGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAt                         <  1:3708640/93‑1 (MQ=255)

                                                                   gatgatGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgc                       <  1:6845458/93‑1 (MQ=255)

                                                                   gatgatGCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgc                       <  2:4147036/93‑1 (MQ=255)

                                                                    atgatgCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgct                      <  2:7744588/93‑1 (MQ=255)

                                                                    atgatgCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgct                      <  2:780836/93‑1 (MQ=255)

                                                                     tgatgCCAGATCTGGAGCATGTACGGGGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgct                      <  2:2463006/93‑1 (MQ=255)

                                                                     tgatgCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctg                     >  2:2097617/1‑93 (MQ=255)

                                                                       atgCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCg                   >  1:5177526/1‑93 (MQ=255)

                                                                        tgCCAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGt                  >  2:4505194/1‑93 (MQ=255)

                                                                          ccAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                >  2:1751893/1‑93 (MQ=255)

                                                                          ccAGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                <  1:4565552/93‑1 (MQ=255)

                                                                            aGATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGa              >  2:1915812/1‑93 (MQ=255)

                                                                             gATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGAt             <  2:1781981/93‑1 (MQ=255)

                                                                             gATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGAt             <  2:5056354/93‑1 (MQ=255)

                                                                             gATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGAt             >  2:1411570/1‑93 (MQ=255)

                                                                             gATCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGAt             >  1:7722312/1‑93 (MQ=255)

                                                                              aTCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATg            >  1:4461225/1‑93 (MQ=255)

                                                                              aTCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATg            >  2:969775/1‑93 (MQ=255)

                                                                              aTCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATg            <  2:3142072/93‑1 (MQ=255)

                                                                              aTCTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATg            <  2:4471225/93‑1 (MQ=255)

                                                                                cTGGAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCt          >  1:6610121/1‑93 (MQ=255)

                                                                                  ggAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCttt        >  1:753147/1‑93 (MQ=255)

                                                                                  ggAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCttt        <  2:5358425/93‑1 (MQ=255)

                                                                                  ggAGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCtta        <  2:4794425/93‑2 (MQ=255)

                                                                                    aGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTga      <  1:7610038/93‑1 (MQ=255)

                                                                                    aGCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTga      >  1:1105177/1‑93 (MQ=255)

                                                                                     gCATGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaa     <  2:5258085/93‑1 (MQ=255)

                                                                                       aTGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaag   <  1:948484/93‑1 (MQ=255)

                                                                                       aTGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaag   <  2:7411709/93‑1 (MQ=255)

                                                                                        tGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga  <  1:7780737/93‑1 (MQ=255)

                                                                                        tGTACG‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga  <  2:3793032/93‑1 (MQ=255)

GAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGA  >  NC_000913/526204‑526383

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