BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	239,326	C→T	pseudogene (137/189 nt)	yafF →	putative uncharacterized protein YafF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	239,326	0	C	T	100.0%	1327.5 / NA	387	pseudogene (137/189 nt)	yafF	putative uncharacterized protein YafF
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (196/191); total (196/191)

Only 100 of 389 total aligned reads displayed.

AGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACC  >  NC_000913/239228‑239423

agcagAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCa                                                                                                 <  2:467102/101‑1 (MQ=17)

 gcagAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCAt                                                                                                >  1:468303/1‑101 (MQ=17)

  cagAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATg                                                                                               >  2:468305/1‑101 (MQ=17)

   agAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATgg                                                                                              >  1:468308/1‑101 (MQ=18)

     aaTTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGAc                                                                                            <  2:467630/101‑1 (MQ=18)

     aaTTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGAc                                                                                            >  2:468312/1‑101 (MQ=18)

      attattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACa                                                                                           <  1:467769/101‑1 (MQ=18)

       ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAg                                                                                          <  1:467570/101‑1 (MQ=18)

       ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAg                                                                                          >  2:468317/1‑101 (MQ=18)

         attTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGaa                                                                                        >  1:468321/1‑101 (MQ=18)

         attTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGaa                                                                                        >  2:468324/1‑101 (MQ=18)

          ttttCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGaaa                                                                                       <  1:467647/101‑1 (MQ=18)

           tttCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATg                                                                                               >  2:468330/1‑92 (MQ=18)

            ttCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACt                                                                                     >  1:468331/1‑101 (MQ=17)

            ttCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACt                                                                                     <  1:467599/101‑1 (MQ=18)

             tCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACa                                                                                           >  1:468336/1‑94 (MQ=18)

                 ggATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTg                                                                                >  2:468344/1‑101 (MQ=17)

                 ggATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTg                                                                                >  1:468345/1‑101 (MQ=18)

                   aTACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGc                                                                              <  1:467704/101‑1 (MQ=18)

                    tACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCg                                                                             >  1:468352/1‑101 (MQ=18)

                     aCGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGt                                                                            >  2:468353/1‑101 (MQ=18)

                      cGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTc                                                                           <  2:467699/101‑1 (MQ=17)

                         cacaTTGCTATTAATATTTTGACGAATGGGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGt                                                                        >  2:468364/1‑101 (MQ=255)

                         cacaTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGt                                                                        <  1:467816/101‑1 (MQ=17)

                         cacaTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGt                                                                        <  2:467989/101‑1 (MQ=17)

                         cacaTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGt                                                                        >  1:468362/1‑101 (MQ=18)

                          acaTTGCTATTAATATTTTGACGAATGAGAAGGTATTGAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCa                                                                          <  1:467726/98‑1 (MQ=16)

                          acaTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTc                                                                       <  1:467503/101‑1 (MQ=17)

                           caTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTcc                                                                      <  2:467831/101‑1 (MQ=17)

                           caTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTcc                                                                      <  2:468015/101‑1 (MQ=17)

                            aTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCt                                                                     <  2:467939/101‑1 (MQ=17)

                             ttGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCtt                                                                    <  2:467556/101‑1 (MQ=17)

                             ttGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCtt                                                                    >  1:468369/1‑101 (MQ=18)

                              tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCtt                                                                    <  2:467880/100‑1 (MQ=18)

                              tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                                                   >  2:468374/1‑101 (MQ=25)

                              tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                                                   <  2:467910/101‑1 (MQ=18)

                               gCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGc                                                                  <  1:467369/101‑1 (MQ=25)

                                cTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCg                                                                 >  2:468378/1‑101 (MQ=25)

                                 tATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTCCCTGGCGTCAGTCCTTGCgg                                                                <  1:467793/101‑1 (MQ=255)

                                 tATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                >  1:468383/1‑101 (MQ=25)

                                 tATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                <  2:467806/101‑1 (MQ=25)

                                 tATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                >  2:468382/1‑101 (MQ=38)

                                  aTTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCggg                                                               <  2:467999/101‑1 (MQ=25)

                                   ttAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgggg                                                              <  2:467742/101‑1 (MQ=25)

                                    tAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGa                                                             <  2:467895/101‑1 (MQ=25)

                                     aaTATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAg                                                            >  1:468391/1‑101 (MQ=38)

                                      atatTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGc                                                           >  2:468393/1‑101 (MQ=25)

                                        atTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCgg                                                         <  1:468126/101‑1 (MQ=25)

                                         ttttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGt                                                                            >  2:468400/1‑81 (MQ=18)

                                         ttttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCggg                                                        >  1:468399/1‑101 (MQ=38)

                                          tttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGc                                                       >  1:468404/1‑101 (MQ=25)

                                           ttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                                                      <  1:467690/101‑1 (MQ=255)

                                           ttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                                                      <  2:467610/101‑1 (MQ=25)

                                             gACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCttt                                                    <  2:467516/101‑1 (MQ=38)

                                              actaaTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTc                                                   >  1:468411/4‑101 (MQ=255)

                                              aCGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTc                                                   <  2:468248/101‑1 (MQ=38)

                                              aCGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTc                                                   >  2:468410/1‑101 (MQ=25)

                                                gAATGAGAAAGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCg                                                                 >  1:468418/1‑85 (MQ=255)

                                                 aaTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAg                                                            <  2:467667/89‑1 (MQ=255)

                                                 aaTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgta                                                <  1:467848/101‑1 (MQ=38)

                                                 aaTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgta                                                >  2:468419/1‑101 (MQ=25)

                                                    gagaAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagta                                             <  1:467875/101‑1 (MQ=255)

                                                    gagaAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagta                                             <  2:467926/101‑1 (MQ=255)

                                                     agaAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTaa                                            <  1:467885/101‑1 (MQ=255)

                                                      gaAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAAt                                           >  1:468424/1‑101 (MQ=255)

                                                        aGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAc                                         >  2:468428/1‑101 (MQ=255)

                                                          gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAccc                                       >  2:468433/1‑101 (MQ=255)

                                                           tATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAcccc                                      >  1:468436/1‑101 (MQ=255)

                                                            aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCg                                     <  2:467760/101‑1 (MQ=255)

                                                             ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGa                                    >  1:468442/1‑101 (MQ=255)

                                                             ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGa                                    >  2:468443/1‑101 (MQ=255)

                                                               cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGAct                                  <  1:467964/101‑1 (MQ=255)

                                                                aaGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGActc                                 >  1:468451/1‑101 (MQ=255)

                                                                aaGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGActc                                 >  2:468450/1‑101 (MQ=255)

                                                                   gCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTcc                              <  2:468086/101‑1 (MQ=255)

                                                                   gCAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTcc                              <  1:467859/101‑1 (MQ=255)

                                                                    cAGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTccc                             >  2:468458/1‑101 (MQ=255)

                                                                     aGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTcccc                            >  1:468459/1‑101 (MQ=255)

                                                                         ttAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTcc                        >  2:468466/1‑101 (MQ=255)

                                                                         ttAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTcc                        >  1:468467/1‑101 (MQ=255)

                                                                           aaGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCtt                      <  2:468121/101‑1 (MQ=255)

                                                                            aGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTa                     >  2:468473/1‑101 (MQ=255)

                                                                            aGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTa                     <  1:467919/101‑1 (MQ=255)

                                                                              aCGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTaaa                   >  1:468475/1‑101 (MQ=255)

                                                                               cGTAAGATGCGAAAAGCAGTAATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTaaa                   <  1:467901/100‑1 (MQ=255)

                                                                                gTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAaca                 >  1:468480/1‑101 (MQ=255)

                                                                                gTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAaca                 <  1:468001/101‑1 (MQ=255)

                                                                                 tAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAacac                <  1:468225/101‑1 (MQ=255)

                                                                                  aaGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAAcaca               <  2:467852/101‑1 (MQ=255)

                                                                                  aaGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAAcaca               >  2:468484/1‑101 (MQ=255)

                                                                                  aaGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAAcaca               >  1:468486/1‑101 (MQ=255)

                                                                                   aGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACaa              <  2:467591/101‑1 (MQ=255)

                                                                                     aTGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAAcc            <  1:468029/101‑1 (MQ=255)

                                                                                     aTGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAAcc            <  1:468052/101‑1 (MQ=255)

                                                                                       gCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAAcccc          >  1:468493/1‑101 (MQ=255)

                                                                                         gAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCa        >  2:468495/1‑101 (MQ=255)

                                                                                         gAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCa        <  1:467829/101‑1 (MQ=255)

                                                                                         gAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTACCCGt                          <  1:467784/83‑1 (MQ=255)

                                                                                          aaaaGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGt                          <  2:467785/82‑1 (MQ=255)

                                                                                          aaaaGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCAc       >  1:468500/1‑101 (MQ=255)

AGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACC  >  NC_000913/239228‑239423

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`