BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,615,859	T→G	100%	T403P (ACC→CCC)	yjjI ←	DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,615,859	0	T	G	100.0%	347.8 / ‑6.7	101	T403P (ACC→CCC)	yjjI	DUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand): ref base T (0/2); new base G (54/45); total (54/47)

Only 100 of 103 total aligned reads displayed.

CTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCG  >  NC_000913/4615764‑4615952

cTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATc                                                                                           >  1:2089383/1‑100 (MQ=255)

     ggTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGggatcggatc                                                                                      >  2:1659294/1‑100 (MQ=255)

     ggTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAAGAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGggatcggatc                                                                                      <  1:1143198/100‑1 (MQ=255)

      gTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCt                                                                                     >  1:1955851/1‑100 (MQ=255)

       ttACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTg                                                                                    <  2:1648844/100‑1 (MQ=255)

       ttACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTg                                                                                    <  2:1153815/100‑1 (MQ=255)

        tACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTgg                                                                                   >  1:2090298/1‑100 (MQ=255)

         aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                  <  2:534981/100‑1 (MQ=255)

         aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                  >  1:778751/1‑100 (MQ=255)

         aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                  >  1:957167/1‑100 (MQ=255)

          ngTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                 >  2:706398/2‑100 (MQ=255)

          cGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                 <  1:1938280/100‑1 (MQ=255)

          cGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                 >  2:2399476/1‑100 (MQ=255)

          cGTTTGATGGTTTCGTCGAGCGACAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                 >  2:1038042/1‑100 (MQ=255)

            tttGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCa                                                                               >  1:2277289/1‑100 (MQ=255)

              tGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATc                                                                             <  1:697117/100‑1 (MQ=255)

               gATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCg                                                                            <  1:821401/100‑1 (MQ=255)

               gATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCg                                                                            <  2:2072763/100‑1 (MQ=255)

                  ggTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCa                                                                         >  1:1274876/1‑100 (MQ=255)

                   gTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGAGGATCGGATCTGGCTCATCGCCAt                                                                        >  1:641952/1‑100 (MQ=255)

                       cgtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACgg                                                                    >  2:1335861/1‑100 (MQ=255)

                       cgtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACgg                                                                    <  1:591642/100‑1 (MQ=255)

                       cgtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACgg                                                                    >  2:289827/1‑100 (MQ=255)

                        gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                   <  1:2031377/100‑1 (MQ=255)

                        gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                   >  2:874658/1‑100 (MQ=255)

                         tcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                   <  1:1300505/99‑1 (MQ=255)

                           gAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCaaa                                                                >  1:475374/1‑100 (MQ=255)

                            aGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAc                                                               >  1:746396/1‑100 (MQ=255)

                            aGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAc                                                               <  1:2175530/100‑1 (MQ=255)

                              cGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcg                                                              <  1:3680/99‑1 (MQ=255)

                               gTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcg                                                            <  1:1463237/100‑1 (MQ=255)

                               gTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcg                                                            <  1:2315817/100‑1 (MQ=255)

                                tctgAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcg                                                            >  2:358004/4‑99 (MQ=255)

                                  aGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACgcgcgc                                                         <  1:1129911/100‑1 (MQ=255)

                                   gAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCAGGCTAATCGCCATACGGCAAAcgc                                                             >  2:2264412/1‑95 (MQ=255)

                                       gTCGCTGATGCCGGAATAATAATAAGAATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCGCATAGCCATTCGGCAAACGCGCGCCCg                                                      >  2:2093030/1‑98 (MQ=255)

                                         cGCTGATGCCGGAGTAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgt                                                  >  1:1610010/1‑100 (MQ=255)

                                           cTGATGTCGGAATAATAATAAGCATGATGCGGGGTGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATAGCCATACTGCAAACGCGCGCCCGGCgtgg                                                <  2:707885/100‑1 (MQ=255)

                                           cTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtgg                                                <  2:110545/100‑1 (MQ=255)

                                            tGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtggt                                               >  1:1781154/1‑100 (MQ=255)

                                             gATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtggtg                                              <  2:685586/100‑1 (MQ=255)

                                             gATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtggtg                                              >  1:1704054/1‑100 (MQ=255)

                                              aTGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGc                                             <  1:1935145/100‑1 (MQ=255)

                                               tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                            <  1:113975/100‑1 (MQ=255)

                                                gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                           >  1:1410832/1‑100 (MQ=255)

                                                gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                           <  2:50503/100‑1 (MQ=255)

                                                gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                           <  2:1783257/100‑1 (MQ=255)

                                                 ccGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGa                                          >  1:1148982/1‑100 (MQ=255)

                                                  cGGAATAATAATAAGCATGATGCGGGGCGGCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGat                                         <  2:1907965/100‑1 (MQ=255)

                                                  cGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGat                                         <  1:1174067/100‑1 (MQ=255)

                                                  cGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGat                                         <  2:1746445/100‑1 (MQ=255)

                                                  cGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGa                                          >  1:496557/1‑99 (MQ=255)

                                                   ggAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGata                                        >  2:316861/1‑100 (MQ=255)

                                                   ggAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGata                                        >  1:727540/1‑100 (MQ=255)

                                                   ggAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGata                                        <  1:2337110/100‑1 (MQ=255)

                                                       taataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                    >  2:32450/1‑100 (MQ=255)

                                                       taataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                    >  2:1539163/1‑100 (MQ=255)

                                                        aataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                    <  2:321297/99‑1 (MQ=255)

                                                        aataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGa                                   <  2:1264290/100‑1 (MQ=255)

                                                        aataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGa                                   >  1:1731570/1‑100 (MQ=255)

                                                          taataaGCATGATGCGGGGGGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGTCCGGCGTGGTGCAAATATCGGAAc                                 <  2:46865/100‑1 (MQ=255)

                                                          taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                 <  2:1579462/100‑1 (MQ=255)

                                                           aataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACt                                >  2:2147048/1‑100 (MQ=255)

                                                            ataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTAGTGCCGATATCGGAACTg                               >  1:154296/1‑100 (MQ=255)

                                                             taaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGa                              >  1:2064413/1‑100 (MQ=255)

                                                              aaGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGAt                             <  2:2171170/100‑1 (MQ=255)

                                                               aGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGAt                             >  2:1814119/1‑99 (MQ=255)

                                                                gtatgatgCGGGGCGAAAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGGCCGGCGTGGTGCCGATATCGGAACTGa                              <  1:835444/95‑1 (MQ=255)

                                                                gCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcc                           <  1:1911451/100‑1 (MQ=255)

                                                                gCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcc                           <  1:1439906/100‑1 (MQ=255)

                                                                  atgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcccc                         <  1:1118668/100‑1 (MQ=255)

                                                                  atgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcccc                         <  2:2064168/100‑1 (MQ=255)

                                                                   tgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcccct                        <  2:1741541/100‑2 (MQ=255)

                                                                      tgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGTAACTGATCCCCGACt                     >  2:1008484/1‑100 (MQ=255)

                                                                      tgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACt                     >  2:383495/1‑100 (MQ=255)

                                                                       gCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtg                    <  1:947267/100‑1 (MQ=255)

                                                                          gggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGc                 >  2:1710134/1‑100 (MQ=255)

                                                                           ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtg                    >  1:741034/1‑96 (MQ=255)

                                                                           ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                <  2:697117/100‑1 (MQ=255)

                                                                           ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                <  2:2150655/100‑1 (MQ=255)

                                                                           ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                <  1:1753419/100‑1 (MQ=255)

                                                                            gCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgt               >  1:52309/1‑100 (MQ=255)

                                                                            gCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgt               >  1:971378/1‑100 (MQ=255)

                                                                            gCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCAACTGTGCgt               >  1:1157753/1‑100 (MQ=255)

                                                                             cGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg              <  1:1067486/100‑1 (MQ=255)

                                                                              gACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt             <  2:2228253/100‑1 (MQ=255)

                                                                              gACAGTTTGCAGATGGGGGATCGGATCTGTCTCATCGCCATACGGGAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg              <  2:385750/99‑1 (MQ=255)

                                                                              gACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt             >  1:899863/1‑100 (MQ=255)

                                                                              gACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt             >  2:1639682/1‑100 (MQ=255)

                                                                              gACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt             <  2:778742/100‑1 (MQ=255)

                                                                                cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa           >  1:603326/1‑100 (MQ=255)

                                                                                cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa           >  1:1868590/1‑100 (MQ=255)

                                                                                 aGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAAc          <  2:2291359/100‑1 (MQ=255)

                                                                                  gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa         >  1:1892932/1‑100 (MQ=255)

                                                                                   tttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACAt        >  2:1993296/1‑100 (MQ=255)

                                                                                     tGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATgg      >  1:1584498/1‑100 (MQ=255)

                                                                                       ccgATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCACGGCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcg    >  2:21720/3‑100 (MQ=255)

                                                                                       cAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcg    >  2:1727555/1‑100 (MQ=255)

                                                                                         gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGcgcg  >  2:2039187/1‑100 (MQ=255)

                                                                                         gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACAAGGcgcg  >  2:2270786/1‑100 (MQ=255)

CTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCG  >  NC_000913/4615764‑4615952

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