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.6%	1683.2 / ‑4.1	523	A53A (GCA→GCG)	yjfY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (254/267); total (256/267)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 538 total aligned reads displayed.

CATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAT  >  NC_000913/4422121‑4422300

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

 atatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCt                                                                                        <  1:4254023/93‑1 (MQ=255)

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

  tatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCtt                                                                                       <  2:5762044/93‑1 (MQ=255)

  tatTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCtt                                                                                       <  1:4608890/93‑1 (MQ=255)

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

    ttGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCg                                                                                     >  2:7772442/1‑93 (MQ=255)

     tGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGc                                                                                    <  1:9806321/93‑1 (MQ=255)

     tGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGc                                                                                    <  2:7280326/93‑1 (MQ=255)

        cGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCttt                                                                                 <  1:6576781/93‑1 (MQ=255)

         gTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTc                                                                                >  1:1765306/1‑93 (MQ=255)

         gTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTc                                                                                <  1:5727898/93‑1 (MQ=255)

          ttACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCa                                                                               >  2:4991243/1‑93 (MQ=255)

           tccTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAg                                                                              <  2:706124/91‑1 (MQ=255)

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

            aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGt                                                                             >  2:6458371/1‑93 (MQ=255)

            aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGt                                                                             <  2:5176437/93‑1 (MQ=255)

                   ggTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGaaa                                                                      <  2:3373268/93‑1 (MQ=255)

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

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

                     ttCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAtt                                                                    >  1:8509683/1‑93 (MQ=255)

                      tCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTa                                                                   <  2:3038374/93‑1 (MQ=255)

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

                          gCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAAGTAtgat                                                               >  1:4890671/1‑93 (MQ=255)

                           cATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgatt                                                              >  2:1749541/1‑93 (MQ=255)

                            aTCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgattg                                                             >  1:3011912/1‑93 (MQ=255)

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

                             tCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTAtgactga                                                            >  2:1312694/1‑93 (MQ=255)

                               aGAATTACGTGGTAGTACGTTGCGCCTTGCTCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGata                                                          <  1:2598897/93‑1 (MQ=255)

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

                               aGAATTACGTGGTAGTACGTTGCGCCTTGCGCAACTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGata                                                          <  1:8937069/93‑1 (MQ=255)

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

                                  aTTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaa                                                       <  2:8794494/93‑1 (MQ=255)

                                  aTTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaa                                                       >  2:3253866/1‑93 (MQ=255)

                                  aTTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaa                                                       <  2:4401035/93‑1 (MQ=255)

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

                                    tACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt                                                     >  2:956855/1‑93 (MQ=255)

                                    tACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt                                                     <  1:6145401/93‑1 (MQ=255)

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

                                      cGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCa                                                   >  2:4457947/1‑93 (MQ=255)

                                      cGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCa                                                   <  1:6001135/93‑1 (MQ=255)

                                      cGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCa                                                   <  1:1859161/93‑1 (MQ=255)

                                      cGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCa                                                   <  2:9357194/93‑1 (MQ=255)

                                       gTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCAc                                                  >  2:2419504/1‑93 (MQ=255)

                                        tGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACg                                                 >  1:2821423/1‑93 (MQ=255)

                                        tGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACg                                                 >  1:9602206/1‑93 (MQ=255)

                                        tGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACg                                                 >  2:4787487/1‑93 (MQ=255)

                                          gtagtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGcc                                               >  2:4283518/1‑93 (MQ=255)

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

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

                                             gtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTaa                                            <  2:9246495/93‑1 (MQ=255)

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

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

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

                                                  ttGCGCGTTGCGCATCTGTTTCTTCATTGAATGCCTGACGCGCTTCGCGTTCAGTGGCGAAATTATGATTGAGATAAATCACGCCTAAGCttt                                       <  1:6744872/93‑1 (MQ=255)

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

                                                     cgcCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGca                                    <  1:6359948/93‑1 (MQ=255)

                                                      gcCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGcac                                   >  2:7587507/1‑93 (MQ=255)

                                                        cTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACAt                                 >  2:5410529/1‑93 (MQ=255)

                                                        cTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACAt                                 >  2:5554435/1‑93 (MQ=255)

                                                         ttGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATc                                <  1:5240110/93‑1 (MQ=255)

                                                         ttGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATc                                <  2:2657210/93‑1 (MQ=255)

                                                         ttGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATc                                >  2:6674548/1‑93 (MQ=255)

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

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

                                                              cATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCa                           <  2:833256/93‑1 (MQ=255)

                                                              cATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCa                           >  2:5107875/1‑93 (MQ=255)

                                                               aTCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCat                          <  1:7802358/93‑1 (MQ=255)

                                                                tCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCata                         >  1:2490633/1‑93 (MQ=255)

                                                                tCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCata                         >  2:7374341/1‑93 (MQ=255)

                                                                tCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTACatt                         >  1:2283999/1‑92 (MQ=255)

                                                                tCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGACTGATATAAATCACGCCTAAGCTTTGCACATCGTCCata                         >  1:7875311/1‑93 (MQ=255)

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

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

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

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

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

                                                                    tttCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTc                     <  1:10010274/93‑1 (MQ=255)

                                                                    tttCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTc                     <  1:5431896/93‑1 (MQ=255)

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

                                                                     ttcttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCt                    >  2:1029152/1‑93 (MQ=255)

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

                                                                        ttcATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGc                 >  1:9523220/1‑93 (MQ=255)

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

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

                                                                          cATTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCt               >  1:8281619/1‑93 (MQ=255)

                                                                          cATTTAAGGCCTGACGCGCTTCGCTTTTAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCTCATCGTCCATATTTCTGGCCt               <  2:7682612/93‑1 (MQ=255)

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

                                                                           aTTTAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtg              >  1:7449616/1‑93 (MQ=255)

                                                                            tttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtgg             <  1:9037987/93‑1 (MQ=255)

                                                                            tttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtgg             >  1:7152188/1‑93 (MQ=255)

                                                                             ttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggt            >  2:7923810/1‑93 (MQ=255)

                                                                             ttAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggt            >  2:1481966/1‑93 (MQ=255)

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

                                                                              tAATGCCTGACGCGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCtggtg           <  2:4921962/93‑1 (MQ=255)

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

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

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

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

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

CATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGCCTAAGCTTTGCACATCGTCCATATTTCTGGCCTGGTGGTTATTAAT  >  NC_000913/4422121‑4422300

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