BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,799	T→C	V462A (GTA→GCA)	metL →	aspartokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,799	0	T	C	99.6%	773.6 / ‑5.3	245	V462A (GTA→GCA)	metL	aspartokinase II
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (114/130); total (115/130)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 255 total aligned reads displayed.

TGGTGGCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGT  >  NC_000913/4128709‑4128888

tggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAt                                                                                         <  1:765263/93‑1 (MQ=255)

tggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAt                                                                                         >  1:1929566/1‑93 (MQ=255)

  gtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTg                                                                                       <  1:4075695/93‑1 (MQ=255)

    tgCATTACTGCGCGCCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGtt                                                                                     <  2:1031699/92‑1 (MQ=255)

     gCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTc                                                                                    <  1:3034277/93‑1 (MQ=255)

     gCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTc                                                                                    <  1:1279252/93‑1 (MQ=255)

     gCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTc                                                                                    <  2:67368/93‑1 (MQ=255)

      cAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAGAACGCATCGGCCTGGCATTGTTCg                                                                                   >  1:2707606/1‑93 (MQ=255)

      cAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCg                                                                                   <  1:1664883/93‑1 (MQ=255)

        gTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGt                                                                                 <  2:86958/93‑1 (MQ=255)

        gTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGt                                                                                 >  2:3580132/1‑93 (MQ=255)

        gTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGt                                                                                 <  2:4662619/93‑1 (MQ=255)

           cTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAg                                                                              <  2:188363/93‑1 (MQ=255)

           cTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAg                                                                              <  1:2644242/93‑1 (MQ=255)

           cTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAg                                                                              <  2:2429336/93‑1 (MQ=255)

            tGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAgg                                                                             <  1:3812122/93‑1 (MQ=255)

            tGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAgg                                                                             >  1:2244963/1‑93 (MQ=255)

            tGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTT‑CGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAggg                                                                            <  1:919707/93‑1 (MQ=255)

              cgcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGc                                                                           <  1:4580420/93‑1 (MQ=255)

               gcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCa                                                                          <  1:241470/93‑1 (MQ=255)

                gaCCGCCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCaa                                                                         <  2:581753/92‑1 (MQ=255)

                cACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCaa                                                                         >  2:2462364/1‑93 (MQ=255)

                  ccGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAata                                                                       <  1:93307/93‑1 (MQ=255)

                  ccGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAata                                                                       <  2:1207493/93‑1 (MQ=255)

                  ccGCCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAata                                                                       <  2:3006416/93‑1 (MQ=255)

                   cGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAatat                                                                      <  2:4265725/93‑1 (MQ=255)

                    ggCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATc                                                                     <  1:4429673/93‑1 (MQ=255)

                     gCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCg                                                                    >  1:291754/1‑93 (MQ=255)

                      cccGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCgg                                                                   <  2:2690900/93‑1 (MQ=255)

                         gaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTc                                                                >  2:1589997/1‑93 (MQ=255)

                           ccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTccc                                                              <  1:1144642/93‑1 (MQ=255)

                           ccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTccc                                                              <  2:1068432/93‑1 (MQ=255)

                           ccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTccc                                                              >  2:3876577/1‑93 (MQ=255)

                                 tccTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGCGCAATATCGGTTCCCGTTGGc                                                        <  2:1821371/92‑1 (MQ=255)

                                 gCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGc                                                        >  1:3052237/1‑93 (MQ=255)

                                  ccTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCt                                                       <  2:2067850/93‑1 (MQ=255)

                                  ccTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCt                                                       <  1:359625/93‑1 (MQ=255)

                                  ccTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCt                                                       <  2:2490699/93‑1 (MQ=255)

                                  ccTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCt                                                       <  2:2734564/93‑1 (MQ=255)

                                      aTTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGaa                                                   <  2:3271852/93‑1 (MQ=255)

                                        tCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACt                                                 <  1:3889623/93‑1 (MQ=255)

                                         cAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTg                                                >  1:3744649/1‑93 (MQ=255)

                                         cAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTg                                                >  2:3681500/1‑93 (MQ=255)

                                         cAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTg                                                <  2:1472668/93‑1 (MQ=255)

                                           ggggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGtt                                              >  2:1365566/1‑93 (MQ=255)

                                             ggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCg                                            <  1:2750430/93‑1 (MQ=255)

                                               cTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAAATGTTCgac                                          >  2:4179722/1‑91 (MQ=255)

                                                tGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGccc                                         >  1:1772834/1‑93 (MQ=255)

                                                tGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGccc                                         >  2:2783464/1‑93 (MQ=255)

                                                 gCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCg                                        <  1:2001156/93‑1 (MQ=255)

                                                  cATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGt                                       >  1:559026/1‑93 (MQ=255)

                                                   aTCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGGACTGTTCGCCCGTg                                      >  2:3941366/1‑93 (MQ=255)

                                                    taaGTCCGTCTTCCGCGCAGAAAAGCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTga                                     <  2:3864365/91‑1 (MQ=255)

                                                    tCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTga                                     <  1:205067/93‑1 (MQ=255)

                                                     cAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgag                                    >  1:3274403/1‑93 (MQ=255)

                                                      aGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagc                                   <  1:4384930/93‑1 (MQ=255)

                                                      aGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagc                                   >  1:1409623/1‑93 (MQ=255)

                                                       gTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagca                                  <  2:3342811/93‑1 (MQ=255)

                                                       gTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagca                                  >  1:4195132/1‑93 (MQ=255)

                                                        tCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcag                                 >  1:2995344/1‑93 (MQ=255)

                                                         ccGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcaga                                <  1:706211/93‑1 (MQ=255)

                                                         ccGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcaga                                >  1:2792806/1‑93 (MQ=255)

                                                          cGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagag                               >  1:1041332/1‑93 (MQ=255)

                                                          cGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagag                               <  1:248905/93‑1 (MQ=255)

                                                          cGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagag                               <  2:2313705/93‑1 (MQ=255)

                                                           gTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagagc                              >  2:2591553/1‑93 (MQ=255)

                                                            tCTTCCTCGCAGAAAAACGCATCGGCCTAGCATTGTTCGGGAAGGGCAATATCGGTTTCCGTTGGCTGGAACTGTTCGGGCGTgaacagagca                             >  2:2861536/1‑93 (MQ=255)

                                                            tCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagagca                             <  2:4037188/93‑1 (MQ=255)

                                                            tCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagagca                             <  2:2181437/93‑1 (MQ=255)

                                                            tCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagagca                             <  1:2279787/93‑1 (MQ=255)

                                                             cTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCAc                            <  2:808868/93‑1 (MQ=255)

                                                             cTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTACGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCAc                            <  2:1518624/93‑1 (MQ=255)

                                                              ttCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACg                           >  1:3366629/1‑93 (MQ=255)

                                                                ccGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCt                         >  2:1337308/1‑93 (MQ=255)

                                                                   cgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGGCCGTGAGCAGAGCACGCTTTc                      >  2:37013/1‑93 (MQ=255)

                                                                   cgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTc                      <  2:2238367/93‑1 (MQ=255)

                                                                   cgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTc                      <  2:3571804/93‑1 (MQ=255)

                                                                   cacAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTc                      <  1:4512520/91‑1 (MQ=255)

                                                                       gAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCa                  >  1:1888891/1‑93 (MQ=255)

                                                                         aaaaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACg                >  2:2035999/1‑93 (MQ=255)

                                                                         aaaaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACg                >  1:3200135/1‑93 (MQ=255)

                                                                         aaaaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACg                >  2:3956763/1‑93 (MQ=255)

                                                                         aaaaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACg                <  2:1597930/93‑1 (MQ=255)

                                                                           aaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTa              >  2:2547963/1‑93 (MQ=255)

                                                                           aaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTa              >  2:441182/1‑93 (MQ=255)

                                                                           aaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACCTa              >  1:1447584/1‑93 (MQ=255)

                                                                            aCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTAc             <  2:288781/93‑1 (MQ=255)

                                                                              gCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCg           >  2:520726/1‑93 (MQ=255)

                                                                              gCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCg           <  1:1506996/93‑1 (MQ=255)

                                                                               cATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCgg          >  1:3904990/1‑93 (MQ=255)

                                                                               cATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCgg          >  2:97845/1‑93 (MQ=255)

                                                                                aTCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGc         >  2:1199065/1‑93 (MQ=255)

                                                                                 tCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCt        >  2:346175/1‑93 (MQ=255)

                                                                                  cGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCtt       <  2:1367097/93‑1 (MQ=255)

                                                                                   ggCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCttt      <  1:1583911/93‑1 (MQ=255)

                                                                                    gCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTg     <  2:1980558/93‑1 (MQ=255)

                                                                                    gCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTg     >  2:4078647/1‑93 (MQ=255)

                                                                                      cTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAg   >  2:3073488/1‑93 (MQ=255)

                                                                                       tGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGt  >  1:2153935/1‑93 (MQ=255)

                                                                                       tGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGt  >  1:1229982/1‑93 (MQ=255)

TGGTGGCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGT  >  NC_000913/4128709‑4128888

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