BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,262,865	G→A	G9D (GGC→GAC)	pspG →	phage shock protein G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,262,865	0	G	A	100.0%	1043.7 / NA	255	G9D (GGC→GAC)	pspG	phage shock protein G
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (151/104); total (151/104)

Only 100 of 255 total aligned reads displayed.

TGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTG  >  NC_000913/4262767‑4262962

tGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACt                                                                                                 <  2:5625620/101‑1 (MQ=255)

 gTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACtt                                                                                                >  2:5625879/1‑101 (MQ=255)

  ttAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACttt                                                                                               >  1:5625880/1‑101 (MQ=255)

   tAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACtttt                                                                                              <  1:5625455/101‑1 (MQ=255)

   tAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACtttt                                                                                              <  2:5625449/101‑1 (MQ=255)

    aaaTCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACttttt                                                                                             >  2:5625884/1‑101 (MQ=255)

     aaTCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACtttttt                                                                                            >  1:5625885/1‑101 (MQ=255)

       tCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGt                                                                                          >  1:5625891/1‑101 (MQ=255)

        cAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTc                                                                                         <  2:5625549/101‑1 (MQ=255)

        cAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTc                                                                                         >  2:5625893/1‑101 (MQ=255)

        cAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTc                                                                                         <  2:5625492/101‑1 (MQ=255)

         aaTAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCa                                                                                        >  1:5625895/1‑101 (MQ=255)

           tAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATg                                                                                      >  2:5625897/1‑101 (MQ=255)

            aGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGc                                                                                     >  1:5625899/1‑101 (MQ=255)

            aGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGc                                                                                     >  1:5625901/1‑101 (MQ=255)

             gATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCt                                                                                    >  1:5625904/1‑101 (MQ=255)

              attattTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTg                                                                                   <  2:5625502/101‑1 (MQ=255)

              attattTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTg                                                                                   <  1:5625554/101‑1 (MQ=255)

              attattTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTg                                                                                   <  2:5625517/101‑1 (MQ=255)

               ttattTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGa                                                                                  >  1:5625908/1‑101 (MQ=255)

                tattTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGAt                                                                                 >  1:5625912/1‑101 (MQ=255)

                  tttttGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATgg                                                                               >  2:5625915/1‑101 (MQ=255)

                  tttttGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATgg                                                                               >  2:5625913/1‑101 (MQ=255)

                   ttttGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGt                                                                              >  1:5625919/1‑101 (MQ=255)

                   ttttGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGt                                                                              <  2:5625552/101‑1 (MQ=255)

                    tttGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTc                                                                             <  1:5625299/101‑1 (MQ=255)

                     ttGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCa                                                                            >  1:5625922/1‑101 (MQ=255)

                      tGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCAc                                                                           >  1:5625924/1‑101 (MQ=255)

                      tGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCAc                                                                           >  2:5625925/1‑101 (MQ=255)

                         cATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCgg                                                                        >  2:5625929/1‑101 (MQ=255)

                          aTGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGc                                                                       <  2:5625535/101‑1 (MQ=255)

                          aTGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGc                                                                       >  2:5625931/1‑101 (MQ=255)

                             aTTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGtt                                                                    >  1:5625934/1‑101 (MQ=255)

                               tCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTc                                                                  <  1:5625676/101‑1 (MQ=255)

                                cTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCg                                                                 <  2:5625367/101‑1 (MQ=255)

                                cTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCg                                                                 >  1:5625938/1‑101 (MQ=255)

                                 ttGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGt                                                                >  1:5625940/1‑101 (MQ=255)

                                 ttGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGt                                                                >  2:5625942/1‑101 (MQ=255)

                                   gTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTg                                                              <  2:5625525/101‑1 (MQ=255)

                                   gTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTg                                                              >  1:5625946/1‑101 (MQ=255)

                                   gTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTg                                                              >  2:5625945/1‑101 (MQ=255)

                                     aaTGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCt                                                            >  2:5625950/1‑101 (MQ=255)

                                      aTGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTg                                                           >  1:5625951/1‑101 (MQ=255)

                                      aTGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTg                                                           <  2:5625568/101‑1 (MQ=255)

                                      aTGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTg                                                           <  1:5625578/101‑1 (MQ=255)

                                       tGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTgg                                                          <  2:5625590/101‑1 (MQ=255)

                                        gCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTggg                                                         >  2:5625955/1‑101 (MQ=255)

                                         ccAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGc                                                        >  1:5625957/1‑101 (MQ=255)

                                           aGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCAt                                                      >  2:5625959/1‑101 (MQ=255)

                                           aGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCAt                                                      <  2:5625636/101‑1 (MQ=255)

                                            gCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCAtt                                                     >  2:5625961/1‑101 (MQ=255)

                                             cAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTa                                                    <  1:5625488/101‑1 (MQ=255)

                                             cAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTa                                                    >  2:5625966/1‑101 (MQ=255)

                                             cAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTa                                                    >  1:5625964/1‑101 (MQ=255)

                                             cAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTa                                                    <  2:5625598/101‑1 (MQ=255)

                                              aaGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTAt                                                   >  2:5625971/1‑101 (MQ=255)

                                              aaGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTAt                                                   >  1:5625970/1‑101 (MQ=255)

                                              aaGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTAt                                                   <  1:5625707/101‑1 (MQ=255)

                                               agagATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATc                                                  <  2:5625713/101‑1 (MQ=255)

                                               agagATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATc                                                  <  2:5625558/101‑1 (MQ=255)

                                                  gaTTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATcgcc                                               <  1:5625498/101‑1 (MQ=255)

                                                  gaTTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATcgcc                                               >  1:5625974/1‑101 (MQ=255)

                                                  gaTTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATcgcc                                               <  2:5625614/101‑1 (MQ=255)

                                                       cATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCTGCGCt                                          <  1:5625651/101‑1 (MQ=255)

                                                       cATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCt                                          <  2:5625683/101‑1 (MQ=255)

                                                       cATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCt                                          >  2:5625977/1‑101 (MQ=255)

                                                        atatTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTg                                         <  1:5625373/101‑1 (MQ=255)

                                                        atatTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTg                                         >  1:5625979/1‑101 (MQ=255)

                                                         tatTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTgg                                        >  2:5625985/1‑101 (MQ=255)

                                                          atTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGt                                       >  1:5625986/1‑101 (MQ=255)

                                                          atTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGt                                       >  1:5625988/1‑101 (MQ=255)

                                                             tGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTtgt                                    <  1:5625746/101‑1 (MQ=255)

                                                             tGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTtgt                                    >  1:5625991/1‑101 (MQ=255)

                                                              gggAGAGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTtgtg                                   <  1:5625446/101‑1 (MQ=255)

                                                                gagaGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGc                                 <  1:5625597/101‑1 (MQ=255)

                                                                gagaGCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGc                                 >  1:5625993/1‑101 (MQ=255)

                                                                   agCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGAc                              <  2:5625604/101‑1 (MQ=255)

                                                                   agCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGAc                              <  2:5625420/101‑1 (MQ=255)

                                                                    gCATCATGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACg                             >  2:5625997/1‑101 (MQ=255)

                                                                     catcatGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTgg                                  >  1:5625999/1‑95 (MQ=255)

                                                                      atcatGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGc                           <  2:5625634/101‑1 (MQ=255)

                                                                      atcatGCTGGAACTACTTTTTGTGATTGACTTTTTTCTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGc                           >  2:5626002/1‑101 (MQ=255)

                                                                        catGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCa                         >  1:5626004/1‑101 (MQ=255)

                                                                         atGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCAt                        >  2:5626005/1‑101 (MQ=255)

                                                                          tGCTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCAtt                       <  1:5625570/101‑1 (MQ=255)

                                                                            cTGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTAt                     <  2:5625738/101‑1 (MQ=255)

                                                                             tGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGc                                                        >  2:5626013/1‑65 (MQ=255)

                                                                             tGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATg                    >  1:5626010/1‑101 (MQ=255)

                                                                             tGGAACTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATg                    >  2:5626011/1‑101 (MQ=255)

                                                                                aaCTACTTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTc                 >  1:5626016/1‑101 (MQ=255)

                                                                                   tactTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTc              >  2:5626019/1‑101 (MQ=255)

                                                                                     ctTTTTGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTcgg            >  2:5626021/1‑101 (MQ=255)

                                                                                      tttttGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTcggc           <  1:5625472/101‑1 (MQ=255)

                                                                                      tttttGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTcggc           <  2:5625586/101‑1 (MQ=255)

                                                                                       ttttGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTcggcg          <  1:5625606/101‑1 (MQ=255)

                                                                                         ttGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGt        >  2:5626025/1‑101 (MQ=255)

                                                                                         ttGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGt        <  1:5625411/101‑1 (MQ=255)

                                                                                         ttGTGATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGt        <  1:5625538/101‑1 (MQ=255)

                                                                                           gtgATTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTAt      >  1:5626027/1‑101 (MQ=255)

                                                                                              aTTGACTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCt   >  1:5626030/1‑101 (MQ=255)

TGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTG  >  NC_000913/4262767‑4262962

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