BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,127,778	C→T	G239G (GGC→GGT)	mviN →	putative virulence factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,127,778	0	C	T	100.0%	528.0 / NA	148	G239G (GGC→GGT)	mviN	putative virulence factor
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (57/91); total (57/91)

Only 100 of 152 total aligned reads displayed.

GAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCGA  >  NC_000913/1127688‑1127867

gaagaTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGgtgt                                                                                         >  1:89130/1‑93 (MQ=255)

   gaTCGGCATGCTGTTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGtctc                                                                                      <  2:4247610/93‑1 (MQ=255)

   gaTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGtctc                                                                                      >  2:1591786/1‑93 (MQ=255)

   gaTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGtctc                                                                                      <  1:3642323/93‑1 (MQ=255)

   gaTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGtctc                                                                                      <  1:2823435/93‑1 (MQ=255)

           tGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGcc                                                                              >  2:4678590/1‑93 (MQ=255)

           tGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGcc                                                                              <  2:2210310/93‑1 (MQ=255)

           tGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGcc                                                                              >  2:1562825/1‑93 (MQ=255)

           tGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGcc                                                                              >  1:852874/1‑93 (MQ=255)

              tGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGTCAGa                                                                           <  2:4519331/93‑1 (MQ=255)

              tGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGa                                                                           >  2:3176299/1‑93 (MQ=255)

              tGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGa                                                                           >  1:3345327/1‑93 (MQ=255)

               ggTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAt                                                                          <  1:1276209/93‑1 (MQ=255)

                gTCCTGCCGCGCATTATCTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAtc                                                                         <  2:4446034/93‑1 (MQ=255)

                gTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAtc                                                                         <  1:1129749/93‑1 (MQ=255)

                gTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAtc                                                                         >  1:996798/1‑93 (MQ=255)

                gTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAtc                                                                         >  2:4641440/1‑93 (MQ=255)

                 tCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAtct                                                                        >  1:539631/1‑93 (MQ=255)

                 tCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAtct                                                                        <  2:1820941/93‑1 (MQ=255)

                   cTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTcc                                                                      <  2:3952922/93‑1 (MQ=255)

                   cTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTcc                                                                      <  2:317422/93‑1 (MQ=255)

                   cTGCCGCGCATTAACTTCCACGATGCCGGAGCAAGGCGCGTGGTGAAACAGATGGGTCCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTcc                                                                      <  2:1151394/93‑1 (MQ=255)

                    tGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCt                                                                     <  1:4499767/93‑1 (MQ=255)

                    tGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCt                                                                     <  2:2945799/93‑1 (MQ=255)

                     gCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGCCCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCtt                                                                    <  2:71894/93‑1 (MQ=255)

                       cgcgCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTaa                                                                  >  1:231914/1‑93 (MQ=255)

                       cgcgCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTaa                                                                  <  1:4512320/93‑1 (MQ=255)

                         cgcATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatc                                                                <  2:2707066/93‑1 (MQ=255)

                          gcATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatca                                                               >  2:3701543/1‑93 (MQ=255)

                           cATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatcat                                                              >  1:4624703/1‑93 (MQ=255)

                           cATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatcat                                                              >  2:1155810/1‑93 (MQ=255)

                            aTTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatcatc                                                             <  2:3431317/93‑1 (MQ=255)

                            aTTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatcatc                                                             >  2:27060/1‑93 (MQ=255)

                            aTTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatcatc                                                             >  1:4654886/1‑93 (MQ=255)

                             ttAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTtcatcatca                                                            >  2:621988/1‑93 (MQ=255)

                             ttAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAAtcatca                                                            >  1:327908/1‑93 (MQ=255)

                             ttAACGTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAAtcatca                                                            <  2:2425476/93‑1 (MQ=255)

                               aaCTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAac                                                          >  2:3388552/1‑93 (MQ=255)

                                 cTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAacac                                                        <  1:1810024/93‑1 (MQ=255)

                                    ccACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAt                                                     <  2:2257517/93‑1 (MQ=255)

                                    ccACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAt                                                     >  2:1883599/1‑93 (MQ=255)

                                    ccACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAt                                                     <  2:3206229/93‑1 (MQ=255)

                                     cACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAtt                                                    <  2:2633768/93‑1 (MQ=255)

                                      aCGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttt                                                   >  2:3431284/1‑93 (MQ=255)

                                      aCGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttt                                                   <  1:1964293/93‑1 (MQ=255)

                                      aCGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttt                                                   <  1:1552995/93‑1 (MQ=255)

                                        gATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttttt                                                 <  2:3279893/93‑1 (MQ=255)

                                        gATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttttt                                                 >  1:1521437/1‑93 (MQ=255)

                                         aTGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTg                                                >  1:3987627/1‑93 (MQ=255)

                                         aTGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTg                                                <  1:3805459/93‑1 (MQ=255)

                                          tGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGc                                               <  2:4491608/93‑1 (MQ=255)

                                           gCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGcc                                              <  2:2396504/93‑1 (MQ=255)

                                           gCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGcc                                              <  2:4536264/93‑1 (MQ=255)

                                            ccGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCt                                             <  1:992706/93‑1 (MQ=255)

                                            ccGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCt                                             <  2:2991738/93‑1 (MQ=255)

                                              ggAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCg                                           <  2:1338266/93‑1 (MQ=255)

                                              ggAGCAATGCGCGAGGTGAAACAGATGGGACCGGCGATTCTTGGTGGCTCTGTGAGCCAGATCTCCTTAATCATCAACCCCATTTTTGCCTCg                                           <  2:2066854/93‑1 (MQ=255)

                                               gAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGt                                          >  2:1940892/1‑93 (MQ=255)

                                               gAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGt                                          <  1:1823382/93‑1 (MQ=255)

                                                aGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGtt                                         <  2:3836645/93‑1 (MQ=255)

                                                 gCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGttt                                        >  1:3539147/1‑93 (MQ=255)

                                                  cAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTc                                       <  1:575354/93‑1 (MQ=255)

                                                   aaTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCt                                      >  2:4650460/1‑93 (MQ=255)

                                                   aaTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCt                                      <  1:2558983/93‑1 (MQ=255)

                                                   aaTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCt                                      <  1:3558697/93‑1 (MQ=255)

                                                   aaTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCt                                      <  1:4229259/93‑1 (MQ=255)

                                                    aTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCtt                                     <  1:373990/93‑1 (MQ=255)

                                                    aTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCtt                                     <  1:4477071/93‑1 (MQ=255)

                                                     tGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTg                                    >  1:3169427/1‑93 (MQ=255)

                                                     tGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTg                                    >  1:658613/1‑93 (MQ=255)

                                                        gcgTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCtt                                 >  2:432374/1‑93 (MQ=255)

                                                              tgAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGtt                           >  1:4198498/1‑93 (MQ=255)

                                                               gAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTc                          <  2:1595237/93‑1 (MQ=255)

                                                               gAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTc                          <  2:207429/93‑1 (MQ=255)

                                                                 aaCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCgg                        <  2:3115243/93‑1 (MQ=255)

                                                                  aCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgt                       <  2:1371550/93‑1 (MQ=255)

                                                                  aCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgt                       >  1:442370/1‑93 (MQ=255)

                                                                  aCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgt                       >  2:4618935/1‑93 (MQ=255)

                                                                   cAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgtg                      >  2:3592724/1‑93 (MQ=255)

                                                                    aGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgtgt                     >  1:547481/1‑93 (MQ=255)

                                                                     gATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTc                    <  2:2599808/93‑1 (MQ=255)

                                                                     gATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTc                    >  1:1644515/1‑93 (MQ=255)

                                                                     gATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTc                    >  1:3264274/1‑93 (MQ=255)

                                                                       tGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCtt                  <  2:3698837/93‑1 (MQ=255)

                                                                        gggACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTg                 <  1:3716950/93‑1 (MQ=255)

                                                                          gACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGa               <  1:2642675/93‑1 (MQ=255)

                                                                          gACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGa               <  1:1649233/93‑1 (MQ=255)

                                                                              ggCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTa           >  1:1963980/1‑93 (MQ=255)

                                                                               gCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTAt          <  1:385103/93‑1 (MQ=255)

                                                                                 gATCCTTGGTGTCTCTGTGAGCCGGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTa        <  2:672915/93‑1 (MQ=255)

                                                                                 gATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTa        <  2:3581478/93‑1 (MQ=255)

                                                                                 gATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTa        <  2:2038197/93‑1 (MQ=255)

                                                                                 gATCCTTGGGGTCTCTGTGAGTCACATCTCCTTAATCATCAATACCATTTTTGTCTCGTTTCTTGCTTTCGGTTCGGTGTCTTGGGGGTAttt        >  2:4318046/1‑92 (MQ=255)

                                                                                   tCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACg      >  2:1996568/1‑93 (MQ=255)

                                                                                    ccTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGc     >  2:2098873/1‑93 (MQ=255)

                                                                                     cTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGcc    >  1:3407078/1‑93 (MQ=255)

                                                                                      ttGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCg   <  1:405042/93‑1 (MQ=255)

                                                                                      ttGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCg   >  2:4208075/1‑93 (MQ=255)

                                                                                      ttGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCg   >  1:2374136/1‑93 (MQ=255)

                                                                                       tGGTGTCACTATGAGCCAGATCTCCTTAATCATCAACTCCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCGa  <  2:3091031/93‑1 (MQ=255)

GAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCGA  >  NC_000913/1127688‑1127867

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