BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,207,235	C→T	A108V (GCC→GTC)	ymfQ →	e14 prophage; conserved protein YmfQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,207,235	0	C	T	100.0%	959.5 / NA	234	A108V (GCC→GTC)	ymfQ	e14 prophage; conserved protein YmfQ
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (107/127); total (107/127)

Only 100 of 248 total aligned reads displayed.

GGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGG  >  NC_000913/1207137‑1207332

gggAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGtcc                                                                                                 >  1:1576584/1‑101 (MQ=255)

 ggAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCt                                                                                                <  1:1576344/101‑1 (MQ=255)

  gAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTg                                                                                               >  1:1576586/1‑101 (MQ=255)

    aCGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTggg                                                                                             >  2:1576588/1‑101 (MQ=255)

     cGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGc                                                                                            <  1:1576416/101‑1 (MQ=255)

      gCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCa                                                                                           >  1:1576591/1‑101 (MQ=255)

      gCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCa                                                                                           >  2:1576592/1‑101 (MQ=255)

       cAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAg                                                                                          >  1:1576595/1‑101 (MQ=255)

       cAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAg                                                                                          <  2:1576074/101‑1 (MQ=255)

        aGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGa                                                                                         <  2:1576287/101‑1 (MQ=255)

          aCCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGAcc                                                                                       <  1:1576382/101‑1 (MQ=255)

           cccTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCa                                                                                      <  2:1576151/101‑1 (MQ=255)

             cTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGa                                                                                    <  2:1576229/101‑1 (MQ=255)

             cTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGa                                                                                    >  2:1576600/1‑101 (MQ=255)

               tCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGAt                                                                                   <  1:1576325/100‑1 (MQ=255)

                  ccAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGcca                                                                               >  1:1576603/1‑101 (MQ=255)

                  ccAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGcca                                                                               <  1:1575987/101‑1 (MQ=255)

                   cAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGccac                                                                              >  2:1576604/1‑101 (MQ=255)

                    aGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGccacc                                                                             >  1:1576607/1‑101 (MQ=255)

                    aGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGccacc                                                                             <  1:1576269/101‑1 (MQ=255)

                       gTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTg                                                                                               <  2:1576169/80‑1 (MQ=255)

                       gTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATc                                                                          >  2:1576611/1‑101 (MQ=255)

                          aGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCAcg                                                                       >  2:1576614/1‑101 (MQ=255)

                           gCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCAcgc                                                                      <  1:1576234/101‑1 (MQ=255)

                            cAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCAcgcg                                                                     <  1:1576215/101‑1 (MQ=255)

                               cGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATa                                                                  <  1:1576392/101‑1 (MQ=255)

                               cGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATa                                                                  <  2:1576324/101‑1 (MQ=255)

                               cGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATa                                                                  >  2:1576621/1‑101 (MQ=255)

                                ggCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATAc                                                                 <  2:1576135/101‑1 (MQ=255)

                                 gCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACg                                                                <  2:1576289/101‑1 (MQ=255)

                                   tGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGAc                                                              >  1:1576622/1‑101 (MQ=255)

                                     gATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCAcg                                                                       >  2:1576625/1‑90 (MQ=255)

                                      aTGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACaaa                                                           >  1:1576626/1‑101 (MQ=255)

                                       tGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACaaaa                                                          <  2:1576350/101‑1 (MQ=255)

                                          gAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCa                                                       >  2:1576632/1‑101 (MQ=255)

                                          gAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCa                                                       >  1:1576631/1‑101 (MQ=255)

                                          gAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCa                                                       <  2:1576342/101‑1 (MQ=255)

                                                tAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCa                                                 <  2:1576263/101‑1 (MQ=255)

                                                tAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTCCa                                                 <  1:1576222/101‑1 (MQ=255)

                                                 aaCCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCAc                                                >  1:1576635/1‑101 (MQ=255)

                                                 aaCCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCAc                                                <  1:1576384/101‑1 (MQ=255)

                                                  aCCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCAcc                                               >  2:1576636/1‑101 (MQ=255)

                                                   ccTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCt                                              <  1:1576362/101‑1 (MQ=255)

                                                   ccTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCt                                              <  1:1576264/101‑1 (MQ=255)

                                                    cTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTg                                             >  1:1576640/1‑101 (MQ=255)

                                                    cTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTg                                             >  1:1576642/1‑101 (MQ=255)

                                                    cTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTg                                             <  1:1576205/101‑1 (MQ=255)

                                                     tGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGc                                            <  2:1576297/101‑1 (MQ=255)

                                                      ggcgggcgGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCt                                           <  1:1576252/101‑1 (MQ=255)

                                                       gcgggcggCATCAATGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTc                                          <  2:1576334/101‑1 (MQ=255)

                                                        cgggcggCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCa                                         >  2:1576646/1‑101 (MQ=255)

                                                         gggcggCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCAt                                        >  1:1576648/1‑101 (MQ=255)

                                                         gggcggCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCAt                                        >  2:1576649/1‑101 (MQ=255)

                                                           gcggcATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCg                                      <  2:1576498/101‑1 (MQ=255)

                                                           gcggcATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCg                                      >  1:1576651/1‑101 (MQ=255)

                                                           gcggcATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCg                                      <  2:1576355/101‑1 (MQ=255)

                                                           gcggcATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCg                                      <  1:1576254/101‑1 (MQ=255)

                                                             ggcATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGATCATCGGc                                    >  2:1576654/1‑101 (MQ=255)

                                                              gcATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCCCGCGATACGACAAAAGCACTTTCACCTGCTCATc                                       <  2:1576445/97‑1 (MQ=255)

                                                              gcATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGcc                                   <  1:1576176/101‑1 (MQ=255)

                                                              gcATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGcc                                   >  1:1576657/1‑101 (MQ=255)

                                                               cATCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCt                                  <  2:1576389/101‑1 (MQ=255)

                                                                aTCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTg                                 >  1:1576659/1‑101 (MQ=255)

                                                                 tCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGt                                >  2:1576662/1‑101 (MQ=255)

                                                                 tCAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGt                                <  2:1576394/101‑1 (MQ=255)

                                                                  cAACGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTa                               <  2:1576419/101‑1 (MQ=255)

                                                                   aaCGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCa                                                                                           >  2:1576666/1‑40 (MQ=255)

                                                                    aCGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACt                             >  1:1576667/1‑101 (MQ=255)

                                                                     cGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTg                            >  2:1576669/1‑101 (MQ=255)

                                                                     cGAGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTg                            <  1:1576140/101‑1 (MQ=255)

                                                                       aGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGAc                          >  2:1576672/1‑101 (MQ=255)

                                                                       aGGATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGAc                          <  1:1576066/101‑1 (MQ=255)

                                                                         gATTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGAcgc                        >  2:1576674/1‑101 (MQ=255)

                                                                         gATCTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGAcgc                        <  1:1576283/101‑1 (MQ=255)

                                                                          aTTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGAcgcg                       <  1:1576158/101‑1 (MQ=255)

                                                                          aTTTTTATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGAcgcg                       <  2:1576427/101‑1 (MQ=255)

                                                                           tttttATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCgg                      <  1:1576314/101‑1 (MQ=255)

                                                                           tttttATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCgg                      <  1:1576354/101‑1 (MQ=255)

                                                                            ttttATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGt                     >  1:1576677/1‑101 (MQ=255)

                                                                             tttATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTg                    <  2:1576189/101‑1 (MQ=255)

                                                                             tttATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTg                    >  2:1576679/1‑101 (MQ=255)

                                                                             tttATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTg                    >  1:1576682/1‑101 (MQ=255)

                                                                              ttATCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGa                   >  2:1576684/1‑101 (MQ=255)

                                                                                aTCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAAt                 <  2:1576458/101‑1 (MQ=255)

                                                                                 tCTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATg                >  2:1576687/1‑101 (MQ=255)

                                                                                  cTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATgc               <  2:1576313/101‑1 (MQ=255)

                                                                                  cTTGCACAGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATgc               >  1:1576690/1‑101 (MQ=255)

                                                                                      cacaGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCg           <  2:1576261/101‑1 (MQ=255)

                                                                                        caGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGa         >  2:1576693/1‑101 (MQ=255)

                                                                                        caGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGa         <  1:1576292/101‑1 (MQ=255)

                                                                                         aGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACg                                                                <  1:1576467/45‑1 (MQ=255)

                                                                                         aGCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCg           <  2:1576216/98‑1 (MQ=255)

                                                                                          gCTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAAt       <  2:1576302/101‑1 (MQ=255)

                                                                                           cTTGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATg      >  1:1576695/1‑101 (MQ=255)

                                                                                            ttGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATgcgc             >  2:1576698/1‑93 (MQ=255)

                                                                                            ttGCTGTCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATgg     <  2:1576245/101‑1 (MQ=255)

                                                                                             tgctgtCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGc    >  1:1576699/1‑101 (MQ=255)

                                                                                             tgctgtCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGc    >  1:1576704/1‑101 (MQ=255)

                                                                                             tgctgtCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGc    >  2:1576701/1‑101 (MQ=255)

                                                                                              gctgtCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCg   <  1:1576347/101‑1 (MQ=255)

GGGAACGCAGACCCTTCGCCAGCGTCAGCAACGGCTGGATGCGAAGGTTAACCTGGCGGGCGGCATCAACGAGGATTTTTATCTTGCACAGCTTGCTGCCCTGGGCAGACCAGATGCCACCATCACGCGATACGACAAAAGCACTTTCACCTGCTCATCGGCCTGTACTGACGCGGTGAATGCGCCGGAATGGCGG  >  NC_000913/1207137‑1207332

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