BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,715,953	T→C	27.4%	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	27.4%	171.9 / 83.3	120	K424E (AAA→GAA)	tyrS	tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base T (42/45); new base C (17/16); total (59/61)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 123 total aligned reads displayed.

AACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTG  >  NC_000913/1715854‑1716052

aacGTGAGACTGGATAGCGCGAATATTTTTCATCATCTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCGCTCCCCCTTTCCACTTAATGCATTAttt                                                                                                     <  2:1956561/100‑1 (MQ=255)

aacGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTAttt                                                                                                     <  2:3939964/100‑1 (MQ=255)

  cGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTcc                                                                                                   <  1:4098338/100‑1 (MQ=255)

  cGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTcc                                                                                                   <  1:480025/100‑1 (MQ=255)

  cGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTcc                                                                                                   <  1:1520183/100‑1 (MQ=255)

   gTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCa                                                                                                  <  2:2702171/100‑1 (MQ=255)

    tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAg                                                                                                 >  1:59840/1‑100 (MQ=255)

    tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAg                                                                                                 >  1:135356/1‑100 (MQ=255)

     gagaCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGc                                                                                                <  1:187791/100‑1 (MQ=255)

     gagaCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGc                                                                                                >  1:3229592/1‑100 (MQ=255)

         cTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAAt                                                                                            <  2:2131264/100‑1 (MQ=255)

         cTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAAt                                                                                            >  2:3024969/1‑100 (MQ=255)

          tGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCGCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATc                                                                                           <  2:2320666/100‑1 (MQ=255)

           ggATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCAACTTAATGCATTATTTCCAGCAAATCa                                                                                          >  2:1294083/1‑100 (MQ=255)

           ggATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAATTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATc                                                                                           >  1:3366662/1‑99 (MQ=255)

              tAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGAc                                                                                       <  2:3266840/100‑1 (MQ=255)

              tAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGAc                                                                                       <  1:524310/100‑1 (MQ=255)

                gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAg                                                                                     >  2:3316117/1‑100 (MQ=255)

                  gagaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCAACTTAATGCATTATTTCCAGCAAATCAGACAGTa                                                                                   >  2:2397929/1‑100 (MQ=255)

                   agaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTa                                                                                   <  2:126318/99‑1 (MQ=255)

                      atatTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTc                                                                               <  1:3476096/100‑1 (MQ=255)

                       tatTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTc                                                                               <  1:3554283/99‑1 (MQ=255)

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

                          ttttCATTATGTTCCCTGTATAGAAACCAGGGGGAGTGATTTCTCACGCCGCCTTTCCACGTAAGGCATTATTCCCAGCAAATCAGACAGTAATTc                                                                               <  2:3900228/96‑1 (MQ=255)

                           tttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCttttt                                                                          <  1:2512183/100‑1 (MQ=255)

                            ttCATCATGTTCCCTGTATAAAACCCAGGGGGAGTGATTTCTCACTCCCTCTTTCCTCTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTa                                                                         <  2:1724217/100‑1 (MQ=255)

                            ttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTa                                                                         >  1:3624398/1‑100 (MQ=255)

                               atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGAAAGTAATTCTTTTTAcc                                                                       >  2:4126287/1‑99 (MQ=255)

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

                                 catGTTCCCTGTATAAAAACCAGGGGGAGCGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACcgcg                                                                    <  1:2317060/100‑1 (MQ=255)

                                   tGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGAc                                                                  >  1:4080867/1‑100 (MQ=255)

                                    gTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACg                                                                 >  1:1726626/1‑100 (MQ=255)

                                     ttCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACg                                                                 <  2:862529/99‑1 (MQ=255)

                                      tCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCa                                                               <  1:3306985/100‑1 (MQ=255)

                                       cccTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAg                                                              >  2:3138612/1‑100 (MQ=255)

                                       cccTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTAccttactgcgt                                                              >  2:1153018/1‑91 (MQ=255)

                                       cccTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAg                                                              <  2:2494997/100‑1 (MQ=255)

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

                                          tGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTaa                                                           <  1:3635230/100‑1 (MQ=255)

                                            tataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAgg                                                         <  2:32589/100‑1 (MQ=255)

                                            tataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAgg                                                         >  1:121672/1‑100 (MQ=255)

                                            tataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAgg                                                         >  1:2858077/1‑100 (MQ=255)

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

                                              taAAAACCAGGGGGAGTGATTTCTCACTCCCCTTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTa                                                       >  1:3960816/1‑100 (MQ=255)

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

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

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

                                               aaaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaa                                                      >  1:1673956/1‑100 (MQ=255)

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

                                               aaaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCAACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaa                                                      >  2:4052576/1‑100 (MQ=255)

                                                aaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaa                                                     >  1:2749030/1‑100 (MQ=255)

                                                  aaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAAc                                                   <  1:3951262/100‑1 (MQ=255)

                                                  aaCCAGGGGGAGTGATTTCTCACGCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaaa                                                    <  1:220773/99‑1 (MQ=255)

                                                   aCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACg                                                  <  1:709501/100‑1 (MQ=255)

                                                    ccAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGa                                                 >  2:3357246/1‑100 (MQ=255)

                                                     cAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTCAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAc                                                <  2:2191983/100‑1 (MQ=255)

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

                                                       gggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCa                                              >  2:733310/1‑100 (MQ=255)

                                                       gggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCa                                              <  1:1538277/100‑1 (MQ=255)

                                                          ggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAAc                                           <  1:3940874/100‑1 (MQ=255)

                                                           gAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACa                                          >  2:482201/1‑100 (MQ=255)

                                                             gTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGa                                        >  1:868737/1‑100 (MQ=255)

                                                             gTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGa                                        >  1:2626807/1‑100 (MQ=255)

                                                             atGATTTTTCACTCCCCCTGTCCACTTATTGCATGATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAg                                         >  2:1810216/2‑99 (MQ=255)

                                                              tGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCCGCCAGTAATTCTTTTTACCGCGCCGCAGTAAGGTAAAACGACCAACCAGAc                                       >  2:1010123/1‑100 (MQ=255)

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

                                                                   tctcACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtct                                  <  2:3876690/100‑1 (MQ=255)

                                                                   tctcACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtct                                  <  2:2098109/100‑1 (MQ=255)

                                                                   tctcACTCCCCCGTTCCACTTAATGCATTATTCCCAGCAAATCAGACATTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtct                                  <  1:2852686/100‑1 (MQ=255)

                                                                    ctcACTCCCCCTTTCCACTTAATGCAATATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtctt                                 >  1:2539761/1‑100 (MQ=255)

                                                                      cACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttct                               <  2:3366636/100‑1 (MQ=255)

                                                                       aCTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttctt                              <  1:2065258/100‑1 (MQ=255)

                                                                        cTCCCCCTTTCCACTTAATGCCTTATTCCCCGCAAATAAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttct                               <  1:1231339/98‑1 (MQ=255)

                                                                        cTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttc                             >  2:4113180/1‑100 (MQ=255)

                                                                        cTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttc                             <  1:3631292/100‑1 (MQ=255)

                                                                        cTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttc                             >  2:129020/1‑100 (MQ=255)

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

                                                                            cccTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTa                         <  2:3095058/100‑1 (MQ=255)

                                                                             ccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAACCAGACGATCTTATTCTTTaa                        >  1:2425674/1‑100 (MQ=255)

                                                                             ccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaa                        <  1:3063680/100‑1 (MQ=255)

                                                                              cTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaaa                       >  1:2352299/1‑100 (MQ=255)

                                                                              cTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaaa                       >  1:1309137/1‑100 (MQ=255)

                                                                                 tCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa                    <  2:30362/100‑1 (MQ=255)

                                                                                 tCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa                    <  2:4108803/100‑1 (MQ=255)

                                                                                  ccACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaag                   <  2:1432975/100‑1 (MQ=255)

                                                                                  ccACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaag                   <  1:3208344/100‑1 (MQ=255)

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

                                                                                   cACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGt                  >  2:1169830/1‑100 (MQ=255)

                                                                                      ttAATGCATTATTCCCAGCAAATCAGACAGAAATTCATTTTACCGCGACGCTGTAAGATAAACCAACCAAGCAGACGATCTTCTTCTTTAGAGAAGTAtt               >  1:1316376/1‑100 (MQ=255)

                                                                                       tAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTc              <  1:3139453/100‑1 (MQ=255)

                                                                                         aTGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAg            >  1:121325/1‑100 (MQ=255)

                                                                                          tGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAgg           >  1:4075745/1‑100 (MQ=255)

                                                                                           gCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGa          >  2:1659995/1‑100 (MQ=255)

                                                                                           gCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGa          <  1:3761535/100‑1 (MQ=255)

                                                                                            cATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACTGTGACGTAGAAAAGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGa          >  2:888906/1‑99 (MQ=255)

                                                                                            cATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGAt         >  2:2041560/1‑100 (MQ=255)

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

                                                                                                attTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCg       <  2:731322/98‑1 (MQ=255)

                                                                                                attTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGa     >  2:2967457/1‑100 (MQ=255)

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

AACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTG  >  NC_000913/1715854‑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: `‑`