BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,048,676	A→G	S104P (TCC→CCC)	ymcD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,048,676	0	A	G	99.5%	603.6 / ‑4.7	190	S104P (TCC→CCC)	ymcD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (103/86); total (104/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 196 total aligned reads displayed.

GCAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGACGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTGCC  >  NC_000913/1048587‑1048764

gCAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgt                                                                                       <  1:4377308/93‑1 (MQ=255)

 cAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgtg                                                                                      >  1:5672218/1‑93 (MQ=255)

  aaaTCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgtgg                                                                                     >  1:497146/1‑93 (MQ=255)

  aaaTCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgtgg                                                                                     <  1:3225681/93‑1 (MQ=255)

   aaTCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgtggt                                                                                    >  1:1115473/1‑93 (MQ=255)

   aaTCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgtggt                                                                                    <  1:6021149/93‑1 (MQ=255)

    aTCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgtggtg                                                                                   <  2:4907671/93‑1 (MQ=255)

    aTCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgtggtg                                                                                   >  1:3430343/1‑93 (MQ=255)

     tctcTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGtggtgg                                                                                  <  2:513625/93‑1 (MQ=255)

       tctAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGtt                                                                                >  2:5256244/1‑93 (MQ=255)

        ctAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTg                                                                               >  1:5345985/1‑93 (MQ=255)

         tAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTACTTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGGTGa                                                                              >  1:2194841/1‑93 (MQ=255)

          aaaGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAg                                                                             >  1:69184/1‑93 (MQ=255)

          aaaGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAg                                                                             >  2:3512155/1‑93 (MQ=255)

            aGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGgt                                                                           >  2:4927071/1‑93 (MQ=255)

            aGCCCGGTACGAATATCACAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGgt                                                                           >  2:1165838/1‑93 (MQ=255)

             gCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGgtg                                                                          <  2:4341998/93‑1 (MQ=255)

             acccGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGgtg                                                                          <  1:3130139/92‑1 (MQ=255)

               ccGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGtt                                                                        >  2:6022130/1‑93 (MQ=255)

                cGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTa                                                                       >  1:3444226/1‑93 (MQ=255)

                 ggTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGGTGAGGTGTTAt                                                                      >  1:2787247/1‑93 (MQ=255)

                  gTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATg                                                                     <  2:6418/93‑1 (MQ=255)

                     cGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATgaga                                                                  <  2:1340680/93‑1 (MQ=255)

                     cGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATgaga                                                                  <  1:988545/93‑1 (MQ=255)

                      gAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGaa                                                                 <  2:728180/93‑1 (MQ=255)

                      gAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGaa                                                                 <  1:4044372/93‑1 (MQ=255)

                       aaNATCNCAACCGGGTATCATNNAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAAc                                                                <  1:2323327/93‑1 (MQ=255)

                         tatCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGGTGAGGTGTTATGAGAACCg                                                              >  2:2247558/1‑93 (MQ=255)

                          atCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGc                                                             <  1:4598615/93‑1 (MQ=255)

                          atCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGACGGTTGAGGTGTTATGAGAACCGc                                                             <  2:880872/93‑1 (MQ=255)

                             gCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCAt                                                          <  2:3367697/93‑1 (MQ=255)

                              cAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATc                                                         <  1:428671/93‑1 (MQ=255)

                              cAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATc                                                         <  2:6119013/93‑1 (MQ=255)

                              cAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATc                                                         <  1:3967852/93‑1 (MQ=255)

                              cAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATc                                                         <  2:3631225/93‑1 (MQ=255)

                                aCCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGc                                                       <  2:2163952/93‑1 (MQ=255)

                                 ccGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGcc                                                      >  2:1079148/1‑93 (MQ=255)

                                  cGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCa                                                     >  2:3827286/1‑93 (MQ=255)

                                  cGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCa                                                     <  1:5644633/93‑1 (MQ=255)

                                  cGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCa                                                     <  2:2703824/93‑1 (MQ=255)

                                   gggTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCAc                                                    >  1:5076646/1‑93 (MQ=255)

                                   gggTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCAc                                                    >  1:2733674/1‑93 (MQ=255)

                                   gggTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCAc                                                    >  2:158734/1‑93 (MQ=255)

                                    ggTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACt                                                   <  2:4202161/93‑1 (MQ=255)

                                    ggTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACt                                                   >  1:1284608/1‑93 (MQ=255)

                                    ggTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACt                                                   <  1:2976846/93‑1 (MQ=255)

                                     gTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTg                                                  <  2:4512014/93‑1 (MQ=255)

                                     gTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTg                                                  >  2:4578616/1‑93 (MQ=255)

                                     gTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTg                                                  <  1:351888/93‑1 (MQ=255)

                                     gTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGGTGAGGTGTTATGAGGACCGCCATCGCCACTg                                                  >  2:2613347/1‑93 (MQ=255)

                                      tATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGt                                                 <  2:2495080/93‑1 (MQ=255)

                                        tCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCa                                               <  2:279914/93‑1 (MQ=255)

                                        tCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCa                                               >  1:4436649/1‑93 (MQ=255)

                                          aTTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCAtt                                             >  2:598749/1‑93 (MQ=255)

                                          aTTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCAtt                                             >  2:3082774/1‑93 (MQ=255)

                                           tttAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCAttt                                            >  1:4919062/1‑93 (MQ=255)

                                           tttAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCAttt                                            <  1:6053413/93‑1 (MQ=255)

                                           tttAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCAttt                                            >  2:2381453/1‑93 (MQ=255)

                                              aaCTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTtgct                                         >  2:2016623/1‑93 (MQ=255)

                                                 tttAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCTCTGTACTTTGCTGca                                      >  2:5652149/1‑93 (MQ=255)

                                                     aTTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGcaacaa                                  <  1:2836208/93‑1 (MQ=255)

                                                      tttCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCaacaac                                 <  2:2266482/93‑1 (MQ=255)

                                                        tCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCCACAACTg                               >  2:5588894/1‑93 (MQ=255)

                                                        tCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTg                               >  2:1919345/1‑93 (MQ=255)

                                                          gACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGcc                             >  2:5960093/1‑93 (MQ=255)

                                                          gACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGcc                             <  1:3325150/93‑1 (MQ=255)

                                                            cgcgTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCaa                           <  2:5772895/93‑1 (MQ=255)

                                                             gcgTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAAc                          <  1:4734251/93‑1 (MQ=255)

                                                               gTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGt                        <  1:361938/93‑1 (MQ=255)

                                                               gTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGt                        <  1:1621734/93‑1 (MQ=255)

                                                               gTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGGGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGt                        >  1:3503822/1‑93 (MQ=255)

                                                                  aaaGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACTGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGc                     >  1:5021760/1‑93 (MQ=255)

                                                                  aaaGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGc                     <  1:1669979/93‑1 (MQ=255)

                                                                  aaaGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGc                     >  1:945539/1‑93 (MQ=255)

                                                                  aaaGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGc                     >  1:6104591/1‑93 (MQ=255)

                                                                   aaGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGcc                    >  1:794005/1‑93 (MQ=255)

                                                                   aaGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGcc                    >  1:2104474/1‑93 (MQ=255)

                                                                   aaGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGcc                    <  2:44419/93‑1 (MQ=255)

                                                                    aGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGGCATTTGCTGCAACAACTGCCAACGTCGCca                   >  1:1468238/1‑93 (MQ=255)

                                                                     gCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCcac                  >  2:5028285/1‑93 (MQ=255)

                                                                     gCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCcac                  >  1:3286777/1‑93 (MQ=255)

                                                                     gCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCcac                  <  2:2589383/93‑1 (MQ=255)

                                                                      cAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCcaca                 >  2:1467070/1‑93 (MQ=255)

                                                                      cAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCcaca                 <  1:158/93‑1 (MQ=255)

                                                                       aGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCacac                <  1:1761673/93‑1 (MQ=255)

                                                                        gTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACAcc               >  1:5445944/1‑93 (MQ=255)

                                                                        gTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACAcc               <  1:1888255/93‑1 (MQ=255)

                                                                        gTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACAcc               >  1:5806258/1‑93 (MQ=255)

                                                                        gTATTAGCGTGTAGTGGGCGTGGTGGTTGACGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACAcc               >  2:2947165/1‑93 (MQ=255)

                                                                         tATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCa              >  2:1855015/1‑93 (MQ=255)

                                                                          aTTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCaa             >  2:3760794/1‑93 (MQ=255)

                                                                           ttAGGGTGTAGTGGGCGTGGTTGTTGAGGTGTTATGAGAACCGCCATCGCCCCTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAAc            <  2:2656404/93‑1 (MQ=255)

                                                                             aGCGTGTGGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCTACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAg          <  1:103480/93‑1 (MQ=255)

                                                                              gCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGc         <  2:4706667/93‑1 (MQ=255)

                                                                               ggtgtAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGcc        <  1:4204411/92‑1 (MQ=255)

                                                                                gtgtAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCg       >  1:538513/1‑93 (MQ=255)

                                                                                    aGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTGc   <  1:3692923/93‑1 (MQ=255)

                                                                                     gTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGTCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTGcc  >  1:6143964/1‑93 (MQ=255)

                                                                                     gTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTGcc  >  2:4064365/1‑93 (MQ=255)

                                                                                     gTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTGcc  <  1:5595583/93‑1 (MQ=255)

GCAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGACGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTGCC  >  NC_000913/1048587‑1048764

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