BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,651,556	C→T	37.4%	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	37.4%	179.3 / 189.3	199	pseudogene (266/741 nt)	yhiS	pseudogene
Reads supporting (aligned to +/- strand): ref base C (67/57); new base T (34/40); total (101/98)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 203 total aligned reads displayed.

TCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAGAT  >  NC_000913/3651457‑3651654

tctcTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACCTATACAAATGCAGACTACTCAAAAGTTGTGGc                                                                                                    <  1:2062223/100‑1 (MQ=255)

 ctctGATCAATATGATGGAGCGCAACACCAACTTCAACAAATTAATGCGTGCAACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGc                                                                                                    <  1:624793/99‑1 (MQ=255)

    tGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCaa                                                                                                <  1:935033/100‑1 (MQ=255)

     gATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAAc                                                                                               >  2:2118003/1‑100 (MQ=255)

       tCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACtt                                                                                             >  2:678190/1‑100 (MQ=255)

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

          atatGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTa                                                                                          >  1:1729733/1‑100 (MQ=255)

           tatGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTaa                                                                                         <  2:2296854/100‑1 (MQ=255)

            atGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAg                                                                                        <  2:366476/100‑1 (MQ=255)

            atGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAg                                                                                        <  1:3391891/100‑1 (MQ=255)

             tGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGt                                                                                       >  2:1460695/1‑100 (MQ=255)

             tGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGt                                                                                       >  2:875797/1‑100 (MQ=255)

               aTTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCt                                                                                     >  1:1000603/1‑100 (MQ=255)

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

                tttAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCtt                                                                                    >  1:813227/1‑100 (MQ=255)

                  tAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCttat                                                                                  >  2:2279439/1‑100 (MQ=255)

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

                    gcgcAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCttatta                                                                                >  2:3142267/1‑100 (MQ=255)

                    gcgcAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCttatta                                                                                <  2:813229/100‑1 (MQ=255)

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

                     cgaaaCACCAGCTTAAAGAACGTACTCAGTGCCAGCCGAAACCTGCATATACAGTTGCAGCCTACTGTCAAGTTGTGGGTCAACTTGTAAGTCttaagag                                                                               <  1:3719612/97‑5 (MQ=255)

                      gcAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACt                                                                              <  2:1505570/100‑1 (MQ=255)

                      gcAACACCAACTTCAACAAATTAATGCGTGCCACAAGAACACTACATATACAAATGCAGACTACTCAAAAGTTCTAGCTCAACTTGTAAGTCTGATTAc                                                                               >  1:3620831/1‑99 (MQ=255)

                        aaCACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTaa                                                                            <  2:3978419/100‑1 (MQ=255)

                         acacCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAAc                                                                           >  2:3055858/1‑100 (MQ=255)

                          cacCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACa                                                                          <  2:110686/100‑1 (MQ=255)

                          cacCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACa                                                                          <  1:2022657/100‑1 (MQ=255)

                           acCAACTTCAACAAATTAATGCGTGCCACACAAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACAt                                                                         <  2:3650908/100‑1 (MQ=255)

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

                             cAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTg                                                                       >  1:1673134/1‑100 (MQ=255)

                                 ttCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAAc                                                                   <  2:1865615/100‑1 (MQ=255)

                                   caacaaATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGa                                                                 <  2:1651287/100‑1 (MQ=255)

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

                                      caaATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATaa                                                              >  1:1202161/1‑100 (MQ=255)

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

                                       aaaTTAATGCGTGCAACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTAAAACGATaa                                                              >  1:732313/1‑99 (MQ=255)

                                         aTTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGtt                                                           <  1:1459007/100‑1 (MQ=255)

                                          ttAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTc                                                          <  2:2080300/100‑1 (MQ=255)

                                          ttAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTc                                                          <  2:132218/100‑1 (MQ=255)

                                           tAATGCGTGCCACACGAACACTACATATACAAATGAAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGtt                                                           >  2:1748893/1‑98 (MQ=255)

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

                                              tGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGaca                                                      >  1:3708969/1‑100 (MQ=255)

                                              tGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGaca                                                      <  1:1533746/100‑1 (MQ=255)

                                              tGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGaca                                                      >  2:907091/1‑100 (MQ=255)

                                              tGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGaca                                                      <  2:2785570/100‑1 (MQ=255)

                                                cGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAATTCTTATTACTAACATTGAAACGATAAGTTCGAcaca                                                    <  2:3008198/100‑1 (MQ=255)

                                                    ccACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACtt                                                >  2:2148366/1‑100 (MQ=255)

                                                    ccACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACtt                                                >  2:2778919/1‑100 (MQ=255)

                                                      acacGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTAc                                              <  1:2721943/100‑1 (MQ=255)

                                                       cacGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTACGTATTATAACTAACATTGAAACGATACGTTCGacgc                                                     >  1:3330691/1‑90 (MQ=255)

                                                        acGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTa                                               <  2:4242418/97‑1 (MQ=255)

                                                         cGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCa                                           >  1:2980777/1‑100 (MQ=255)

                                                             cacTACATATACAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACa                                        >  2:442624/1‑99 (MQ=255)

                                                             cacTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAg                                       <  1:2223911/100‑1 (MQ=255)

                                                               cTACATATACAAATGTAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGac                                     <  1:3286223/100‑1 (MQ=255)

                                                               atACATATACAAAAGCAGACTACTAACAAGTTGTGGTTCAACTTGTAAGTATTATTACTAACATTAAAACGATAAGTTAGACACAACATACGCAACa                                        >  1:179291/2‑97 (MQ=255)

                                                                 aCATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGacac                                   >  1:3081055/1‑100 (MQ=255)

                                                                 aCATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGacac                                   >  2:3482231/1‑100 (MQ=255)

                                                                    tataCAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaat                                <  1:2343913/100‑1 (MQ=255)

                                                                     ataCAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatc                               <  2:2022262/100‑1 (MQ=255)

                                                                     ataCAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatc                               <  1:1065245/100‑1 (MQ=255)

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

                                                                      taCAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatca                              <  1:2778908/100‑1 (MQ=255)

                                                                      taCAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatca                              >  1:1643549/1‑100 (MQ=255)

                                                                      taCAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatca                              >  2:2870116/1‑100 (MQ=255)

                                                                       aCAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatcaa                             <  2:2488296/100‑1 (MQ=255)

                                                                        cAAATGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACAcaatcaat                            <  1:1363818/100‑1 (MQ=255)

                                                                          aaTGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATAc                          <  1:4149/100‑1 (MQ=255)

                                                                          aaTGCAGACTACTCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATAc                          >  2:959761/1‑100 (MQ=255)

                                                                          aaTGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATAc                          <  1:3482218/100‑1 (MQ=255)

                                                                          aaTGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATAc                          <  1:773906/100‑1 (MQ=255)

                                                                             gCAGACTACTCAAAAGTTGTGGCTCAGTGTATATCTCGGATTTCTCAGATTGATTCAATAAAATAGACACAACTTACGCAACAGACACAATCAATACTTa                       >  1:712242/1‑100 (MQ=255)

                                                                              cAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACTCAATCAATACTTaa                      <  2:2542632/100‑1 (MQ=255)

                                                                                 actactCAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCa                   <  2:2689027/100‑1 (MQ=255)

                                                                                  ctactCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAg                  <  2:96207/100‑1 (MQ=255)

                                                                                   tactCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTACGACACAACTTACGCAACAGACACAATCAATACTTAATCAg                  <  1:1248085/99‑1 (MQ=255)

                                                                                     cttacaaGTTGAGGCACAACTTGTAAGTCTTATTACTAAAATTGAAAGGATAAGTTCGACACAACTTACGAAACAAACACAATCAATACTTAATCAGa                 >  1:3123801/6‑98 (MQ=255)

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

                                                                                     ctCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGAtt               <  1:1890536/100‑1 (MQ=255)

                                                                                     ctCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACCAACCTTGAAACGATAAGTTAGACACAACTTACGCAACAGACACCATTCATACTTTATCAGa                 >  1:952715/1‑98 (MQ=255)

                                                                                       cAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGAtt               >  1:3886842/1‑98 (MQ=255)

                                                                                       cAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATtaa             >  1:3428193/1‑100 (MQ=255)

                                                                                       cAAAAGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATtaa             >  2:2360016/1‑100 (MQ=255)

                                                                                       cAAAAGTTGTGGCTCAACTTGTAAGTCTTATTCCTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGAt                >  1:2584928/1‑97 (MQ=255)

                                                                                       cAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATtaa             >  1:497018/1‑100 (MQ=255)

                                                                                       cAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATta              >  1:2667467/1‑99 (MQ=255)

                                                                                        aaaaGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATtaat            >  2:1832853/1‑100 (MQ=255)

                                                                                         aaaGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATtaata           <  1:4191575/100‑1 (MQ=255)

                                                                                         aaaGTTGTGGTTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATtaata           <  1:3440360/100‑1 (MQ=255)

                                                                                         aaaGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATtaata           <  1:4061727/100‑1 (MQ=255)

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

                                                                                            gTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACa        >  2:514954/1‑100 (MQ=255)

                                                                                            gTTGTGGCTCAACGTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACa        >  2:2903957/1‑100 (MQ=255)

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

                                                                                               gtgGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTa     <  2:402603/100‑1 (MQ=255)

                                                                                               gtgGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTa     <  1:2127228/100‑1 (MQ=255)

                                                                                               gtgGCTCAACGTGTAAGTCTTATGACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACATACACAATCAATACTTAATCAGATTAATAACATTa     <  2:1911443/100‑1 (MQ=255)

                                                                                                  gCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAGat  >  2:3360741/1‑100 (MQ=255)

                                                                                                  gCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAGat  <  2:147692/100‑1 (MQ=255)

                                                                                                  gCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAGat  >  1:1155056/1‑100 (MQ=255)

TCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAGAT  >  NC_000913/3651457‑3651654

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