BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,840,760	G→A	D432N (GAC→AAC)	ynjE →	molybdopterin synthase sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,840,760	0	G	A	100.0%	1241.3 / NA	293	D432N (GAC→AAC)	ynjE	molybdopterin synthase sulfurtransferase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (132/161); total (132/161)

Only 100 of 293 total aligned reads displayed.

GCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGC  >  NC_000913/1840662‑1840857

gcgcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAAc                                                                                                 >  2:3315409/1‑101 (MQ=255)

gcgcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAAc                                                                                                 >  1:3315408/1‑101 (MQ=255)

  gcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAg                                                                                               <  2:3314946/101‑1 (MQ=255)

  gcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAg                                                                                               >  1:3315413/1‑101 (MQ=255)

   ccATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGt                                                                                              <  2:3315016/101‑1 (MQ=255)

   ccATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGt                                                                                              <  2:3314896/101‑1 (MQ=255)

    cATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTa                                                                                             >  1:3315416/1‑101 (MQ=255)

    cATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTa                                                                                             <  1:3315065/101‑1 (MQ=255)

      tGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAg                                                                                               <  1:3314760/97‑1 (MQ=255)

       gggTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCa                                                                                          <  2:3314823/101‑1 (MQ=255)

       gggTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCa                                                                                          >  1:3315421/1‑101 (MQ=255)

       gggTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCa                                                                                          <  1:3315184/101‑1 (MQ=255)

        ggTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCaa                                                                                         <  2:3314941/101‑1 (MQ=255)

         gTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCaaa                                                                                        <  1:3314736/101‑1 (MQ=255)

           tGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATa                                                                                      <  2:3314783/101‑1 (MQ=255)

           tGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATa                                                                                      >  1:3315426/1‑101 (MQ=255)

           tGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATa                                                                                      >  2:3315425/1‑101 (MQ=255)

            ggAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATaa                                                                                     <  1:3314913/101‑1 (MQ=255)

            ggAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATaa                                                                                     <  1:3314843/101‑1 (MQ=255)

             gaagaaTGTTTCCGTGTATGACGTCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAAc                                                                                    >  1:3315428/1‑101 (MQ=255)

              aagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAg                                                                                               <  2:3315056/89‑1 (MQ=255)

              aagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACa                                                                                   >  2:3315431/1‑101 (MQ=255)

                gaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACAtt                                                                                 >  2:3315435/1‑101 (MQ=255)

                  aTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGc                                                                               >  1:3315439/1‑101 (MQ=255)

                      ttCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAAtt                                                                           >  1:3315442/1‑101 (MQ=255)

                        ccGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTAc                                                                         >  2:3315444/1‑101 (MQ=255)

                        ccGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTAc                                                                         <  2:3315093/101‑1 (MQ=255)

                         cGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACt                                                                        <  1:3314938/101‑1 (MQ=255)

                         cGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACt                                                                        <  2:3314879/101‑1 (MQ=255)

                          gtgtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTg                                                                       <  1:3315091/101‑1 (MQ=255)

                          gtgtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTg                                                                       >  1:3315450/1‑101 (MQ=255)

                          gtgtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTg                                                                       >  2:3315449/1‑101 (MQ=255)

                            gtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGAc                                                                     >  1:3315455/1‑101 (MQ=255)

                             tATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGa                                                                      >  2:3315457/1‑99 (MQ=255)

                              aTGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGc                                                                   <  2:3315050/101‑1 (MQ=255)

                              aTGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGc                                                                   >  2:3315459/1‑101 (MQ=255)

                               tGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCt                                                                  <  2:3314979/101‑1 (MQ=255)

                               tGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCt                                                                  >  1:3315463/1‑101 (MQ=255)

                                  cggcggCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGa                                                               >  2:3315465/1‑101 (MQ=255)

                                    gcggcTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGc                                                             <  1:3315010/101‑1 (MQ=255)

                                     cggcTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCg                                                            <  1:3315152/101‑1 (MQ=255)

                                     cggcTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCg                                                            >  1:3315468/1‑101 (MQ=255)

                                         tGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACtt                                                        >  2:3315472/1‑101 (MQ=255)

                                          ggTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTc                                                       <  1:3315255/101‑1 (MQ=255)

                                          ggTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTc                                                       <  1:3314834/101‑1 (MQ=255)

                                          ggTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTc                                                       <  2:3315154/101‑1 (MQ=255)

                                            tACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAg                                                     <  2:3315025/101‑1 (MQ=255)

                                            tACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAg                                                     <  2:3315127/101‑1 (MQ=255)

                                             aCGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                                                    >  1:3315475/1‑101 (MQ=255)

                                              cGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCg                                                   >  2:3315477/1‑101 (MQ=255)

                                               gAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGt                                                  <  2:3315068/101‑1 (MQ=255)

                                                aaTGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                 <  2:3315252/101‑1 (MQ=255)

                                                aaTGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                 >  1:3315480/1‑101 (MQ=255)

                                                  tGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAg                                               >  2:3315485/1‑101 (MQ=255)

                                                  tGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAg                                               <  2:3315036/101‑1 (MQ=255)

                                                   ggAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                                              <  2:3314851/101‑1 (MQ=255)

                                                    gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGAt                                             >  1:3315486/1‑101 (MQ=255)

                                                     agcagcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATa                                            <  2:3314868/101‑1 (MQ=255)

                                                      gcagcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAc                                           >  1:3315490/1‑101 (MQ=255)

                                                      gcagcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAc                                           <  1:3314801/101‑1 (MQ=255)

                                                        agcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAccc                                         >  2:3315493/1‑101 (MQ=255)

                                                        agcGATCAAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACt                                                                        <  2:3314959/70‑1 (MQ=255)

                                                          cGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCAc                                       >  1:3315496/1‑101 (MQ=255)

                                                          cGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCAc                                       >  2:3315498/1‑101 (MQ=255)

                                                           gATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCg                                                            >  2:3315503/1‑79 (MQ=255)

                                                           gATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACt                                      <  1:3315170/101‑1 (MQ=255)

                                                            aTCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTc                                     <  2:3314935/101‑1 (MQ=255)

                                                             tCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTcc                                    >  1:3315505/1‑101 (MQ=255)

                                                                aaaaaaTCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcaca                                 >  2:3315509/1‑101 (MQ=255)

                                                                 aaaaaTCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCacac                                <  1:3314999/101‑1 (MQ=255)

                                                                 aaaaaTCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCacac                                <  1:3314987/101‑1 (MQ=255)

                                                                 aaaaaTCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCacac                                <  2:3314997/101‑1 (MQ=255)

                                                                  aaaaTCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACAcc                               <  2:3314986/101‑1 (MQ=255)

                                                                  aaaaTCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACAcc                               >  1:3315511/1‑101 (MQ=255)

                                                                   aaaTCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACAccc                              <  2:3315101/101‑1 (MQ=255)

                                                                     aTCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgc                            >  1:3315514/1‑101 (MQ=255)

                                                                      tCCGGTAGCAACCGGTGAACGCGGCCCTAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgcg                           <  2:3315007/101‑1 (MQ=255)

                                                                      tCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                                                    <  1:3315105/76‑1 (MQ=255)

                                                                      tCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgcg                           <  1:3314923/101‑1 (MQ=255)

                                                                      tCCGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgcg                           >  2:3315515/1‑101 (MQ=255)

                                                                       ccGGTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGa                                                                      <  2:3314656/57‑1 (MQ=255)

                                                                         ggTAGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTaa                        >  2:3315519/1‑101 (MQ=255)

                                                                            aGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAgg                     <  2:3315079/101‑1 (MQ=255)

                                                                            aGCAACCGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAgg                     >  2:3315523/1‑101 (MQ=255)

                                                                                 ccGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcgtc                <  1:3315029/101‑1 (MQ=255)

                                                                                 ccGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcgtc                <  1:3314880/101‑1 (MQ=255)

                                                                                 ccGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcgtc                >  1:3315525/1‑101 (MQ=255)

                                                                                  cGGTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAg                      >  1:3315529/1‑94 (MQ=255)

                                                                                   ggTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCAt              <  2:3314965/101‑1 (MQ=255)

                                                                                    gTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtc                   <  1:3314857/95‑1 (MQ=255)

                                                                                    gTGAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATc             >  1:3315532/1‑101 (MQ=255)

                                                                                      gAACGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCa           >  2:3315533/1‑101 (MQ=255)

                                                                                       aaCGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAg          <  2:3315166/101‑1 (MQ=255)

                                                                                        aCGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGc         <  2:3315181/101‑1 (MQ=255)

                                                                                        aCGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGc         <  2:3315109/101‑1 (MQ=255)

                                                                                        aCGCGGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGc         <  1:3314947/101‑1 (MQ=255)

                                                                                          gcgGCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGcac       >  1:3315536/1‑101 (MQ=255)

                                                                                             gCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAAc    <  1:3315039/101‑1 (MQ=255)

                                                                                             gCCCGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAAc    >  2:3315540/1‑101 (MQ=255)

                                                                                              cccGAACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACg   <  1:3315128/101‑1 (MQ=255)

GCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGC  >  NC_000913/1840662‑1840857

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