BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	404,504	A→G	I268V (ATT→GTT)	dgcC →	diguanylate cyclase DgcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	404,504	0	A	G	100.0%	1141.9 / NA	268	I268V (ATT→GTT)	dgcC	diguanylate cyclase DgcC
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (151/117); total (151/117)

Only 100 of 268 total aligned reads displayed.

GTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTAC  >  NC_000913/404406‑404601

gTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGtt                                                                                                 <  1:605165/101‑1 (MQ=255)

gTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGtt                                                                                                 >  1:605642/1‑101 (MQ=255)

 tCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTtg                                                                                                >  1:605648/1‑101 (MQ=255)

 tCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTtg                                                                                                >  2:605649/1‑101 (MQ=255)

  cggcggCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTtgt                                                                                               >  2:605653/1‑101 (MQ=255)

  cggcggCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTtgt                                                                                               <  1:605370/101‑1 (MQ=255)

  cggcggCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTtgt                                                                                               <  1:605285/101‑1 (MQ=255)

   ggcggcATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTtgtg                                                                                              >  1:605654/1‑101 (MQ=255)

   ggcggcATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTtgtg                                                                                              <  2:605146/101‑1 (MQ=255)

    gcggcATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTgg                                                                                             <  1:605204/101‑1 (MQ=255)

      ggcATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGgcg                                                                                           <  2:605290/101‑1 (MQ=255)

       gcATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTtg                                                                                                <  1:605306/95‑1 (MQ=255)

       gcATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGgcgc                                                                                          >  1:605657/1‑101 (MQ=255)

        cATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCt                                                                                         >  1:605659/1‑101 (MQ=255)

        cATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCt                                                                                         >  1:605662/1‑101 (MQ=255)

         aTAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCtt                                                                                        <  2:605334/101‑1 (MQ=255)

           aaTCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTAc                                                                                      >  1:605666/1‑101 (MQ=255)

              cgcgATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCg                                                                                   <  2:605206/101‑1 (MQ=255)

               gcgATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGa                                                                                  >  1:605669/1‑101 (MQ=255)

                cgATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGAc                                                                                 >  2:605671/1‑101 (MQ=255)

                  aTGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAg                                                                               >  1:605673/1‑101 (MQ=255)

                  aTGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAg                                                                               <  1:605242/101‑1 (MQ=255)

                    gCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGtt                                                                             >  1:605678/1‑101 (MQ=255)

                    gCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGtt                                                                             >  2:605679/1‑101 (MQ=255)

                      aaCGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGAc                                                                                 >  2:605683/1‑95 (MQ=255)

                      aaCGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTAc                                                                           >  1:605682/1‑101 (MQ=255)

                           tACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAAtt                                                                      >  2:605690/1‑101 (MQ=255)

                           tACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAAtt                                                                      >  1:605688/1‑101 (MQ=255)

                           tACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAAtt                                                                      <  2:605259/101‑1 (MQ=255)

                            aCTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTa                                                                     <  1:605347/101‑1 (MQ=255)

                            aCTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTa                                                                     <  1:605323/101‑1 (MQ=255)

                            aCTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTa                                                                     >  2:605693/1‑101 (MQ=255)

                              tGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTAcc                                                                   >  1:605696/1‑101 (MQ=255)

                               gATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTAccc                                                                  >  1:605698/1‑101 (MQ=255)

                                aTTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCt                                                                 <  1:605274/101‑1 (MQ=255)

                                 ttATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTg                                                                <  2:605317/101‑1 (MQ=255)

                                 ttATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTg                                                                <  2:605338/101‑1 (MQ=255)

                                 ttATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTg                                                                >  2:605700/1‑101 (MQ=255)

                                 ttATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTg                                                                <  1:605053/101‑1 (MQ=255)

                                 ttATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTg                                                                >  1:605701/1‑101 (MQ=255)

                                  tatcgatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTgc                                                               <  1:605411/101‑1 (MQ=255)

                                   atcgatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTgcg                                                              >  1:605704/1‑101 (MQ=255)

                                    tcgatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTgcgc                                                             <  1:605266/101‑1 (MQ=255)

                                      gatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCgg                                                           <  2:605199/101‑1 (MQ=255)

                                      gatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCgg                                                           >  1:605708/1‑101 (MQ=255)

                                      gatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCgg                                                           >  2:605709/1‑101 (MQ=255)

                                       atatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGt                                                          <  1:605310/101‑1 (MQ=255)

                                        tatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTa                                                         >  2:605713/1‑101 (MQ=255)

                                         atcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAg                                                        <  2:605522/101‑1 (MQ=255)

                                          tcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGc                                                       >  2:605716/1‑101 (MQ=255)

                                           cgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTgcg                                                              >  2:605720/1‑93 (MQ=255)

                                           cgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCg                                                      >  1:605718/1‑101 (MQ=255)

                                            gaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGa                                                     >  2:605723/1‑101 (MQ=255)

                                             aCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGAt                                                    >  2:605725/1‑101 (MQ=255)

                                               cATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGAtgt                                                  >  1:605729/1‑101 (MQ=255)

                                               cATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGAtgt                                                  <  2:605592/101‑1 (MQ=255)

                                                aTTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGAtgtg                                                 <  2:605481/101‑1 (MQ=255)

                                                 tttCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGa                                                >  1:605732/1‑101 (MQ=255)

                                                   tCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGAtt                                              >  2:605733/1‑101 (MQ=255)

                                                    cAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTg                                             >  2:605738/1‑101 (MQ=255)

                                                     aaGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTgg                                            >  2:605741/1‑101 (MQ=255)

                                                     aaGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTgg                                            >  1:605740/1‑101 (MQ=255)

                                                         gCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCgg                                        <  1:605420/101‑1 (MQ=255)

                                                         gCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCgg                                        <  2:605094/101‑1 (MQ=255)

                                                         gCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCgg                                        <  1:605145/101‑1 (MQ=255)

                                                          cATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGt                                       >  2:605744/1‑101 (MQ=255)

                                                            tCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGttt                                     <  1:605445/101‑1 (MQ=255)

                                                            tCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGttt                                     <  2:605474/101‑1 (MQ=255)

                                                             cAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTg                                    <  1:605172/101‑1 (MQ=255)

                                                              aaCGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTgg                                   >  1:605749/1‑101 (MQ=255)

                                                               aCGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCg                                         >  1:605751/1‑94 (MQ=255)

                                                                cGATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTggcg                                 >  1:605754/1‑101 (MQ=255)

                                                                 gATACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTggcgg                                >  1:605757/1‑101 (MQ=255)

                                                                  aTACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTggcggc                               <  2:605277/101‑1 (MQ=255)

                                                                  aTACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTggcggc                               <  2:605080/101‑1 (MQ=255)

                                                                   tACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTg                                    <  1:605188/95‑1 (MQ=255)

                                                                   tACCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGgcggcg                              >  2:605758/1‑101 (MQ=255)

                                                                    aCCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTggcg                                 >  2:605760/1‑97 (MQ=255)

                                                                    aCCTGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGa                             <  2:605367/101‑1 (MQ=255)

                                                                       tGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGgcggcg                              >  2:605766/1‑97 (MQ=255)

                                                                       tGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGa                          >  1:605765/1‑101 (MQ=255)

                                                                       tGGGGCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGa                          <  2:605452/101‑1 (MQ=255)

                                                                        ggggCCATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAg                         >  1:605768/1‑101 (MQ=255)

                                                                            ccATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTg                     <  2:605423/101‑1 (MQ=255)

                                                                             cATGATGTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGc                    >  2:605772/1‑101 (MQ=255)

                                                                              atgatgTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCa                   >  1:605774/1‑101 (MQ=255)

                                                                               tgatgTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGttt                      <  1:604941/97‑1 (MQ=255)

                                                                               tgatgTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAg                  <  1:605380/101‑1 (MQ=255)

                                                                                 atgTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTa                <  1:605484/101‑1 (MQ=255)

                                                                                 atgTGGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTa                >  2:605778/1‑101 (MQ=255)

                                                                                   gtgGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAAt              >  1:605779/1‑101 (MQ=255)

                                                                                   gtgGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAAt              <  1:605340/101‑1 (MQ=255)

                                                                                    tgGGCGATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATc             >  1:605781/1‑101 (MQ=255)

                                                                                         gATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTc        >  2:605786/1‑101 (MQ=255)

                                                                                         gATGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTc        >  1:605787/1‑101 (MQ=255)

                                                                                          aTGAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTcc       <  1:605387/101‑1 (MQ=255)

                                                                                            gAAGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCgg     <  1:605520/101‑1 (MQ=255)

                                                                                             aaGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGt    >  2:605790/1‑101 (MQ=255)

                                                                                              aGCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTa   <  1:605216/101‑1 (MQ=255)

                                                                                               gCGGTTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTAc  <  2:605417/101‑1 (MQ=255)

GTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTAC  >  NC_000913/404406‑404601

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