BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,715,953	T→C	53.9%	K424E (AAA→GAA)	tyrS ←	tyrosyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,715,953	0	T	C	53.9%	39.1 / 195.6	155	K424E (AAA→GAA)	tyrS	tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base T (40/31); new base C (47/36); total (87/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

Only 100 of 158 total aligned reads displayed.

CGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTG  >  NC_000913/1715856‑1716052

cGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTcc                                                                                                   >  2:1002104/1‑100 (MQ=255)

cGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTACCAGAATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTcc                                                                                                   <  1:460976/100‑1 (MQ=255)

  tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAg                                                                                                 <  1:86951/100‑1 (MQ=255)

  tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTATCAATCCCCCTTTCCACATAATGCATTATTTCCAg                                                                                                 >  2:512706/1‑100 (MQ=255)

    agaCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCa                                                                                               >  2:3592286/1‑100 (MQ=255)

     gaCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCaa                                                                                              >  2:3887495/1‑100 (MQ=255)

           aTAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACACCCCCTTTACACTTAATGCATTATTTCCAGCAAATCAGa                                                                                        >  2:2759406/1‑100 (MQ=255)

            tAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGAc                                                                                       <  2:3123939/100‑1 (MQ=255)

             aGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACa                                                                                      <  2:3399587/100‑1 (MQ=255)

              gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAg                                                                                     <  1:135313/100‑1 (MQ=255)

              gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAg                                                                                     >  1:2181532/1‑100 (MQ=255)

               cGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGt                                                                                    <  1:253542/100‑1 (MQ=255)

               cGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAg                                                                                     >  2:888172/1‑99 (MQ=255)

                 agaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTaa                                                                                  <  1:3663588/100‑1 (MQ=255)

                 agaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTaa                                                                                  <  2:1869928/100‑1 (MQ=255)

                  gaATATTTTTCATCATGTTCCCCGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAAt                                                                                 >  2:3246886/1‑100 (MQ=255)

                       tttttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCttt                                                                            <  1:1061222/100‑1 (MQ=255)

                       tttttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCttt                                                                            <  2:3197183/100‑1 (MQ=255)

                        ttttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCtttt                                                                           >  1:4165882/1‑100 (MQ=255)

                         tttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCttttt                                                                          <  1:1291368/100‑1 (MQ=255)

                          ttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTa                                                                         <  2:6626/100‑1 (MQ=255)

                          ttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTATCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTa                                                                         >  1:1579973/1‑100 (MQ=255)

                            catcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTAcc                                                                       >  1:3489267/1‑100 (MQ=255)

                            catcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTAcc                                                                       >  1:1706579/1‑100 (MQ=255)

                              tcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACcgc                                                                     <  1:1002101/100‑1 (MQ=255)

                               catGTTCCCTGTATAAAAACCAGGGGGAGTGAATTCTCACTCCCCCTTTCCACTTAATGCATTATTCTCAGCTAATCAGACAGTAATTCTTTTTACcg                                                                      >  1:871123/1‑98 (MQ=255)

                                   ttCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGc                                                                <  2:198081/100‑1 (MQ=255)

                                   ttCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGc                                                                <  1:1040303/100‑1 (MQ=255)

                                    tCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCa                                                               >  1:136327/1‑100 (MQ=255)

                                    tCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCa                                                               >  1:2531896/1‑100 (MQ=255)

                                      ccTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGt                                                             >  2:4118476/1‑100 (MQ=255)

                                      ccTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGt                                                             >  2:331734/1‑100 (MQ=255)

                                       cTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAAACAGACAGTAATTCTTTTTACCGCGACGCAGTa                                                            >  1:3251470/1‑100 (MQ=255)

                                        tGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTaa                                                           <  2:773471/100‑1 (MQ=255)

                                        tGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTaa                                                           >  2:1654275/1‑100 (MQ=255)

                                        tGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTa                                                            >  1:2793857/1‑99 (MQ=255)

                                         gTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                          <  1:2534466/100‑1 (MQ=255)

                                          tataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCGACGTAATGTATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTa                                                            <  1:4069813/97‑1 (MQ=255)

                                           ataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGt                                                        >  2:1889849/1‑100 (MQ=255)

                                           ataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGt                                                        <  2:2841371/100‑1 (MQ=255)

                                            taAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTa                                                       >  1:1534200/1‑100 (MQ=255)

                                            taAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTa                                                       >  2:3914791/1‑100 (MQ=255)

                                            taAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTa                                                       >  2:654743/1‑100 (MQ=255)

                                             aaaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaa                                                      >  2:2044399/1‑100 (MQ=255)

                                              aaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCCTTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaa                                                     >  1:3008964/1‑100 (MQ=255)

                                              aaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGCCAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaa                                                     >  1:645051/1‑100 (MQ=255)

                                               aaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAgagatcggaag                                                    >  2:1333542/1‑90 (MQ=255)

                                               aaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaaa                                                    <  1:1956398/100‑1 (MQ=255)

                                                aaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAAc                                                   >  1:2392931/1‑100 (MQ=255)

                                                   cAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAc                                                <  2:2542451/100‑1 (MQ=255)

                                                    aGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGaac                                               <  2:2898925/100‑3 (MQ=255)

                                                    aGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAcc                                               <  1:2197830/100‑1 (MQ=255)

                                                    aGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATCCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAcc                                               <  2:3224019/100‑1 (MQ=255)

                                                     gggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGCCAGTAATTCTTTTTACCGCGACGCAGTACGGTAAAACGAcc                                               >  1:2718296/1‑99 (MQ=255)

                                                     gggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCa                                              >  2:1307881/1‑100 (MQ=255)

                                                     gggggAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCa                                              >  2:2564713/1‑100 (MQ=255)

                                                      ggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCaa                                             >  2:500716/1‑100 (MQ=255)

                                                       gggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCaaa                                            >  2:2261711/1‑100 (MQ=255)

                                                        ggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAAc                                           >  2:3686514/1‑100 (MQ=255)

                                                         gAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACa                                          >  1:2682611/1‑100 (MQ=255)

                                                         gAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACa                                          <  2:2425215/100‑1 (MQ=255)

                                                          aGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAg                                         <  2:2198248/100‑1 (MQ=255)

                                                             gATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACg                                      >  2:1276930/1‑100 (MQ=255)

                                                              agaaatcACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGa                                     >  2:1240667/6‑100 (MQ=255)

                                                              aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGa                                     >  1:3819557/1‑100 (MQ=255)

                                                              aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAATCGACCAAACAGACGa                                     >  2:4200687/1‑100 (MQ=255)

                                                               tttCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAt                                    <  2:2392937/100‑1 (MQ=255)

                                                               tttCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAt                                    >  1:3228183/1‑100 (MQ=255)

                                                               tttCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAt                                    >  1:1497924/1‑100 (MQ=255)

                                                                ttCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtc                                   >  1:1875348/1‑100 (MQ=255)

                                                                 tctcACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtct                                  <  2:464212/100‑1 (MQ=255)

                                                                  ctcACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtctt                                 >  1:3757315/1‑100 (MQ=255)

                                                                    cACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttct                               >  2:4013385/1‑100 (MQ=255)

                                                                        cccccTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATcttcttctt                           >  1:2871256/1‑100 (MQ=255)

                                                                        cccccTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATcttcttctt                           <  2:2445987/100‑1 (MQ=255)

                                                                          cccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTa                         <  2:1801442/100‑1 (MQ=255)

                                                                            cTTTCCATTTAATCCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCTACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaa                        <  1:2001359/99‑1 (MQ=255)

                                                                             tttCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaaan                      <  2:745653/100‑2 (MQ=255)

                                                                              ttCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaaga                     <  1:2595234/100‑1 (MQ=255)

                                                                               tCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGCTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa                    <  1:2294033/100‑1 (MQ=255)

                                                                                ccACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaag                   <  1:2901653/100‑1 (MQ=255)

                                                                                 cACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGt                  >  2:2360029/1‑100 (MQ=255)

                                                                                  aCTTAATGCATCATTCCCAGCAAATCAGACAGGAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTa                 >  1:1564313/1‑100 (MQ=255)

                                                                                     tAATGCATTATTTCCAGCAAATAAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTc              >  1:2484948/1‑100 (MQ=255)

                                                                                     tAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTc              >  2:1591723/1‑100 (MQ=255)

                                                                                     tAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTc              >  2:1268319/1‑100 (MQ=255)

                                                                                       aTGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAg            <  1:2705388/100‑1 (MQ=255)

                                                                                       aTGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAg            <  1:2072908/100‑1 (MQ=255)

                                                                                        tGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCa             >  2:3483009/1‑98 (MQ=255)

                                                                                        tGCATTATGCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGACGTATTCAg            <  1:2933616/99‑1 (MQ=255)

                                                                                         gCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaagaagatcgtc          <  2:2181892/100‑10 (MQ=255)

                                                                                         gCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGa          >  2:3646104/1‑100 (MQ=255)

                                                                                          cATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGAt         <  2:278702/100‑1 (MQ=255)

                                                                                           attattTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGAt         <  2:462013/99‑1 (MQ=255)

                                                                                           attattCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATc        <  1:240986/100‑1 (MQ=255)

                                                                                            ttattTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCg       <  1:474138/100‑1 (MQ=255)

                                                                                            ttattCCCAGCAAATCAGACCGTAATTCTTTTTACAGCGACGCAGTAATGTAAAACGACCAAACAGACGATCTTCTTCTTTAACGAAGTATTCCGGATCg       <  1:1536802/100‑1 (MQ=255)

                                                                                             tattCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCgg      <  2:2516686/100‑1 (MQ=255)

                                                                                               tttCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGAc    >  2:475431/1‑100 (MQ=255)

                                                                                                 tCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTg  >  2:1189305/1‑100 (MQ=255)

CGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTG  >  NC_000913/1715856‑1716052

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