BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,422,211	T→C	A53A (GCA→GCG)	yjfY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,422,211	0	T	C	99.2%	830.8 / ‑3.6	257	A53A (GCA→GCG)	yjfY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (133/122); total (135/122)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 265 total aligned reads displayed.

CATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAT  >  NC_000913/4422121‑4422300

cATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACgcgc                                                                                         <  2:533666/93‑1 (MQ=255)

cATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACgcgc                                                                                         <  2:1364344/93‑1 (MQ=255)

cATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACgcgc                                                                                         >  1:800811/1‑93 (MQ=255)

 atatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCt                                                                                        >  2:3459652/1‑93 (MQ=255)

 atatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCt                                                                                        >  1:1896535/1‑93 (MQ=255)

  tatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCtt                                                                                       >  2:2943952/1‑93 (MQ=255)

  tatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCtt                                                                                       >  2:1813990/1‑93 (MQ=255)

  tatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCtt                                                                                       >  1:2427707/1‑93 (MQ=255)

           tACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAg                                                                              <  2:613773/93‑1 (MQ=255)

           tACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAg                                                                              >  1:3860995/1‑93 (MQ=255)

            tcTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGt                                                                             <  2:1490098/92‑1 (MQ=255)

             cTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTg                                                                            <  2:4736521/93‑1 (MQ=255)

             cTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTg                                                                            <  1:3449302/93‑1 (MQ=255)

              tCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTgg                                                                           >  1:4116328/1‑93 (MQ=255)

               ccccGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGc                                                                          <  1:2531074/93‑1 (MQ=255)

                cccGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCg                                                                         >  1:1092160/1‑93 (MQ=255)

                 ccGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGTGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGa                                                                        >  2:2010240/1‑93 (MQ=255)

                  cGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGaa                                                                       >  1:4671883/1‑93 (MQ=255)

                    gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAt                                                                     >  2:1706728/1‑93 (MQ=255)

                    gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAt                                                                     <  1:984926/93‑1 (MQ=255)

                    gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAt                                                                     >  1:753992/1‑93 (MQ=255)

                    gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAt                                                                     >  1:473040/1‑93 (MQ=255)

                       cccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAt                                                                  <  1:2519465/93‑1 (MQ=255)

                       cccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAt                                                                  >  2:4296708/1‑93 (MQ=255)

                         cGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtga                                                                <  1:2566458/93‑1 (MQ=255)

                         cGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtga                                                                <  1:51473/93‑1 (MQ=255)

                          gCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgat                                                               <  1:4293320/93‑1 (MQ=255)

                           cATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgatt                                                              <  1:1462481/93‑1 (MQ=255)

                            aTCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtgattg                                                             >  1:2133982/1‑93 (MQ=255)

                             tCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattga                                                            <  2:982476/93‑1 (MQ=255)

                             tCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattga                                                            <  2:3498444/93‑1 (MQ=255)

                             tCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattga                                                            <  1:4604027/93‑1 (MQ=255)

                             tCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattga                                                            >  1:272918/1‑93 (MQ=255)

                              cAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattgat                                                           <  1:1542782/93‑1 (MQ=255)

                               aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGTGAAATTATGATTGata                                                          <  1:1254371/93‑1 (MQ=255)

                               aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGata                                                          <  2:3688368/93‑1 (MQ=255)

                                gAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGatat                                                         >  2:2851902/1‑93 (MQ=255)

                                gAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGatat                                                         >  2:2304567/1‑93 (MQ=255)

                                 aaTTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGAtata                                                        >  2:1688108/1‑93 (MQ=255)

                                   ttACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaaa                                                      <  1:2302864/93‑1 (MQ=255)

                                   ttACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaaa                                                      >  1:1353769/1‑93 (MQ=255)

                                   ttACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaaa                                                      <  1:3300509/93‑1 (MQ=255)

                                     aCGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATc                                                    <  2:2196272/93‑1 (MQ=255)

                                        tGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACg                                                 <  1:2679729/93‑1 (MQ=255)

                                         ggTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc                                                >  2:81283/1‑93 (MQ=255)

                                          gtagtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGcc                                               <  2:3724210/93‑1 (MQ=255)

                                           tagtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCt                                              <  2:752935/93‑1 (MQ=255)

                                            agtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGCAATTATGATTGATATAAATCACGCCTa                                             >  1:3168505/1‑93 (MQ=255)

                                            agtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTa                                             <  1:4488830/93‑1 (MQ=255)

                                              taCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAg                                           <  2:3243769/93‑1 (MQ=255)

                                               aCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGc                                          >  1:3057357/1‑93 (MQ=255)

                                                cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCt                                         <  2:156074/93‑1 (MQ=255)

                                                cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCt                                         >  2:2036254/1‑93 (MQ=255)

                                                 gTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCtt                                        <  2:2560291/93‑1 (MQ=255)

                                                 gTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCtt                                        <  2:110686/93‑1 (MQ=255)

                                                  ttGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCttt                                       <  1:2879607/93‑1 (MQ=255)

                                                  ttGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCttt                                       >  1:2755673/1‑93 (MQ=255)

                                                  ttGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCttt                                       >  1:4037226/1‑93 (MQ=255)

                                                  ttGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCttt                                       <  2:3310449/93‑1 (MQ=255)

                                                    gcgcCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGc                                     >  1:392816/1‑93 (MQ=255)

                                                       ccTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGcaca                                  >  2:3813738/1‑93 (MQ=255)

                                                       ccTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGcaca                                  >  1:840843/1‑93 (MQ=255)

                                                        cTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACAt                                 <  2:2650795/93‑1 (MQ=255)

                                                           gcgcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGt                              <  2:105810/93‑1 (MQ=255)

                                                            cgcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTc                             <  1:655693/93‑1 (MQ=255)

                                                            cgcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTc                             <  1:3185277/93‑1 (MQ=255)

                                                             gcATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTcc                            >  2:2592258/1‑93 (MQ=255)

                                                               aTCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCat                          >  1:3461919/1‑93 (MQ=255)

                                                               aTCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCat                          >  1:3580009/1‑93 (MQ=255)

                                                                tCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCatt                         >  1:2217817/1‑92 (MQ=255)

                                                                  tGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATAtt                       >  1:69465/1‑93 (MQ=255)

                                                                  tGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATAtt                       <  1:250696/93‑1 (MQ=255)

                                                                  tGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATAtt                       <  1:4082096/93‑1 (MQ=255)

                                                                   gTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATAttt                      >  1:2645069/1‑93 (MQ=255)

                                                                   gTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATAttt                      <  2:4520653/93‑1 (MQ=255)

                                                                    tttCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTc                     >  2:3979164/1‑93 (MQ=255)

                                                                     ttcttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCt                    <  2:211643/93‑1 (MQ=255)

                                                                      tcttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTg                   <  2:3541864/93‑1 (MQ=255)

                                                                      tcttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTg                   >  1:2963657/1‑93 (MQ=255)

                                                                      tcttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGATTTGCACATCGTCCATATTTCTg                   <  2:1860966/93‑1 (MQ=255)

                                                                         tcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGcc                <  2:240050/93‑1 (MQ=255)

                                                                         tcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCATCCATATTTCTGGcc                >  1:3226039/1‑93 (MQ=255)

                                                                          cATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCt               <  2:1163694/93‑1 (MQ=255)

                                                                           aTTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtg              <  1:94635/93‑1 (MQ=255)

                                                                             ttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggt            <  1:4445630/93‑1 (MQ=255)

                                                                               aaTGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggtgg          >  1:4631822/1‑93 (MQ=255)

                                                                               aaTGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggtgg          >  1:1051021/1‑93 (MQ=255)

                                                                                aTGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTggtggt         >  1:1710398/1‑93 (MQ=255)

                                                                                 tGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGtt        <  2:3237658/93‑1 (MQ=255)

                                                                                 tGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGtt        >  1:2612064/1‑93 (MQ=255)

                                                                                 tGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGtt        >  1:1176620/1‑93 (MQ=255)

                                                                                  gCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGtta       >  2:3023996/1‑93 (MQ=255)

                                                                                  gCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGtta       <  1:1848163/93‑1 (MQ=255)

                                                                                   ccTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGttat      <  2:851306/93‑1 (MQ=255)

                                                                                    cTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGttatt     <  1:2781277/93‑1 (MQ=255)

                                                                                    cTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGttatt     >  2:3621168/1‑93 (MQ=255)

                                                                                      gACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTaa   >  2:2251266/1‑93 (MQ=255)

                                                                                       aCGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAt  >  2:4601111/1‑93 (MQ=255)

                                                                                       aCGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAt  <  2:402611/93‑1 (MQ=255)

                                                                                       aCGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAt  >  2:3994513/1‑93 (MQ=255)

CATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAT  >  NC_000913/4422121‑4422300

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