BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,560,020	T→C	I388T (ATT→ACT)	intZ →	CPZ‑55 prophage; putative phage integrase IntZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,560,020	0	T	C	100.0%	1261.8 / NA	305	I388T (ATT→ACT)	intZ	CPZ‑55 prophage; putative phage integrase IntZ
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (151/154); total (151/154)

Only 100 of 305 total aligned reads displayed.

CCAACTAGCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTAT  >  NC_000913/2559922‑2560117

 cAACTAGCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGa                                                                                                <  2:3314246/101‑1 (MQ=255)

 cAACTAGCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGa                                                                                                >  2:3314746/1‑101 (MQ=255)

  aaCTAGCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGAt                                                                                               <  2:3314348/101‑1 (MQ=255)

  aaCTAGCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGAt                                                                                               >  2:3314750/1‑101 (MQ=255)

  aaCTAGCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGAt                                                                                               <  1:3314425/101‑1 (MQ=255)

   aCTAGCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGAtt                                                                                              >  1:3314754/1‑101 (MQ=255)

      aGCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCaa                                                                                           >  1:3314758/1‑101 (MQ=255)

      aGCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCaa                                                                                           <  1:3314298/101‑1 (MQ=255)

       gCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAat                                                                                          >  2:3314760/1‑101 (MQ=255)

        cccACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCGATGATGCAGGATTACGCCGACGCTACTGATTCAata                                                                                         <  2:3314027/101‑1 (MQ=255)

        cccACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAata                                                                                         >  1:3314765/1‑101 (MQ=255)

          cACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATAtt                                                                                       >  1:3314768/1‑101 (MQ=255)

            cGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTgg                                                                                     >  2:3314771/1‑101 (MQ=255)

             gTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGc                                                                                    >  1:3314773/1‑101 (MQ=255)

                ccGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAgg                                                                                 >  1:3314778/1‑101 (MQ=255)

                ccGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAgg                                                                                 >  2:3314777/1‑101 (MQ=255)

                  gtaGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTa                                                                               >  2:3314783/3‑101 (MQ=255)

                  gAAGGGGATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTa                                                                               >  1:3314784/1‑101 (MQ=255)

                   aaGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAgg                                                                                 <  1:3314193/98‑1 (MQ=255)

                   aaGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTaa                                                                              >  1:3314788/1‑101 (MQ=255)

                   aaGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTaa                                                                              <  2:3314378/101‑1 (MQ=255)

                    aGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAAc                                                                             <  1:3314385/101‑1 (MQ=255)

                      ggAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACgg                                                                           >  2:3314793/1‑101 (MQ=255)

                      ggAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACgg                                                                           >  1:3314792/1‑101 (MQ=255)

                         aTAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCaa                                                                        >  2:3314798/1‑101 (MQ=255)

                         aTAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCaa                                                                        >  1:3314800/1‑101 (MQ=255)

                          tAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAAt                                                                       >  1:3314805/1‑101 (MQ=255)

                          tAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAAt                                                                       >  2:3314806/1‑101 (MQ=255)

                           aaGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATc                                                                      <  1:3314437/101‑1 (MQ=255)

                              gggCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTg                                                                                      <  1:3314120/82‑1 (MQ=255)

                               ggCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCctct                                                                  >  1:3314812/1‑101 (MQ=255)

                               ggCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCctct                                                                  <  2:3314316/101‑1 (MQ=255)

                               ggCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCctct                                                                  >  2:3314811/1‑101 (MQ=255)

                                gCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTa                                                                 <  1:3314606/101‑1 (MQ=255)

                                gCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTa                                                                 <  1:3314247/101‑1 (MQ=255)

                                 cGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAg                                                                <  2:3314397/101‑1 (MQ=255)

                                  gTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGa                                                               <  1:3314498/101‑1 (MQ=255)

                                   tttATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGaa                                                              >  1:3314816/1‑101 (MQ=255)

                                    ttATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAAc                                                             >  2:3314819/1‑101 (MQ=255)

                                     tataACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAAcc                                                            <  1:3314337/101‑1 (MQ=255)

                                     tataACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAAcc                                                            <  2:3314258/101‑1 (MQ=255)

                                      ataACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCg                                                           >  1:3314821/1‑101 (MQ=255)

                                       taACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCgg                                                          >  1:3314824/1‑101 (MQ=255)

                                        aaCCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGga                                                         <  1:3314380/101‑1 (MQ=255)

                                         aCCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaa                                                        >  1:3314829/1‑101 (MQ=255)

                                           cAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatg                                                      >  2:3314830/1‑101 (MQ=255)

                                           cAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatg                                                      <  2:3314521/101‑1 (MQ=255)

                                             ggCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatgaa                                                    <  1:3314160/101‑1 (MQ=255)

                                             ggCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatgaa                                                    <  1:3314516/101‑1 (MQ=255)

                                             ggCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatgaa                                                    <  1:3314314/101‑1 (MQ=255)

                                              gCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatgaat                                                   >  2:3314833/1‑101 (MQ=255)

                                               cTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGaata                                                  <  2:3314290/101‑1 (MQ=255)

                                                  ccTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAAtt                                               >  1:3314837/1‑101 (MQ=255)

                                                    tACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTaa                                             >  1:3314840/1‑101 (MQ=255)

                                                    tACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTaa                                             <  1:3314472/101‑1 (MQ=255)

                                                      ccTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGc                                           <  2:3314570/101‑1 (MQ=255)

                                                      ccTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGc                                           >  2:3314844/1‑101 (MQ=255)

                                                       cTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCg                                          <  2:3314410/101‑1 (MQ=255)

                                                       cTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCg                                          >  1:3314845/1‑101 (MQ=255)

                                                        tGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAgg                                                                                 >  2:3314851/1‑61 (MQ=255)

                                                        tGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCgg                                         >  2:3314848/1‑101 (MQ=255)

                                                        tGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCgg                                         <  1:3314288/101‑1 (MQ=255)

                                                              aaCGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCttttt                                   >  2:3314856/1‑101 (MQ=255)

                                                              aaCGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCttttt                                   >  1:3314857/1‑101 (MQ=255)

                                                               aCGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCtttttt                                  <  2:3314434/101‑1 (MQ=255)

                                                                cGGCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTg                                 <  1:3314177/101‑1 (MQ=255)

                                                                 ggCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGc                                >  1:3314862/1‑101 (MQ=255)

                                                                 ggCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGc                                <  1:3314391/101‑1 (MQ=255)

                                                                 ggCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGc                                <  1:3314203/101‑1 (MQ=255)

                                                                  gCGAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCa                               <  1:3314420/101‑1 (MQ=255)

                                                                    gAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAAtt                                               <  2:3314331/83‑1 (MQ=255)

                                                                    gAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCaaa                             <  1:3314485/101‑1 (MQ=255)

                                                                    gAGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCaaa                             >  1:3314865/1‑101 (MQ=255)

                                                                     aGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGc                                <  2:3314220/97‑1 (MQ=255)

                                                                     aGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAg                            >  2:3314867/1‑101 (MQ=255)

                                                                     aGCAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAg                            <  1:3314350/101‑1 (MQ=255)

                                                                       cAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgcg                          >  2:3314871/1‑101 (MQ=255)

                                                                       cAATGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgcg                          >  1:3314873/1‑101 (MQ=255)

                                                                        aaTGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgcgc                         >  2:3314875/1‑101 (MQ=255)

                                                                        aaTGATGCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgcgc                         <  1:3314568/101‑1 (MQ=255)

                                                                          tgatgCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCaa                       <  2:3314464/101‑1 (MQ=255)

                                                                           gatgCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAAt                      >  2:3314880/1‑101 (MQ=255)

                                                                            atgCAGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAAtt                     >  1:3314883/1‑101 (MQ=255)

                                                                                aGGATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGc                 <  1:3314261/101‑1 (MQ=255)

                                                                                  gATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCAt               >  2:3314886/1‑101 (MQ=255)

                                                                                  gATTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCAt               <  2:3314497/101‑1 (MQ=255)

                                                                                   aTTACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCAtt              <  1:3314623/101‑1 (MQ=255)

                                                                                    ttACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCAttt             <  2:3314469/101‑1 (MQ=255)

                                                                                     tACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCAtttt            <  1:3314360/101‑1 (MQ=255)

                                                                                     tACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCAtttt            >  2:3314891/1‑101 (MQ=255)

                                                                                     tACGCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCAtttt            >  1:3314887/1‑101 (MQ=255)

                                                                                        gCCGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGta         <  1:3314413/101‑1 (MQ=255)

                                                                                        gCCGACGCTACTGATTCAATAGTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAg                            >  1:3314895/1‑82 (MQ=255)

                                                                                         ccGACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGtat        <  2:3314552/101‑1 (MQ=255)

                                                                                           gACGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATaa      <  2:3314668/101‑1 (MQ=255)

                                                                                            aCGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAg     <  1:3314445/101‑1 (MQ=255)

                                                                                            aCGCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAg     <  1:3314681/101‑1 (MQ=255)

                                                                                              gCTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTa   <  2:3314618/101‑1 (MQ=255)

                                                                                               cTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTAt  <  1:3314457/101‑1 (MQ=255)

                                                                                               cTACTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTAt  <  2:3314453/101‑1 (MQ=255)

CCAACTAGCCCACGTCCCGAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTAT  >  NC_000913/2559922‑2560117

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