BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	580,132	C→T	intergenic (+46/+119)	ybcW → / ← ylcI	DLP12 prophage; protein YbcW/DLP12 prophage; DUF3950 domain‑containing protein YlcI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,132	0	C	T	100.0%	1358.2 / NA	320	intergenic (+46/+119)	ybcW/ylcI	DLP12 prophage; protein YbcW/DLP12 prophage; DUF3950 domain‑containing protein YlcI
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (177/143); total (177/143)

Only 100 of 320 total aligned reads displayed.

CGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTC  >  NC_000913/580034‑580229

  ttGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAac                                                                                               >  1:859915/1‑101 (MQ=255)

   tGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAaca                                                                                              >  2:859916/1‑101 (MQ=255)

    ggTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAacac                                                                                             <  1:859451/101‑1 (MQ=255)

      tGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGc                                                                                           >  2:859921/1‑101 (MQ=255)

       gAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCa                                                                                          <  2:859220/101‑1 (MQ=255)

        aaGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGc                                                                                           <  1:859438/99‑1 (MQ=255)

        aaGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAg                                                                                         >  1:859924/1‑101 (MQ=255)

          gCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTc                                                                                       >  2:859928/1‑101 (MQ=255)

          gCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTc                                                                                       <  1:859430/101‑1 (MQ=255)

          gCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTc                                                                                       >  1:859930/1‑101 (MQ=255)

           cACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTct                                                                                      >  2:859933/1‑101 (MQ=255)

           cACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTct                                                                                      <  2:859469/101‑1 (MQ=255)

            aCGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctg                                                                                     >  1:859937/1‑101 (MQ=255)

             cGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctgt                                                                                    <  1:859512/101‑1 (MQ=255)

              gATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctgtc                                                                                   >  1:859941/1‑101 (MQ=255)

               ataataataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctgtca                                                                                  >  2:859942/1‑101 (MQ=255)

               ataataataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctgtca                                                                                  <  1:859474/101‑1 (MQ=255)

                 aataataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctgtcact                                                                                >  1:859947/1‑101 (MQ=255)

                 aataataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctgtcact                                                                                >  2:859950/1‑101 (MQ=255)

                   taataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctgtcactgt                                                                              >  1:859953/1‑101 (MQ=255)

                    aataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctgtcactgtc                                                                             >  2:859956/1‑101 (MQ=255)

                     atatGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctgtcactgtca                                                                            <  1:859379/101‑1 (MQ=255)

                      tatGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAg                                                                           >  2:859962/1‑101 (MQ=255)

                      tatGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAg                                                                           >  1:859961/1‑101 (MQ=255)

                      tatGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAg                                                                           <  2:859579/101‑1 (MQ=255)

                       atGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAgg                                                                          >  2:859966/1‑101 (MQ=255)

                        tGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGa                                                                         <  2:859460/101‑1 (MQ=255)

                        tGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGa                                                                         <  1:858804/101‑1 (MQ=255)

                         gAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTctgtcactgtca                                                                            >  1:859970/1‑97 (MQ=255)

                           aGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAg                                                                      >  1:859974/1‑101 (MQ=255)

                            ggATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGt                                                                     >  2:859976/1‑101 (MQ=255)

                               ttattCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGt                                                                  <  2:859549/101‑1 (MQ=255)

                                tattCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTa                                                                 >  1:859980/1‑101 (MQ=255)

                                tattCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTa                                                                 <  2:859428/101‑1 (MQ=255)

                                tattCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTa                                                                 >  2:859982/1‑101 (MQ=255)

                                 attCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTaa                                                                >  2:859986/1‑101 (MQ=255)

                                  ttCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTaaa                                                               >  1:859990/1‑101 (MQ=255)

                                   tCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTaaaa                                                              >  2:859992/1‑101 (MQ=255)

                                    cccTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAAc                                                             >  2:859996/1‑101 (MQ=255)

                                       tggtggTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGc                                                          <  2:859568/101‑1 (MQ=255)

                                        ggtggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCa                                                         >  2:860000/1‑101 (MQ=255)

                                         gtggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCaa                                                        <  2:859370/101‑1 (MQ=255)

                                         gtggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCaa                                                        >  2:860004/1‑101 (MQ=255)

                                          tggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAAc                                                       <  1:859392/101‑1 (MQ=255)

                                          tggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAAc                                                       >  1:860006/1‑101 (MQ=255)

                                           ggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACt                                                      <  2:859539/101‑1 (MQ=255)

                                           ggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACt                                                      >  2:860008/1‑101 (MQ=255)

                                               gACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAAt                                                  >  1:860012/1‑101 (MQ=255)

                                                aCTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAAtt                                                 >  1:860016/1‑101 (MQ=255)

                                                 cTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTa                                                >  2:860018/1‑101 (MQ=255)

                                                 cTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTa                                                <  1:859343/101‑1 (MQ=255)

                                                  tGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTAc                                               >  2:860022/1‑101 (MQ=255)

                                                   gATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACt                                              >  1:860024/1‑101 (MQ=255)

                                                    aTCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTg                                             <  1:859290/101‑1 (MQ=255)

                                                     tCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACt                                              <  2:859314/99‑1 (MQ=255)

                                                     tCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGc                                            <  2:859453/101‑1 (MQ=255)

                                                      cACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCa                                           >  1:860030/1‑101 (MQ=255)

                                                       aCCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCaa                                          >  2:860032/1‑101 (MQ=255)

                                                       aCCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCaa                                          <  2:859385/101‑1 (MQ=255)

                                                       aCCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCaa                                          <  2:859494/101‑1 (MQ=255)

                                                        ccATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAAt                                         <  2:859621/101‑1 (MQ=255)

                                                        ccATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAAt                                         >  1:860034/1‑101 (MQ=255)

                                                            aaCTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGccc                                     <  2:859652/101‑1 (MQ=255)

                                                            aaCTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGccc                                     >  1:860037/1‑101 (MQ=255)

                                                             aCTGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCt                                    >  2:860041/1‑101 (MQ=255)

                                                               tGCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCg                                  >  2:860046/1‑101 (MQ=255)

                                                                gCTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGt                                 <  1:859605/101‑1 (MQ=255)

                                                                 cTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTa                                >  2:860050/1‑101 (MQ=255)

                                                                 cTAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTa                                >  1:860048/1‑101 (MQ=255)

                                                                  tAATCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTaa                               <  2:859303/101‑1 (MQ=255)

                                                                    aTCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAAtt                             >  1:860055/1‑101 (MQ=255)

                                                                     tCATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTa                            >  2:860056/1‑101 (MQ=255)

                                                                      cATTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTaa                           >  2:860062/1‑101 (MQ=255)

                                                                       aTTCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAg                          <  2:859587/101‑1 (MQ=255)

                                                                        ttCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGt                         <  2:859510/101‑1 (MQ=255)

                                                                         tCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTg                        >  1:860066/1‑101 (MQ=255)

                                                                         tCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTg                        <  1:859422/101‑1 (MQ=255)

                                                                         tCAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTg                        <  1:859872/101‑1 (MQ=255)

                                                                          cAAACTACTTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGa                       >  2:860067/1‑101 (MQ=255)

                                                                             actactTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAAtt                    >  1:860073/1‑101 (MQ=255)

                                                                             actactTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAAtt                    >  2:860072/1‑101 (MQ=255)

                                                                               tactTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTa                  <  1:859548/101‑1 (MQ=255)

                                                                                 ctTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACa                <  2:859743/101‑1 (MQ=255)

                                                                                 ctTAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACa                <  1:859718/101‑1 (MQ=255)

                                                                                   tAACCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAat              >  1:860077/1‑101 (MQ=255)

                                                                                    aaCCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAata             >  2:860080/1‑101 (MQ=255)

                                                                                     aCCTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAatat            <  2:859700/101‑1 (MQ=255)

                                                                                       cTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCg          <  2:859879/101‑1 (MQ=255)

                                                                                       cTGTGACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCg          <  2:859602/101‑1 (MQ=255)

                                                                                        tgtgACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGt         >  1:860082/1‑101 (MQ=255)

                                                                                          tgACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTcc       <  1:859595/101‑1 (MQ=255)

                                                                                          tgACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTcc       >  2:860085/1‑101 (MQ=255)

                                                                                           gACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCt      <  1:859775/101‑1 (MQ=255)

                                                                                           gACAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCt      <  2:859643/101‑1 (MQ=255)

                                                                                            aCAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTg     <  1:859665/101‑1 (MQ=255)

                                                                                            aCAGAGTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTg     <  2:859328/101‑1 (MQ=255)

                                                                                              agagTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGtt   <  1:859687/101‑1 (MQ=255)

                                                                                              agagTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGtt   <  2:859532/101‑1 (MQ=255)

                                                                                              agagTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGtt   >  1:860090/1‑101 (MQ=255)

                                                                                               gagTCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTc  <  2:859708/101‑1 (MQ=255)

CGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTC  >  NC_000913/580034‑580229

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