BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,209,141	T→C	pseudogene (515/537 nt)	stfE ←	e14 prophage; putative side tail fiber protein fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,209,141	0	T	C	100.0%	965.9 / NA	233	pseudogene (515/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (122/111); total (122/111)

Only 100 of 233 total aligned reads displayed.

TCACCATCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGTTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATAATCCGGGGCTGG  >  NC_000913/1209043‑1209238

tCACCATCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTa                                                                                                 <  1:1578801/101‑1 (MQ=255)

 cACCATCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTaa                                                                                                <  2:1578775/101‑1 (MQ=255)

  aCCATCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTaaa                                                                                               >  1:1579260/1‑101 (MQ=255)

  aCCATCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTaaa                                                                                               >  2:1579258/1‑101 (MQ=255)

  aCCATCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTaaa                                                                                               <  1:1578930/101‑1 (MQ=255)

  aCCATCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTaaa                                                                                               <  1:1579007/101‑1 (MQ=255)

   ccATCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAAt                                                                                              <  2:1578757/101‑1 (MQ=255)

    cATCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATg                                                                                             <  1:1578745/101‑1 (MQ=255)

     aTCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGc                                                                                            <  1:1578765/101‑1 (MQ=255)

     aTCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGc                                                                                            >  1:1579262/1‑101 (MQ=255)

      tCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCa                                                                                           >  1:1579264/1‑101 (MQ=255)

       cGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCaa                                                                                          >  2:1579265/1‑101 (MQ=255)

        gTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAat                                                                                         <  1:1578911/101‑1 (MQ=255)

         tAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAata                                                                                        <  1:1579068/101‑1 (MQ=255)

          aGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAatat                                                                                       >  2:1579268/1‑101 (MQ=255)

           gTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATAtt                                                                                      >  1:1579270/1‑101 (MQ=255)

            ttATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATAtt                                                                                      <  1:1578828/100‑1 (MQ=255)

            ttATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTa                                                                                     <  2:1578832/101‑1 (MQ=255)

            ttATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTa                                                                                     >  1:1579272/1‑101 (MQ=255)

             tataAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTAc                                                                                    <  1:1578843/101‑1 (MQ=255)

              ataAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTAcg                                                                                   >  2:1579277/1‑101 (MQ=255)

              ataAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTAcg                                                                                   >  1:1579275/1‑101 (MQ=255)

                aaaCGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTAcgcg                                                                                 <  1:1578961/101‑1 (MQ=255)

                aaaCGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTAcgcg                                                                                 <  2:1578964/101‑1 (MQ=255)

                  aCGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGa                                                                               <  2:1578869/101‑1 (MQ=255)

                   cGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGAc                                                                              <  2:1578898/101‑1 (MQ=255)

                       aaCATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGt                                                                          >  2:1579285/1‑101 (MQ=255)

                       aaCATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGt                                                                          >  2:1579282/1‑101 (MQ=255)

                       aaCATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGt                                                                          >  1:1579283/1‑101 (MQ=255)

                        aCATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGtt                                                                         <  1:1579133/101‑1 (MQ=255)

                        aCATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGtt                                                                         <  2:1578872/101‑1 (MQ=255)

                         cATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGttt                                                                        <  2:1578934/101‑1 (MQ=255)

                           tCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCt                                                                      >  1:1579289/1‑101 (MQ=255)

                            ccccTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTg                                                                     >  1:1579291/1‑101 (MQ=255)

                              ccTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAg                                                                   >  2:1579293/1‑101 (MQ=255)

                              ccTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAg                                                                   <  2:1578999/101‑1 (MQ=255)

                               cTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAgg                                                                  <  2:1579161/101‑1 (MQ=255)

                                tGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGc                                                                 >  1:1579297/1‑101 (MQ=255)

                                tGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGc                                                                 >  1:1579295/1‑101 (MQ=255)

                                tGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGc                                                                 <  1:1578865/101‑1 (MQ=255)

                                 gCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCt                                                                >  1:1579299/1‑101 (MQ=255)

                                   cTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGc                                                              <  2:1578982/101‑1 (MQ=255)

                                   cTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGc                                                              >  2:1579303/1‑101 (MQ=255)

                                   cTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGc                                                              >  2:1579305/1‑101 (MQ=255)

                                   cTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGc                                                              >  1:1579306/1‑101 (MQ=255)

                                         tgCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCAcc                                                        >  2:1579310/1‑101 (MQ=255)

                                         tgCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCAcc                                                        <  2:1578944/101‑1 (MQ=255)

                                          gCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCt                                                       <  1:1578996/101‑1 (MQ=255)

                                            aaTAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTaa                                                     >  1:1579312/1‑101 (MQ=255)

                                            aaTAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTaa                                                     >  2:1579313/1‑101 (MQ=255)

                                             aTAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTaaa                                                    <  1:1578698/101‑1 (MQ=255)

                                             aTAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTaaa                                                    >  1:1579316/1‑101 (MQ=255)

                                             aTAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTaaa                                                    <  1:1579241/101‑1 (MQ=255)

                                               aaGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAAcc                                                  >  2:1579317/1‑101 (MQ=255)

                                                 ggTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCAt                                                >  1:1579321/1‑101 (MQ=255)

                                                    cACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCa                                             >  1:1579323/1‑101 (MQ=255)

                                                      cTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACt                                           >  1:1579326/1‑101 (MQ=255)

                                                       tATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTg                                          <  1:1578977/101‑1 (MQ=255)

                                                       tATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTg                                          >  2:1579329/1‑101 (MQ=255)

                                                       tATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTg                                          <  2:1579073/101‑1 (MQ=255)

                                                       tATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTg                                          <  1:1578903/101‑1 (MQ=255)

                                                       tATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTg                                          >  2:1579327/1‑101 (MQ=255)

                                                         tttAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGAt                                        >  2:1579331/1‑101 (MQ=255)

                                                             atatTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGtt                                    >  2:1579336/1‑101 (MQ=255)

                                                               atTGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTta                                  >  1:1579337/1‑101 (MQ=255)

                                                                ttGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCaa                                                                                          >  2:1579343/1‑44 (MQ=255)

                                                                ttGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTtat                                 >  2:1579340/1‑101 (MQ=255)

                                                                ttGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTtat                                 >  1:1579341/1‑101 (MQ=255)

                                                                 tGCTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTtata                                <  2:1579126/101‑1 (MQ=255)

                                                                   cTTTATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATg                              <  2:1579091/101‑1 (MQ=255)

                                                                    tttATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGt                             >  2:1579346/1‑101 (MQ=255)

                                                                      tATGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGttt                           <  2:1579203/101‑1 (MQ=255)

                                                                       aTGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGtttt                          >  2:1579350/1‑101 (MQ=255)

                                                                       aTGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGtttt                          >  1:1579349/1‑101 (MQ=255)

                                                                       aTGCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGtttt                          <  1:1578925/101‑1 (MQ=255)

                                                                         gCATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTaa                        <  1:1579056/101‑1 (MQ=255)

                                                                          cATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTaaa                       >  2:1579352/1‑101 (MQ=255)

                                                                          cATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTaaa                       <  1:1578993/101‑1 (MQ=255)

                                                                          cATCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTaaa                       <  2:1578921/101‑1 (MQ=255)

                                                                            tCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAgg                     >  1:1579355/1‑101 (MQ=255)

                                                                            tCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAgg                     <  1:1579149/101‑1 (MQ=255)

                                                                            tCAGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAgg                     <  1:1579075/101‑1 (MQ=255)

                                                                              aGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGtt                   >  1:1577357/1‑101 (MQ=255)

                                                                              aGGCTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGtt                   >  1:1579358/1‑101 (MQ=255)

                                                                              aGGCTGCCCTCACGATATAGCTAAAGGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGtt                   >  1:1579361/1‑101 (MQ=255)

                                                                                 cTGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCa                >  2:1579362/1‑101 (MQ=255)

                                                                                  tGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCAt               <  1:1578882/101‑1 (MQ=255)

                                                                                  tGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCAt               <  2:1578972/101‑1 (MQ=255)

                                                                                  tGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCAt               >  1:1579366/1‑101 (MQ=255)

                                                                                  tGCCCTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCAt               <  2:1578743/101‑1 (MQ=255)

                                                                                    cccTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATaa             <  2:1579019/101‑1 (MQ=255)

                                                                                     ccTCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATAAt            >  2:1579368/1‑101 (MQ=255)

                                                                                       tCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCAt               <  2:1579029/96‑1 (MQ=255)

                                                                                       tCACGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATAATcc          >  1:1579369/1‑101 (MQ=255)

                                                                                        cACGATATAGCTAAATGCAATATTACGCGGACGCTTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATAATCCg         <  2:1578913/101‑1 (MQ=255)

                                                                                         aCGATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATAATCCgg        <  1:1578940/101‑1 (MQ=255)

                                                                                           gATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATAATCCgggg      <  2:1578974/101‑1 (MQ=255)

                                                                                           gATATAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATAATCCgggg      >  1:1579372/1‑101 (MQ=255)

                                                                                            atatAGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATAATCCGGGGc     >  2:1579373/1‑101 (MQ=255)

                                                                                              ataGCTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATAATCCGGGGCTg   <  2:1579052/101‑1 (MQ=255)

TCACCATCGTAGTTATAAACGGTAACATCCCCTGCCTTTGTTGCAATAAGGTCACTATTTAATATTGCTTTATGCATCAGGCTGCCCTCACGATATAGTTAAATGCAATATTACGCGGACGCGTTTCTGAGGCTGCGGCACCTAAACCATCCACTGATTGTTTATATGTTTTAAAGGTTCCATAATCCGGGGCTGG  >  NC_000913/1209043‑1209238

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