BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,887,352	T→A	100%	intergenic (‑133/+226)	ygcB ← / ← cysH	Cascade complex anti‑viral R‑loop helicase‑annealase Cas3/phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,887,352	0	T	A	100.0%	618.3 / ‑6.7	181	intergenic (‑133/+226)	ygcB/cysH	Cascade complex anti‑viral R‑loop helicase‑annealase Cas3/phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent
Reads supporting (aligned to +/- strand): ref base T (0/0); major base A (102/78); minor base C (1/0); total (103/78)

Only 100 of 181 total aligned reads displayed.

TATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTC  >  NC_000913/2887255‑2887447

tatCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGCGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGt                                                                                               >  2:2747113/1‑100 (MQ=255)

   cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAAAATTCGAATAGTTTg                                                                                            >  1:4036268/1‑100 (MQ=255)

     gCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTGTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTggg                                                                                          >  2:1145950/1‑100 (MQ=255)

      cGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGt                                                                                         <  2:4152788/100‑1 (MQ=255)

        ggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTc                                                                                       >  2:420512/1‑100 (MQ=255)

        ggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTc                                                                                       <  2:4184247/100‑1 (MQ=255)

          cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCg                                                                                      >  1:3313155/1‑99 (MQ=255)

          cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGa                                                                                     >  2:1422505/1‑100 (MQ=255)

          cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAACAAGAGACTAACATTCGAATAGTTTGGGTTCGa                                                                                     >  1:64695/1‑100 (MQ=255)

             tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACg                                                                                  <  1:71875/100‑1 (MQ=255)

               tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCt                                                                                <  1:869765/100‑1 (MQ=255)

                ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTg                                                                               <  1:612108/100‑1 (MQ=255)

                ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTg                                                                               >  1:821764/1‑100 (MQ=255)

                ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTg                                                                               <  2:3375382/100‑1 (MQ=255)

                ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTg                                                                               <  2:1117196/100‑1 (MQ=255)

                 tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                              <  2:1621017/100‑1 (MQ=255)

                  cGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgga                                                                             >  2:25101/1‑99 (MQ=255)

                   gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            <  2:1081985/100‑1 (MQ=255)

                   gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            <  2:1208351/100‑1 (MQ=255)

                   gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            <  2:3402711/100‑1 (MQ=255)

                   gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            <  2:1574073/100‑1 (MQ=255)

                   gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            >  2:2401091/1‑100 (MQ=255)

                      gctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCa                                                                         >  1:3164733/1‑100 (MQ=255)

                           tCATAAATATCTAATTACAGTTATATACACATGCCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                              <  1:1483428/90‑1 (MQ=255)

                           tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATAGTTGGGGTTCGAACGCTGGCCTCAGGtt                                                                     >  1:345352/1‑99 (MQ=255)

                           tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGtt                                                                     <  1:3870425/99‑1 (MQ=255)

                           tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTg                                                                    <  2:1119333/100‑1 (MQ=255)

                           tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTg                                                                    <  2:232487/100‑1 (MQ=255)

                            cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGa                                                                   <  2:1364147/100‑1 (MQ=255)

                            cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGa                                                                   <  2:4105182/100‑1 (MQ=255)

                            cATAAATATCGAATTACAGTGATATACACAGGTATTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGACTCAGGTTGa                                                                   >  1:2805205/1‑100 (MQ=255)

                             aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                  <  1:2045654/100‑1 (MQ=255)

                                   atCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAat                                                            <  2:2175798/100‑1 (MQ=255)

                                    tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCTAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAata                                                           >  1:2990534/1‑100 (MQ=255)

                                     cGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAatat                                                          >  2:3679340/1‑100 (MQ=255)

                                      gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATc                                                         >  1:751196/1‑100 (MQ=255)

                                      gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATc                                                         <  1:1683877/100‑1 (MQ=255)

                                         ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGcn                                                      <  1:1572374/100‑2 (MQ=255)

                                         ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGcc                                                      <  1:1587761/100‑1 (MQ=255)

                                           aCAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCCGGTTGATAGAAATATCGCCTg                                                    >  2:3524407/1‑100 (MQ=255)

                                            cAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                   <  2:2678834/100‑1 (MQ=255)

                                              gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 >  1:2741937/1‑100 (MQ=255)

                                              gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 <  2:292094/100‑1 (MQ=255)

                                              gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 <  2:293363/100‑1 (MQ=255)

                                              gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 >  2:3136225/1‑100 (MQ=255)

                                              gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGATGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 >  1:3013770/1‑100 (MQ=255)

                                              gTGATATACACATGTCTTGCCTCCCTGCAAAGAACAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 <  1:2247849/100‑1 (MQ=255)

                                               tGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGc                                                <  2:3206887/100‑1 (MQ=255)

                                                 atatACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCtt                                              <  2:2842001/100‑1 (MQ=255)

                                                  tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                             >  2:705981/1‑100 (MQ=255)

                                                    taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           >  1:449730/1‑100 (MQ=255)

                                                    taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           <  2:2152019/100‑1 (MQ=255)

                                                     acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                          >  2:2972841/1‑100 (MQ=255)

                                                     acacATGTCTTGCCTCACTGCAAAGATGAGAGTAACATTCGACTAGTTTGAAATAAAACCCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           >  1:2124863/1‑99 (MQ=255)

                                                      cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGACTAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                         >  2:1728180/1‑100 (MQ=255)

                                                      cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                         <  1:1983753/100‑1 (MQ=255)

                                                      cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                         >  2:2635172/1‑100 (MQ=255)

                                                        caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                       >  1:1148906/1‑100 (MQ=255)

                                                        caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCGCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                                        >  2:2999206/1‑99 (MQ=255)

                                                            tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                   >  2:3042766/1‑100 (MQ=255)

                                                            tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                   >  2:1334736/1‑100 (MQ=255)

                                                             cTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGTGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTgg                                  >  1:1760213/1‑100 (MQ=255)

                                                              ttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGa                                 >  2:1829772/1‑100 (MQ=255)

                                                              ttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGa                                 >  2:2052917/1‑100 (MQ=255)

                                                               tGCCTCCCTGCAAAGAGGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa                                >  1:1839835/1‑100 (MQ=255)

                                                                gCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc                               <  2:4250036/100‑1 (MQ=255)

                                                                 ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              >  1:3781390/1‑100 (MQ=255)

                                                                 ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              >  1:3632019/1‑100 (MQ=255)

                                                                 ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              >  2:674856/1‑100 (MQ=255)

                                                                  cTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt                             >  2:540335/1‑100 (MQ=255)

                                                                   tCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc                            >  1:3302792/1‑100 (MQ=255)

                                                                    cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGCACCTcg                           <  1:1263292/100‑1 (MQ=255)

                                                                    cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg                           <  1:2151207/100‑1 (MQ=255)

                                                                    cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg                           >  1:2377403/1‑100 (MQ=255)

                                                                     ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                          >  1:561687/1‑100 (MQ=255)

                                                                      cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         <  2:537204/100‑1 (MQ=255)

                                                                      cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         <  2:3622366/100‑1 (MQ=255)

                                                                        gCAAAGAAGAGACGAACATTAGAATAGTTTGGGTTCGATCGGTGGACTCAGGTTGATAGAAATATCGCGTGGGGCTTTTGTGCATCTGGGACCTggcg                         >  1:2647544/1‑98 (MQ=255)

                                                                         cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt                      >  2:2731300/1‑100 (MQ=255)

                                                                            agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                   <  2:2221430/100‑1 (MQ=255)

                                                                            agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                   <  1:433625/100‑1 (MQ=255)

                                                                             gaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCtt                  >  2:3071200/1‑100 (MQ=255)

                                                                             gaagaGACTAACATTCGAATAGTTTGGGGTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATATAGAAACTCACGAATGCtt                  >  1:1500968/1‑100 (MQ=255)

                                                                                 agaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACg              <  2:437187/100‑1 (MQ=255)

                                                                                   aCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc            <  1:1099891/100‑1 (MQ=255)

                                                                                    cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa           >  2:1633709/1‑100 (MQ=255)

                                                                                    cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa           >  2:1010168/1‑100 (MQ=255)

                                                                                    cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa           >  1:2193680/1‑100 (MQ=255)

                                                                                     tAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa           >  2:1302923/1‑99 (MQ=255)

                                                                                      aaCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGa         <  2:3913091/100‑1 (MQ=255)

                                                                                      aaCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGa         >  2:154782/1‑100 (MQ=255)

                                                                                        cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa       >  2:3963723/1‑100 (MQ=255)

                                                                                        cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa       >  1:3079682/1‑100 (MQ=255)

                                                                                         aTTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTGTTGTCCATCTGGAACCTCGCGAATGCTTAACgccagcca       >  1:100434/1‑99 (MQ=255)

                                                                                          ttCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc     >  2:793250/1‑100 (MQ=255)

                                                                                           tCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGcc    >  1:3620008/1‑100 (MQ=255)

                                                                                           tCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGAGCTTTTGTCCATCTGGAACCTCGCGAATGCATAACGCCAGACAGcc    >  1:2496869/1‑100 (MQ=255)

                                                                                            cGAATCGATTGGGTTCGAACGCTGGCCTCAGGTTGATAGAACTATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTCACGCCAGACAGcc    >  1:135404/1‑99 (MQ=255)

                                                                                            cGAATAGTTTTGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCt   >  2:597773/1‑100 (MQ=255)

                                                                                             gAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc  >  1:2968434/1‑100 (MQ=255)

TATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTC  >  NC_000913/2887255‑2887447

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