BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,261,695	C→T	A61T (GCC→ACC)	prs ←	ribose‑phosphate diphosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,261,695	0	C	T	100.0%	951.5 / NA	237	A61T (GCC→ACC)	prs	ribose‑phosphate diphosphokinase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (89/148); total (89/148)

Only 100 of 244 total aligned reads displayed.

GCCAAAGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCG  >  NC_000913/1261597‑1261792

  cAAAGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGtaca                                                                                                <  2:1647151/100‑1 (MQ=255)

  cAAAGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACaa                                                                                               >  1:1647792/1‑101 (MQ=255)

   aaaGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAg                                                                                              >  2:1647794/1‑101 (MQ=255)

    aaGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGt                                                                                             >  2:1647797/1‑101 (MQ=255)

     aGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTg                                                                                            <  2:1647446/101‑1 (MQ=255)

       tAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGa                                                                                          >  2:1647800/1‑101 (MQ=255)

         ggggATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACt                                                                                        >  1:1647801/1‑101 (MQ=255)

         ggggATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACt                                                                                        <  1:1647479/101‑1 (MQ=255)

          gggATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTg                                                                                       <  1:1647350/101‑1 (MQ=255)

          gggATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTg                                                                                       <  1:1647456/101‑1 (MQ=255)

          gggATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTg                                                                                       <  1:1647474/101‑1 (MQ=255)

           ggATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTgg                                                                                      <  1:1647443/101‑1 (MQ=255)

            gATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGga                                                                                     >  1:1647804/1‑101 (MQ=255)

            gATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGga                                                                                     <  1:1647500/101‑1 (MQ=255)

               aaCAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGgatga                                                                                  >  1:1647806/1‑101 (MQ=255)

                aCAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGgatgat                                                                                 <  1:1647223/101‑1 (MQ=255)

                aCAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGgatgat                                                                                 >  2:1647807/1‑101 (MQ=255)

                  aGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGAtgatga                                                                               >  2:1647811/1‑101 (MQ=255)

                   gCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGaa                                                                              >  1:1647813/1‑101 (MQ=255)

                   gCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGaa                                                                              <  1:1647340/101‑1 (MQ=255)

                     ggTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGaaaa                                                                            >  1:1647815/1‑101 (MQ=255)

                         aTACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATc                                                                        <  1:1647422/101‑1 (MQ=255)

                          tACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATca                                                                       >  2:1647817/1‑101 (MQ=255)

                            cGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAat                                                                           <  2:1647418/95‑1 (MQ=255)

                            cGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATcacc                                                                     >  2:1647821/1‑101 (MQ=255)

                             gACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATcacca                                                                    >  2:1647824/1‑101 (MQ=255)

                             gACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATcacca                                                                    <  2:1647484/101‑1 (MQ=255)

                               ccTGAGGAAGCACGACGCAGGGCATCAACCATAACGTCTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATca                                                                       <  2:1647499/96‑1 (MQ=255)

                                cTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACcg                                                                 >  2:1647826/1‑101 (MQ=255)

                                 tGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACcgc                                                                >  1:1647827/1‑101 (MQ=255)

                                 tGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACcgc                                                                >  2:1647830/1‑101 (MQ=255)

                                 tGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACcgc                                                                <  1:1647396/101‑1 (MQ=255)

                                  gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACcgcg                                                               >  2:1647833/1‑101 (MQ=255)

                                   cGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGt                                                              <  2:1647489/101‑1 (MQ=255)

                                     gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTAc                                                            >  1:1647835/1‑101 (MQ=255)

                                      aaGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACa                                                           <  2:1647369/101‑1 (MQ=255)

                                        gCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACAtt                                                         <  2:1647523/101‑1 (MQ=255)

                                        gCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACAtt                                                         <  2:1647672/101‑1 (MQ=255)

                                         cACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACAttt                                                        <  1:1647561/101‑1 (MQ=255)

                                         cACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACAttt                                                        >  2:1647837/1‑101 (MQ=255)

                                         cACGACGCAGGCCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGaaa                                                                             <  1:1647387/80‑1 (MQ=255)

                                           cgacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTc                                                      <  2:1647333/101‑1 (MQ=255)

                                            gacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCCCCACCGCGTACATTTTCa                                                     <  1:1647580/101‑1 (MQ=255)

                                            gacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCa                                                     >  2:1647839/1‑101 (MQ=255)

                                            gacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCa                                                     <  2:1647351/101‑1 (MQ=255)

                                             acgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCAt                                                    <  2:1647637/101‑1 (MQ=255)

                                              cgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCAtt                                                   <  2:1647531/101‑1 (MQ=255)

                                              cgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCAtt                                                   <  2:1647405/101‑1 (MQ=255)

                                               gCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACGAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTa                                                  <  1:1647128/101‑1 (MQ=255)

                                               gCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTa                                                  <  2:1647515/101‑1 (MQ=255)

                                               acAGGGCATCAACCATAACGACTAATTCCATCTGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTa                                                  <  2:1647510/100‑1 (MQ=255)

                                                cAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTaa                                                 >  2:1647842/1‑101 (MQ=255)

                                                 aGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAAt                                                >  1:1647846/1‑101 (MQ=255)

                                                  gggCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAAtt                                               <  1:1647606/101‑1 (MQ=255)

                                                    gCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTg                                             <  2:1647542/101‑1 (MQ=255)

                                                    gCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTg                                             <  1:1647311/101‑1 (MQ=255)

                                                      aTCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTa                                           <  2:1647291/101‑1 (MQ=255)

                                                       tCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTAc                                          <  2:1647493/101‑1 (MQ=255)

                                                       tCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTAc                                          >  1:1647848/1‑101 (MQ=255)

                                                         aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGc                                        >  2:1647851/1‑101 (MQ=255)

                                                          aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCt                                       <  1:1647540/101‑1 (MQ=255)

                                                          aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCt                                       >  1:1647853/1‑101 (MQ=255)

                                                          aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCt                                       <  2:1647614/101‑1 (MQ=255)

                                                           ccATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTg                                      <  2:1647466/101‑1 (MQ=255)

                                                           ccATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTg                                      <  1:1647567/101‑1 (MQ=255)

                                                            cATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGa                                     <  2:1647645/101‑1 (MQ=255)

                                                                aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTc                                 <  1:1647512/101‑1 (MQ=255)

                                                                aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTc                                 <  1:1647415/101‑1 (MQ=255)

                                                                  gACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGc                               <  1:1647437/101‑1 (MQ=255)

                                                                  gACTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGGGGACTGGATGATGaaa                                                                             >  2:1647855/1‑55 (MQ=255)

                                                                    cTAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACcgcg                                                               >  1:1647859/1‑67 (MQ=255)

                                                                     tAATTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCAt                            >  2:1647860/1‑101 (MQ=255)

                                                                      aaTTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTGATTTGTACGCTGACTTCGCCATc                           >  1:1647862/1‑101 (MQ=255)

                                                                       aTTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCg                          <  2:1647550/101‑1 (MQ=255)

                                                                       aTTCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCg                          >  2:1647865/1‑101 (MQ=255)

                                                                        ttCCATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGc                         <  2:1647431/101‑1 (MQ=255)

                                                                          ccATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTa                       <  1:1647427/101‑1 (MQ=255)

                                                                          ccATCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTa                       <  1:1647530/101‑1 (MQ=255)

                                                                            aTCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTaaa                     <  2:1647574/101‑1 (MQ=255)

                                                                             tCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAg                    <  2:1647692/101‑1 (MQ=255)

                                                                             tCAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAg                    <  2:1647375/101‑1 (MQ=255)

                                                                              cAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGc                   <  1:1647646/101‑1 (MQ=255)

                                                                              cAGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGc                   <  1:1647554/101‑1 (MQ=255)

                                                                               aGGTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCg                  >  1:1647869/1‑101 (MQ=255)

                                                                                ggTTGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGa                 >  1:1647871/1‑101 (MQ=255)

                                                                                  ttGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGAcc               >  2:1647872/1‑101 (MQ=255)

                                                                                   tGTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCt              >  1:1647876/1‑101 (MQ=255)

                                                                                    gTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCg                  >  2:1647877/1‑96 (MQ=255)

                                                                                    gTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTa             <  1:1647728/101‑1 (MQ=255)

                                                                                    gTCGTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTa             <  2:1647630/101‑1 (MQ=255)

                                                                                       gTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAg          <  1:1647295/101‑1 (MQ=255)

                                                                                       gTTAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAg          <  2:1647459/101‑1 (MQ=255)

                                                                                        ttAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAgc         <  1:1647235/101‑1 (MQ=255)

                                                                                        ttAGTAGGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAgc         <  2:1647254/101‑1 (MQ=255)

                                                                                         tagtagGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAgcg        <  1:1647624/101‑1 (MQ=255)

                                                                                           gtagGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAgcggc      <  1:1647673/101‑1 (MQ=255)

                                                                                             agGGGTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGt    >  1:1647881/1‑101 (MQ=255)

                                                                                              ggggTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTc   <  2:1647592/101‑1 (MQ=255)

                                                                                               gggTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCg  <  1:1647468/101‑1 (MQ=255)

                                                                                               gggTACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCg  <  2:1647578/101‑1 (MQ=255)

GCCAAAGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCG  >  NC_000913/1261597‑1261792

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