BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,137,208	A→G	T27T (ACT→ACC)	talC ←	putative transaldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,208	0	A	G	97.6%	508.0 / 1.6	167	T27T (ACT→ACC)	talC	putative transaldolase
Reads supporting (aligned to +/- strand): ref base A (2/2); new base G (77/86); total (79/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

Only 100 of 172 total aligned reads displayed.

AAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCT  >  NC_000913/4137118‑4137298

aaaCAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGt                                                                                          <  1:126822/93‑1 (MQ=255)

aaaCAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGt                                                                                          <  1:1950513/93‑1 (MQ=255)

 aaCAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGggttggtt                                                                                         <  1:2103357/93‑1 (MQ=255)

   cAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGt                                                                                       <  1:4007979/93‑1 (MQ=255)

   cAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGt                                                                                       <  2:135994/93‑1 (MQ=255)

     gAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCa                                                                                     >  2:2593314/1‑93 (MQ=255)

     gAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGGCa                                                                                     >  1:3689084/1‑93 (MQ=255)

       aTGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCAcc                                                                                   >  1:2354704/1‑93 (MQ=255)

        tGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCAccc                                                                                  <  2:187697/93‑1 (MQ=255)

          cccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCAccctg                                                                                <  2:2769918/93‑3 (MQ=255)

          cccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCAccctg                                                                                <  1:879467/93‑3 (MQ=255)

           ccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTgg                                                                               >  2:1233532/1‑93 (MQ=255)

            cTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGc                                                                              <  1:1742575/93‑1 (MQ=255)

            cTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGc                                                                              >  1:3783907/1‑93 (MQ=255)

              catcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCaa                                                                            <  2:436630/93‑1 (MQ=255)

              catcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCaa                                                                            >  2:1474235/1‑93 (MQ=255)

               atcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAAt                                                                           >  2:2165229/1‑93 (MQ=255)

               atcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAAt                                                                           <  1:3123978/93‑1 (MQ=255)

               atcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAAt                                                                           >  1:3819643/1‑93 (MQ=255)

                tcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATg                                                                          >  2:3074286/1‑93 (MQ=255)

                  aCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATggg                                                                        <  1:2169731/93‑1 (MQ=255)

                   ccAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATgggg                                                                       <  1:4361656/93‑1 (MQ=255)

                   ccAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATgggg                                                                       >  1:2818258/1‑93 (MQ=255)

                    cAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGa                                                                      >  2:2422914/1‑93 (MQ=255)

                    cAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGa                                                                      <  1:3227308/93‑1 (MQ=255)

                    cAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGa                                                                      >  1:4138014/1‑93 (MQ=255)

                     aaTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGaa                                                                     <  2:708393/93‑1 (MQ=255)

                      aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAat                                                                    <  2:209311/93‑1 (MQ=255)

                       tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                   >  1:4439281/1‑93 (MQ=255)

                       tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                   >  2:2211990/1‑93 (MQ=255)

                       tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                   <  1:11285/93‑1 (MQ=255)

                       tCGCTGTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                   <  2:2176266/93‑1 (MQ=255)

                        cGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                  >  1:2726336/1‑93 (MQ=255)

                        cGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                  <  2:3958397/93‑1 (MQ=255)

                        cGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                  <  2:899915/93‑1 (MQ=255)

                        cGCTTTTTGCAGAAGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                  <  1:1583106/93‑1 (MQ=255)

                         gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                 <  2:2524219/93‑1 (MQ=255)

                         gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                 >  2:1802049/1‑93 (MQ=255)

                           tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCg                                                               >  2:3811493/1‑93 (MQ=255)

                           tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCg                                                               >  2:1280992/1‑93 (MQ=255)

                           tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCg                                                               <  1:1569293/93‑1 (MQ=255)

                            ttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGt                                                              <  1:53654/93‑1 (MQ=255)

                            ttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGt                                                              <  1:732410/93‑1 (MQ=255)

                               tGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGcc                                                           <  2:1215840/93‑1 (MQ=255)

                                gCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGccg                                                          >  1:1408929/1‑92 (MQ=255)

                                  aGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGa                                                        <  1:2475258/93‑1 (MQ=255)

                                   gACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGAc                                                       <  1:697536/93‑1 (MQ=255)

                                    aCGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACg                                                      <  2:1130995/93‑1 (MQ=255)

                                     cgcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGt                                                     <  1:228448/93‑1 (MQ=255)

                                     cgcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGt                                                     <  1:1836511/93‑1 (MQ=255)

                                      gcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGtt                                                    <  1:2828304/93‑1 (MQ=255)

                                      gcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGtt                                                    <  1:1987028/93‑1 (MQ=255)

                                       cgtcagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTc                                                   <  2:1265575/90‑1 (MQ=255)

                                        ggCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCg                                                  <  2:212258/93‑1 (MQ=255)

                                          cagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGAc                                                <  2:2104825/93‑1 (MQ=255)

                                            gcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACtt                                              <  1:1022526/93‑1 (MQ=255)

                                              aCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCt                                            >  2:4157846/1‑93 (MQ=255)

                                                ttCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                          <  2:919727/93‑1 (MQ=255)

                                                ttCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                          >  2:790694/1‑93 (MQ=255)

                                                ttCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                          <  2:3337864/93‑1 (MQ=255)

                                                   ccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGAc                                       <  2:1309245/93‑1 (MQ=255)

                                                    cATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACg                                      >  1:1040031/1‑93 (MQ=255)

                                                    cATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACg                                      >  2:2374440/1‑93 (MQ=255)

                                                     atatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTTTGCGACGt                                     <  2:3839589/93‑1 (MQ=255)

                                                      tatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGtt                                    >  1:2016213/1‑93 (MQ=255)

                                                      tatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGtt                                    <  2:4661547/93‑1 (MQ=255)

                                                       atGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGGCGTTa                                   >  2:1709696/1‑93 (MQ=255)

                                                       atGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTa                                   >  1:1367261/1‑93 (MQ=255)

                                                       atGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTa                                   <  2:2932738/93‑1 (MQ=255)

                                                        tGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAg                                  <  2:3365971/93‑1 (MQ=255)

                                                        tGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAg                                  >  1:3324840/1‑93 (MQ=255)

                                                         ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGc                                 >  1:881545/1‑93 (MQ=255)

                                                         ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGc                                 <  1:1984158/93‑1 (MQ=255)

                                                         ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGc                                 <  2:2743588/93‑1 (MQ=255)

                                                                ttGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTcc                          >  2:321783/1‑93 (MQ=255)

                                                                ttGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTcc                          >  2:3132904/1‑93 (MQ=255)

                                                                 tGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCa                         <  1:575689/93‑1 (MQ=255)

                                                                   cTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGa                       >  1:3868998/1‑93 (MQ=255)

                                                                    tGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGAt                      >  1:59875/1‑93 (MQ=255)

                                                                    tGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGAt                      >  1:3489320/1‑93 (MQ=255)

                                                                    tGGCAGCGATAATGCGCGGGTTGGTTGTCACCCTGGCAATGAGGAATATGCGTGCCAGACGTTCGACTTCTTCGACGTTAGCGGTGTCCAgtt                      >  2:1667136/1‑91 (MQ=255)

                                                                      gCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATAc                    >  2:4406430/1‑93 (MQ=255)

                                                                       cAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACa                   <  1:2553566/93‑1 (MQ=255)

                                                                        aGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAg                  <  2:3710117/93‑1 (MQ=255)

                                                                        aGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAg                  >  1:1876549/1‑93 (MQ=255)

                                                                         gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                 <  2:3466618/93‑1 (MQ=255)

                                                                         gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                 >  1:4675514/1‑93 (MQ=255)

                                                                         gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                 >  1:1800684/1‑93 (MQ=255)

                                                                           gATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTc               <  1:502806/93‑1 (MQ=255)

                                                                             tAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCa             >  2:1619191/1‑93 (MQ=255)

                                                                             tAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGATCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTTCa             >  2:2812083/1‑93 (MQ=255)

                                                                              aaTGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCat            >  2:2137492/1‑93 (MQ=255)

                                                                               aTGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatg           <  1:2583035/93‑1 (MQ=255)

                                                                               aTGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatg           >  1:1002453/1‑93 (MQ=255)

                                                                               aTGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatg           >  1:980730/1‑93 (MQ=255)

                                                                               aTGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatg           <  2:4602489/93‑1 (MQ=255)

                                                                                tGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatga          <  2:3517924/93‑1 (MQ=255)

                                                                                   tCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGt       >  2:1187752/1‑93 (MQ=255)

                                                                                       gttggttgTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTcc   <  2:2607346/93‑1 (MQ=255)

                                                                                        ttAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCt  <  2:4021495/93‑1 (MQ=255)

AAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCT  >  NC_000913/4137118‑4137298

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