BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,597,313:1	(C)_7→8	100%	pseudogene (675/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,306	1	.	C	55.1%	199.7 / 210.6	107	pseudogene (682/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand): ref base . (29/19); new base C (35/24); total (64/43)
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.26e-01

Only 100 of 118 total aligned reads displayed.

GAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTT  >  NC_000913/1597208‑1597404

 agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa‑cc                                                                                                  <  1:710534/100‑1 (MQ=255)

   tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑cccc                                                                                                <  2:3296613/100‑1 (MQ=255)

      gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑ccccccc                                                                                             >  1:1073506/1‑100 (MQ=255)

        gTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTc                                                                                           <  1:1211203/100‑1 (MQ=255)

         tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑cccccccctc                                                                                          >  2:3203954/1‑97 (MQ=255)

           ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTcgcg                                                                                        <  2:1259933/100‑1 (MQ=255)

            cGTGATGACAGAGTCTTGATCCTGCACATGCAGGGTTCCGATACCGGGGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTcgcg                                                                                        <  2:789627/100‑1 (MQ=255)

             gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGaa                                                                                      >  2:1123813/1‑100 (MQ=255)

              tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGaa                                                                                      <  1:1313285/100‑1 (MQ=255)

               gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATa                                                                                    >  1:761517/1‑100 (MQ=255)

               gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATa                                                                                    >  1:3393936/1‑100 (MQ=255)

                atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAg                                                                                   >  1:3016842/1‑100 (MQ=255)

                  gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTc                                                                                 >  1:2550925/1‑100 (MQ=255)

                  gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTc                                                                                 >  2:1147419/1‑100 (MQ=255)

                   ncAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAG‑CCCCCCCTCGCGAATAGTCg                                                                                >  2:738712/2‑100 (MQ=255)

                   aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCg                                                                                <  2:2514645/100‑1 (MQ=255)

                    cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTc                                                                                 >  1:1239123/1‑99 (MQ=255)

                    cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCg                                                                                >  1:2364730/1‑100 (MQ=255)

                        gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTc                                                                            >  2:2771453/1‑99 (MQ=255)

                         tCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTccc                                                                          <  2:3445583/100‑1 (MQ=255)

                         tCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTccc                                                                          >  1:2593400/1‑100 (MQ=255)

                         tCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTcc                                                                           >  1:2648414/1‑99 (MQ=255)

                          cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGACTAGTCGCCTccc                                                                          >  2:1123584/1‑100 (MQ=255)

                          cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTccc                                                                          >  2:62493/1‑100 (MQ=255)

                          cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTcccc                                                                         >  2:2998799/1‑100 (MQ=255)

                          cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTcccc                                                                         >  2:2989824/1‑100 (MQ=255)

                          cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTcccc                                                                         >  2:1161694/1‑100 (MQ=255)

                           ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTcccc                                                                         <  1:3201731/100‑1 (MQ=255)

                           ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCa                                                                        >  2:2725277/1‑100 (MQ=255)

                            tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTTAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTcccc                                                                         >  1:3340810/1‑99 (MQ=255)

                              aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGt                                                                      >  2:1151825/1‑100 (MQ=255)

                                  tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                                 >  2:3442084/1‑100 (MQ=255)

                                  tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                                 >  2:3442040/1‑100 (MQ=255)

                                   gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                                <  1:3006281/100‑1 (MQ=255)

                                    cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                                                               <  1:2989800/100‑1 (MQ=255)

                                     acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                                                              <  1:231837/100‑1 (MQ=255)

                                       aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                             <  1:988425/100‑1 (MQ=255)

                                       aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                            <  1:1123803/100‑1 (MQ=255)

                                        ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                            >  2:1016902/1‑100 (MQ=255)

                                         tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                                                           <  2:857189/100‑1 (MQ=255)

                                         tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTAATCCGTAACGAA‑CCACCCCTCGCGAATAGTCGTCTCCCCAGGTCCTTGCCTAcc                                                          >  1:153148/1‑100 (MQ=255)

                                          cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                          <  1:2567245/100‑1 (MQ=255)

                                            gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAt                                                       >  2:133646/1‑100 (MQ=255)

                                             ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAt                                                       >  1:616571/1‑100 (MQ=255)

                                             ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                      <  2:883103/100‑1 (MQ=255)

                                             ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                      <  2:2521772/100‑1 (MQ=255)

                                              gTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                      >  2:1218956/1‑100 (MQ=255)

                                                tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                    <  1:2172511/100‑1 (MQ=255)

                                                 ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                                  <  2:2593419/100‑1 (MQ=255)

                                                 ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                                  >  1:565136/1‑100 (MQ=255)

                                                  cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                                 <  2:1768642/100‑1 (MQ=255)

                                                   gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                <  1:989369/100‑1 (MQ=255)

                                                   gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                <  2:3016866/100‑1 (MQ=255)

                                                    aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                <  1:3235797/100‑1 (MQ=255)

                                                    aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                <  1:3005924/100‑1 (MQ=255)

                                                    aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                <  1:2847804/100‑1 (MQ=255)

                                                    aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                <  1:1317462/100‑1 (MQ=255)

                                                       ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                             >  1:1257317/1‑100 (MQ=255)

                                                       ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa                                            >  2:1505651/1‑100 (MQ=255)

                                                       ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa                                            <  1:133646/100‑1 (MQ=255)

                                                        cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                           <  1:275544/100‑1 (MQ=255)

                                                         ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                           <  2:3322129/100‑1 (MQ=255)

                                                         ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                           >  2:421854/1‑100 (MQ=255)

                                                         ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                          >  2:1546028/1‑100 (MQ=255)

                                                           tGGCATTGCCACCGATAATCGTATATTCAGCCGTAAACAACCCCCCCCTCGCAAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGCAATTCAAt                                           <  1:1180521/98‑1 (MQ=255)

                                                            ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGa                                        >  2:72090/1‑100 (MQ=255)

                                                             gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                       <  2:1549948/100‑1 (MQ=255)

                                                              cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                      <  1:1809901/100‑1 (MQ=255)

                                                              cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                      >  1:719020/1‑100 (MQ=255)

                                                                ttGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGACTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa                                    >  1:2207497/1‑100 (MQ=255)

                                                                ttGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAAt                                   >  2:25329/1‑100 (MQ=255)

                                                                    cACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCACCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                                >  1:885227/1‑100 (MQ=255)

                                                                    cACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                               >  2:2351058/1‑100 (MQ=255)

                                                                       cGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattat                             >  2:929114/1‑100 (MQ=255)

                                                                           aaTCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                        >  1:2775136/1‑100 (MQ=255)

                                                                           aaTCGTATTTCCAGCCGTAACCAA‑CCCCCCCTCGCGTATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                        <  1:1240275/100‑1 (MQ=255)

                                                                             tCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                       >  2:448814/1‑100 (MQ=255)

                                                                              cGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTTAGGAATCATTATTTTttat                      >  2:2415378/1‑100 (MQ=255)

                                                                              cGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                      <  1:30146/100‑1 (MQ=255)

                                                                               ntATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                     >  2:3043779/2‑100 (MQ=255)

                                                                                 aTTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg                   <  1:2800950/100‑1 (MQ=255)

                                                                                 aTTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg                   <  2:2127310/100‑1 (MQ=255)

                                                                                  ttttCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc                  <  2:616575/100‑1 (MQ=255)

                                                                                  ttttCAGCCGTAACCAA‑CCCCCCCTCGCGACTAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc                 >  1:1058971/1‑100 (MQ=255)

                                                                                    ttCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCa                >  2:2134001/1‑100 (MQ=255)

                                                                                    ttCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCAc               >  1:3247347/1‑100 (MQ=255)

                                                                                      cAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACt              >  1:3386150/1‑100 (MQ=255)

                                                                                      cAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc             >  2:1608059/1‑100 (MQ=255)

                                                                                       aGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGCTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCa                >  1:3126531/1‑97 (MQ=255)

                                                                                        gCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc             <  2:1443743/99‑1 (MQ=255)

                                                                                        gCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg            >  2:218353/1‑100 (MQ=255)

                                                                                        gCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGc           <  2:480631/100‑1 (MQ=255)

                                                                                         ccGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGc           >  2:2626660/1‑100 (MQ=255)

                                                                                           gTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc        >  1:263481/1‑100 (MQ=255)

                                                                                            tAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc        >  2:2650073/1‑100 (MQ=255)

                                                                                            tAACCAA‑CCACCCCTCGCGAATAGTCGCCTACCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc        <  1:1013495/99‑1 (MQ=255)

                                                                                             aaccaaccCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCatc       <  2:2346644/100‑3 (MQ=255)

                                                                                              accaaccCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc      >  1:581568/1‑100 (MQ=255)

                                                                                              accaa‑ccCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCt     >  1:46267/1‑100 (MQ=255)

                                                                                                 aa‑ccCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATTATTATTTTTTATTAGCCACTCGCCACCCtttt  <  2:3281869/100‑1 (MQ=255)

GAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTT  >  NC_000913/1597208‑1597404

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