BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,644,167	T→G	E339A (GAA→GCA)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,644,167	0	T	G	97.5%	389.1 / ‑0.7	122	E339A (GAA→GCA)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand): ref base T (3/0); new base G (58/61); total (61/61)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

Only 100 of 145 total aligned reads displayed.

TAACCACTTCAAACAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTTCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTG  >  NC_000913/3644072‑3644263

  aCCACTTCAAACAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTgt                                                                                            >  1:307763/1‑100 (MQ=255)

   ccACTTCAAACAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTgta                                                                                           <  1:1507346/100‑1 (MQ=255)

    cACTTCAAACAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTgtag                                                                                          <  1:276765/100‑1 (MQ=255)

    cACTTCAAACAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTgtag                                                                                          <  2:1108924/100‑1 (MQ=255)

     aCTTCAAACAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTgtagt                                                                                         >  2:2939005/1‑100 (MQ=255)

         cAAACAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGc                                                                                     <  1:783513/100‑1 (MQ=255)

         cAAACAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCCGTAGTAAGc                                                                                     <  1:2090248/100‑1 (MQ=255)

          aaaCAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCg                                                                                    >  1:2281/1‑100 (MQ=255)

             cAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGata                                                                                 <  1:2051333/100‑1 (MQ=255)

               ggCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATc                                                                               >  1:3636947/1‑100 (MQ=255)

               ggCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATc                                                                               <  1:446016/100‑1 (MQ=255)

               ggCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATc                                                                               >  1:2770303/1‑100 (MQ=255)

               ggCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATc                                                                               >  1:1615138/1‑100 (MQ=255)

               ggCCGTTAACCGCTTTGTTATGCGGGAAGTACGGACGCAGCTATTCGTCACTGATGCTGTAGCGGTGATGTTTTTGTTTTGCGCTGTAGTAAGCGatagc                                                                               >  1:3478302/1‑98 (MQ=255)

                 ccGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATccc                                                                             <  2:316928/100‑1 (MQ=255)

                  cGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCa                                                                            >  1:2359919/1‑100 (MQ=255)

                   gTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCAc                                                                           >  1:312027/1‑100 (MQ=255)

                   gTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCAc                                                                           <  2:2446750/100‑1 (MQ=255)

                   gTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCAc                                                                           <  1:2194974/100‑1 (MQ=255)

                   gTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCAc                                                                           <  1:639163/100‑1 (MQ=255)

                   gTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCAc                                                                           <  1:2607293/100‑1 (MQ=255)

                      aaCCGCTTTGTTTTCCGGGAAGTACTGACGCAGCTGTCCGTCACTGATGTTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCCCGGc                                                                        <  2:1309866/100‑1 (MQ=255)

                       aCCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCt                                                                       >  1:1699494/1‑100 (MQ=255)

                       aCCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCt                                                                       >  1:2786782/1‑100 (MQ=255)

                        ccGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTg                                                                      >  2:3889422/1‑100 (MQ=255)

                          gCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCa                                                                    <  2:672074/100‑1 (MQ=255)

                          gCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCa                                                                    <  2:671516/100‑1 (MQ=255)

                          gCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCa                                                                    <  2:2937087/100‑1 (MQ=255)

                          gCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTATTGCGCTGTAGTAAGCGATATCCCACGGCTGCa                                                                    >  1:4038153/1‑100 (MQ=255)

                             ttGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACt                                                                 <  1:909335/100‑1 (MQ=255)

                             ntGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGCGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACt                                                                 >  1:1187573/2‑100 (MQ=255)

                               gTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCa                                                               <  1:3945720/100‑1 (MQ=255)

                                ttttCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCAt                                                              <  1:1497422/100‑1 (MQ=255)

                                ttttCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCAt                                                              >  1:1667381/1‑100 (MQ=255)

                                  ttCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCg                                                            <  1:4001905/100‑1 (MQ=255)

                                     cGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATGTCCCACGGCTGCAACTCATCGACg                                                         <  2:2932226/100‑1 (MQ=255)

                                     cGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACg                                                         <  1:2982479/100‑1 (MQ=255)

                                     cGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACg                                                         <  1:476472/100‑1 (MQ=255)

                                      gggAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGc                                                        >  1:310366/1‑100 (MQ=255)

                                       ggAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGcc                                                       <  1:780081/100‑1 (MQ=255)

                                          aGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCaaa                                                    >  2:804044/1‑100 (MQ=255)

                                           gTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAAt                                                   >  2:1126823/1‑100 (MQ=255)

                                             agggACGCCGCTCTCCGTAGCTGATGCTGTAGAGGTGCTGTTTTTGTTTTGAGCTGTAGTGAGCGCCACCCCACGGCTGCAACTCATCGACGCCAAATTc                                                 <  1:3319653/98‑1 (MQ=255)

                                             aCGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTc                                                 >  2:3959590/1‑100 (MQ=255)

                                             aCGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTACTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTc                                                 >  2:1540508/1‑100 (MQ=255)

                                              cGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCg                                                <  1:134816/100‑1 (MQ=255)

                                              cGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCg                                                >  2:3787128/1‑100 (MQ=255)

                                               ggACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCgg                                               <  2:306436/100‑1 (MQ=255)

                                                gACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGc                                              <  1:1967594/100‑1 (MQ=255)

                                                gACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGc                                              <  1:2581911/100‑1 (MQ=255)

                                                  cGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCtt                                            >  2:1159638/1‑100 (MQ=255)

                                                  cGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCAGCGACGCCAAATTCGGCtt                                            >  1:3823239/1‑100 (MQ=255)

                                                    cAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTg                                          <  1:1988111/100‑1 (MQ=255)

                                                     aGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTgg                                         <  1:2275823/100‑1 (MQ=255)

                                                      gCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGc                                        >  2:178902/1‑100 (MQ=255)

                                                      gCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGc                                        <  2:188822/100‑1 (MQ=255)

                                                      gCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGc                                        >  1:292795/1‑100 (MQ=255)

                                                       cTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCa                                       <  1:2547788/100‑1 (MQ=255)

                                                         gTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCaaa                                     >  2:1259444/1‑100 (MQ=255)

                                                          ttCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAg                                    <  1:524967/100‑1 (MQ=255)

                                                            cGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCGCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGc                                  >  2:692923/1‑100 (MQ=255)

                                                            cGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGc                                  >  1:3985564/1‑100 (MQ=255)

                                                              tCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCAc                                <  1:2549930/100‑1 (MQ=255)

                                                               cACTGATGCTGTAGAGGTGCTGTTTTTTTTTTTTTTTTTTGTAAGCGATAGCCCACGGCTGCAACTCCTCGACGCCAGATTCGGCTTTGGCTAAgccccg                               >  2:2075474/1‑95 (MQ=255)

                                                                  tGATGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCaa                            >  2:3631845/1‑100 (MQ=255)

                                                                   gATGCTGTAGAGGTGCTGTTTCTGGTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCTACGCCAAATTCGGCTCTGGCCAAGGCACGCAAt                           >  2:1556989/1‑100 (MQ=255)

                                                                    aTGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAAtt                          <  2:3743726/100‑1 (MQ=255)

                                                                     tGCTGTAGAGGTGCTGTTTTTGTTTTTCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGAAATTg                         >  1:1934833/1‑100 (MQ=255)

                                                                     tGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCGCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTg                         >  2:2657124/1‑100 (MQ=255)

                                                                     tGCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTg                         <  1:3842478/100‑1 (MQ=255)

                                                                      gCTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTgc                        >  2:373190/1‑100 (MQ=255)

                                                                       cTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTgcg                       <  2:146320/100‑1 (MQ=255)

                                                                       cTGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTgcg                       <  1:1118856/100‑1 (MQ=255)

                                                                        tGTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTgcgc                      >  1:4063464/1‑100 (MQ=255)

                                                                         gTAGAGGTGCTGTTTTTGTTTTTCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGcc                     >  2:29341/1‑100 (MQ=255)

                                                                         gTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGcc                     >  2:2193028/1‑100 (MQ=255)

                                                                         gTAGAGGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGcc                     <  1:897016/100‑1 (MQ=255)

                                                                           agagGTGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAg                   <  2:1424245/100‑1 (MQ=255)

                                                                             gggTGCTGTCTTTGTGTCGCGCCGTAGTAATCGAGATCCCGCGGCGGGAGCTGATCGAGGCCAAGTTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGct                 <  2:2032586/99‑1 (MQ=255)

                                                                                tGCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCtt              <  1:2990723/100‑1 (MQ=255)

                                                                                 gCTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCttt             <  1:2932575/100‑1 (MQ=255)

                                                                                  cTGTTTTTGTTTTGTTTTGTAATAAGAGATATCCGACGGCTGCAACTCATCGACGCCAAAGTGGCCCTTGGCAAAGTCGCGTACTTTCGCTATATCtttt            >  2:1042789/1‑100 (MQ=255)

                                                                                  cTGTTTTTGTTTTGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTAGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCtgtt            >  2:2582208/1‑100 (MQ=255)

                                                                                   tgtttttgttttGCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGCCGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCttttt           >  1:3886913/1‑100 (MQ=255)

                                                                                    gtttttgttttgCGCTGTAGTAAGCGAGATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTc          <  1:1380969/100‑1 (MQ=255)

                                                                                    gtgtttgttttgCGCTGTAGTAAGCGCTATCCCACGGCTGCAACTCATCGACGCCAAACTCGGCTTTGGCAAAGACACGCAATTGCGCCAGCTCTTTTTc          <  2:3855555/97‑1 (MQ=255)

                                                                                      ttttgttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCga        >  2:2954924/1‑99 (MQ=255)

                                                                                      ttttgttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGc        >  1:1947936/1‑100 (MQ=255)

                                                                                      ttttgttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGc        >  2:2047073/1‑100 (MQ=255)

                                                                                      ttttgttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGc        >  1:498044/1‑100 (MQ=255)

                                                                                       tttgtttttCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGcc       <  1:29812/100‑1 (MQ=255)

                                                                                       tttgttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGcc       <  2:1671560/100‑1 (MQ=255)

                                                                                       tttgttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGcc       <  1:1557274/100‑1 (MQ=255)

                                                                                        ttgttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCt      <  1:1575576/100‑1 (MQ=255)

                                                                                         tgttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCtt     >  1:1139763/1‑100 (MQ=255)

                                                                                          gttttgCGCTGTAGTAAGCGATATCCCTCGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTtg    <  2:243184/100‑1 (MQ=255)

                                                                                          gttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTtg    <  2:3954821/100‑1 (MQ=255)

                                                                                          gttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTtg    >  1:2328471/1‑100 (MQ=255)

                                                                                          gttttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTtg    <  2:131572/100‑1 (MQ=255)

                                                                                            tttgCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTtgtg  <  2:3522055/100‑1 (MQ=255)

TAACCACTTCAAACAGGCCGTTAACCGCTTTGTTTTCCGGGAAGTACGGACGCAGCTGTTCGTCACTGATGCTGTAGAGGTGCTGTTTTTGTTTTTCGCTGTAGTAAGCGATATCCCACGGCTGCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTG  >  NC_000913/3644072‑3644263

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: `‑`