BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,651,556	C→T	78.7%	pseudogene (266/741 nt)	yhiS →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,651,556	0	C	T	78.7%	393.1 / 107.9	216	pseudogene (266/741 nt)	yhiS	pseudogene
Reads supporting (aligned to +/- strand): ref base C (15/30); new base T (96/72); total (114/102)
Fisher's exact test for biased strand distribution p-value = 6.77e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

Only 100 of 219 total aligned reads displayed.

TCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAGATA  >  NC_000913/3651457‑3651655

tctcTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGc                                                                                                     <  1:2713694/100‑1 (MQ=255)

  tctGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTc                                                                                                   <  2:4027033/100‑1 (MQ=255)

  tctGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTc                                                                                                   <  1:1564695/100‑1 (MQ=255)

   ctGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCa                                                                                                  >  2:210653/1‑100 (MQ=255)

    tGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCaa                                                                                                 >  2:3812230/1‑100 (MQ=255)

     gATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAAc                                                                                                <  2:420088/100‑1 (MQ=255)

      aTCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACt                                                                                               <  2:1694012/100‑1 (MQ=255)

        cAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTg                                                                                             >  1:2463995/1‑100 (MQ=255)

          atatGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTa                                                                                           >  1:3579413/1‑100 (MQ=255)

          atatGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTa                                                                                           <  2:557319/100‑1 (MQ=255)

          atatGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTa                                                                                           <  1:2116670/100‑1 (MQ=255)

                tttAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCtt                                                                                     <  1:1128847/100‑1 (MQ=255)

                  tAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTTCATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCttat                                                                                   <  1:2662972/100‑1 (MQ=255)

                  tAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCttat                                                                                   >  1:360715/1‑100 (MQ=255)

                  tAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCttat                                                                                   <  1:2599896/100‑1 (MQ=255)

                   aGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCttatt                                                                                  <  2:2303716/100‑1 (MQ=255)

                   aGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCttatt                                                                                  <  2:1426371/100‑1 (MQ=255)

                   aGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCttatt                                                                                  <  1:1385706/100‑1 (MQ=255)

                    gcgcAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCttatta                                                                                 <  2:1189636/100‑1 (MQ=255)

                     cgcAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTAc                                                                                >  2:1892180/1‑100 (MQ=255)

                     cgcAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTAc                                                                                <  1:1162201/100‑1 (MQ=255)

                     cgcAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTAc                                                                                >  1:1050564/1‑100 (MQ=255)

                      gcAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACt                                                                               <  2:4096624/100‑1 (MQ=255)

                      gcAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACt                                                                               <  1:171312/100‑1 (MQ=255)

                      gcAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTAAACTTGTAAGTCTTATTAc                                                                                >  2:3892145/1‑99 (MQ=255)

                       taaCACCAACTTCAACAAATTAATGCGTCTCACACGAACACTACATATACAAATGCATACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTa                                                                              <  2:703247/99‑1 (MQ=255)

                        aaCACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTaa                                                                             >  2:1352873/1‑100 (MQ=255)

                         acacCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAAc                                                                            <  1:3812259/100‑1 (MQ=255)

                          caccacCTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACa                                                                           <  1:3342688/100‑1 (MQ=255)

                           acCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACAt                                                                          >  1:1334175/1‑100 (MQ=255)

                            ccAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACAtt                                                                         >  1:3669237/1‑100 (MQ=255)

                            ccAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACAtt                                                                         >  1:532589/1‑100 (MQ=255)

                             cAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTg                                                                        <  1:3263779/100‑1 (MQ=255)

                             cAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTg                                                                        <  1:3437065/100‑1 (MQ=255)

                                 ttCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAAc                                                                    <  1:1997371/100‑1 (MQ=255)

                                 ttCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAAc                                                                    <  1:425195/100‑1 (MQ=255)

                                   caacaaATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGa                                                                  >  1:2316455/1‑100 (MQ=255)

                                   caacaaATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGa                                                                  <  1:989141/100‑1 (MQ=255)

                                    aacaaATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGAt                                                                 >  2:2944168/1‑100 (MQ=255)

                                    aacaaATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGAt                                                                 <  1:1746195/100‑1 (MQ=255)

                                    aacaaATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGATCAACTTGTAAGTCTTATTACTAACATTGAAACGAt                                                                 >  1:2922778/1‑100 (MQ=255)

                                     acaaATTAATGCGTGCCACACTAACACTAAATAAACAAAAGAAGAATACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATa                                                                <  2:881923/100‑1 (MQ=255)

                                     acaaATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATa                                                                >  1:2862853/1‑100 (MQ=255)

                                     acaaATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATa                                                                >  2:2908375/1‑100 (MQ=255)

                                      caaATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATaa                                                               >  1:1757564/1‑100 (MQ=255)

                                       aaaTTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAg                                                              <  1:1020567/100‑1 (MQ=255)

                                       aaaTTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAg                                                              <  1:1036143/100‑1 (MQ=255)

                                        aaTTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAg                                                              <  2:3236958/99‑1 (MQ=255)

                                             aTGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGac                                                        <  1:3729668/100‑1 (MQ=255)

                                             aTGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGac                                                        <  1:771542/100‑1 (MQ=255)

                                             aTGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGac                                                        <  2:3521346/100‑1 (MQ=255)

                                             aTGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGac                                                        <  2:3852772/100‑1 (MQ=255)

                                             aTGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGac                                                        <  1:3482732/100‑1 (MQ=255)

                                             aTGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGac                                                        <  1:4066149/100‑1 (MQ=255)

                                                cGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTGAGTCTTATTACAAACATTGAAACGATAAGTTCGAcaca                                                     >  2:1664904/1‑100 (MQ=255)

                                                 gTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACaa                                                    >  1:150341/1‑100 (MQ=255)

                                                   gCCACACGAACACTACATATACAAATGCAGACTCCTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACt                                                  >  2:3352843/1‑100 (MQ=255)

                                                     cacaCGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTa                                                >  2:2035773/1‑100 (MQ=255)

                                                      acacGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTAc                                               <  2:3043420/100‑1 (MQ=255)

                                                       cacGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACg                                              >  1:4074438/1‑100 (MQ=255)

                                                        acGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGc                                             <  1:1065926/100‑1 (MQ=255)

                                                        acGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGc                                             <  2:4074407/100‑1 (MQ=255)

                                                           aaCACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAAc                                          >  1:2281406/1‑100 (MQ=255)

                                                           aaCACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAAc                                          <  2:328789/100‑1 (MQ=255)

                                                             cacTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAg                                        >  1:8957/1‑100 (MQ=255)

                                                              acTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGa                                       >  2:2026814/1‑100 (MQ=255)

                                                                tACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGaca                                     <  2:313554/100‑1 (MQ=255)

                                                                tACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACATACGCAACAg                                        <  2:327498/97‑1 (MQ=255)

                                                                 aCATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGacac                                    >  1:2153595/1‑100 (MQ=255)

                                                                 aCATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGacac                                    >  1:1650626/1‑100 (MQ=255)

                                                                  cATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGAcaca                                   <  1:2853943/100‑1 (MQ=255)

                                                                    tataCAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTATGCAACAGACAcaat                                 >  2:698461/1‑100 (MQ=255)

                                                                     ataCAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatc                                <  1:2560710/100‑1 (MQ=255)

                                                                     ataCAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatc                                >  1:1956637/1‑100 (MQ=255)

                                                                      taCAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatca                               >  1:2238604/1‑100 (MQ=255)

                                                                      taCAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatca                               >  1:3977086/1‑100 (MQ=255)

                                                                      taCAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatca                               >  1:290533/1‑100 (MQ=255)

                                                                        cAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAATAGACAcaatcaa                              >  2:2249396/1‑99 (MQ=255)

                                                                        cAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatcaat                             >  2:2646016/1‑100 (MQ=255)

                                                                        cAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatcaat                             >  2:3728533/1‑100 (MQ=255)

                                                                           aTGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACt                          <  1:3675992/100‑1 (MQ=255)

                                                                            tGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACtt                         >  1:3064128/1‑100 (MQ=255)

                                                                               aGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAAt                      <  1:3144032/100‑1 (MQ=255)

                                                                                gACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATc                     >  1:2788225/1‑100 (MQ=255)

                                                                                 actactCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACATAACTTACGCAACAGACACAATCAATACTTAATc                     <  2:1738120/99‑1 (MQ=255)

                                                                                 actactCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAATTTCGACACAACTTACGCAACAGACCCAATCAATACTTGATc                     <  2:273347/99‑1 (MQ=255)

                                                                                     ctCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGAtt                >  2:241453/1‑100 (MQ=255)

                                                                                      tCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATta               <  2:79234/100‑1 (MQ=255)

                                                                                      tCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATta               <  2:1592748/100‑1 (MQ=255)

                                                                                       cAAAAGATGTAGTGCAACTTGTAAGTCTTATTACTAACATTGAAACGAAAAGTTCGACAAAACTTACGCAAGAGAAACAATCAATACTTAATCAGAAtaa              >  2:111301/1‑100 (MQ=255)

                                                                                        aaaaGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATtaat             >  2:1016815/1‑100 (MQ=255)

                                                                                          aaGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATtaataa           <  1:1550965/100‑1 (MQ=255)

                                                                                            gTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACa         <  1:3183575/100‑1 (MQ=255)

                                                                                              tgtgGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACAtt       >  2:1490192/1‑100 (MQ=255)

                                                                                               gtgGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTa      <  2:488839/100‑1 (MQ=255)

                                                                                               gtgGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTa      >  1:3883360/1‑100 (MQ=255)

                                                                                                tgGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAg     <  1:3537775/100‑1 (MQ=255)

                                                                                                 ggCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAGa    <  1:3715939/100‑1 (MQ=255)

                                                                                                  gCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAAGACTTAATCAGATTAATAACATTAGat   <  1:723908/100‑1 (MQ=255)

                                                                                                   cTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAGata  >  2:3638162/1‑100 (MQ=255)

TCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAGATA  >  NC_000913/3651457‑3651655

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