BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,215,052	A→G	E122G (GAG→GGG)	aceA →	isocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,215,052	0	A	G	99.7%	981.2 / NA	307	E122G (GAG→GGG)	aceA	isocitrate lyase
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (137/169); minor base T (1/0); total (138/169)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

GATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCC  >  NC_000913/4214962‑4215141

gATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGc                                                                                         <  2:1898137/93‑1 (MQ=255)

 aTGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCg                                                                                        <  2:8843310/93‑1 (MQ=255)

 aTGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCg                                                                                        <  1:3326023/93‑1 (MQ=255)

     cAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATc                                                                                    <  1:1015089/93‑1 (MQ=255)

     cAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATc                                                                                    <  1:4892445/93‑1 (MQ=255)

     cAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATc                                                                                    <  1:2151517/93‑1 (MQ=255)

      aGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATca                                                                                   >  1:7922126/1‑93 (MQ=255)

       ggTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATcaa                                                                                  <  2:7393633/93‑1 (MQ=255)

       ggTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATcaa                                                                                  >  1:8630061/1‑93 (MQ=255)

       ggTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATcaa                                                                                  >  2:3264899/1‑93 (MQ=255)

        gTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGGCGGATcaa                                                                                  >  2:583930/1‑93 (MQ=255)

         tAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATcaaca                                                                                <  1:8395592/93‑1 (MQ=255)

          aGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATcaacaa                                                                               <  2:130450/93‑1 (MQ=255)

           gcggcgGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCaacaac                                                                              <  2:6611217/93‑1 (MQ=255)

           gcggcgGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCaacaac                                                                              >  1:1378696/1‑93 (MQ=255)

             ggcggACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAacac                                                                            <  2:8342730/93‑1 (MQ=255)

             ggcggACGCTAACCTGGCGGCCAGCATGTAACCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAacac                                                                            <  2:1958843/93‑1 (MQ=255)

              gcggACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACAcc                                                                           <  2:3411039/93‑1 (MQ=255)

               cggACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCt                                                                          >  1:9582032/1‑93 (MQ=255)

                 gACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTc                                                                        <  2:2099268/93‑1 (MQ=255)

                 gACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTc                                                                        >  1:4755283/1‑93 (MQ=255)

                   cGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCcg                                                                      <  1:9932185/93‑1 (MQ=255)

                   cGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCcg                                                                      <  1:5722202/93‑1 (MQ=255)

                    gCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCcgt                                                                     <  2:2239378/93‑1 (MQ=255)

                    gCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCcgt                                                                     <  1:1678412/93‑1 (MQ=255)

                    gCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCGATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCcgt                                                                     <  2:3693245/93‑1 (MQ=255)

                     cTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCcgtc                                                                    >  1:9744960/1‑93 (MQ=255)

                      tAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCcgtcg                                                                   >  1:5992642/1‑93 (MQ=255)

                       aaCCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCcgtcgt                                                                  >  2:5600684/1‑93 (MQ=255)

                        aCCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTg                                                                 >  1:9290197/1‑93 (MQ=255)

                         gcTGGCGGGCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACCCGGTGCCGGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGc                                                                <  1:2649820/92‑1 (MQ=255)

                         ccTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGTGCGGATCAACAACACCTTCCGTCGTGc                                                                >  1:6082960/1‑93 (MQ=255)

                            ggcggcGAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGa                                                             >  2:9880147/1‑93 (MQ=255)

                              cggcCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAATCTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATc                                                           <  2:3070905/93‑1 (MQ=255)

                              cggcCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATc                                                           <  1:2691285/93‑1 (MQ=255)

                               ggcCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCa                                                          >  2:8221132/1‑93 (MQ=255)

                                gcCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAg                                                         >  2:9458729/1‑93 (MQ=255)

                                gcCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAg                                                         <  2:5537319/93‑1 (MQ=255)

                                gcCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAg                                                         <  1:6154368/93‑1 (MQ=255)

                                cccAGCATTTATCCCGATAAGCCGCTCTATAACGAAAACTCGGAGACAGCTGAGGC‑GGACCGGATCAACAACACTTTCCGCCGCGCACATCAg                                                         >  2:4436345/2‑93 (MQ=255)

                                 ccAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGa                                                        >  2:5420633/1‑93 (MQ=255)

                                 ccAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGa                                                        >  2:4181158/1‑93 (MQ=255)

                                 ccAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGa                                                        >  2:4488644/1‑93 (MQ=255)

                                  cAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGAt                                                       >  1:4209164/1‑93 (MQ=255)

                                  cAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGAt                                                       <  2:2732469/93‑1 (MQ=255)

                                  cAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGAt                                                       >  2:9066539/1‑93 (MQ=255)

                                    gCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCATATcc                                                     <  2:4588156/93‑1 (MQ=255)

                                    gCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATcc                                                     >  1:9452755/1‑93 (MQ=255)

                                      aTGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCaa                                                   <  1:6596062/93‑1 (MQ=255)

                                       tGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAAt                                                  <  2:5023716/93‑1 (MQ=255)

                                       tGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAAt                                                  >  2:394815/1‑93 (MQ=255)

                                       tGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAAt                                                  >  2:2920150/1‑93 (MQ=255)

                                       tGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAAt                                                  >  1:3199620/1‑93 (MQ=255)

                                       tGTATCCGGATCAGGCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTACCGTCGTGCCGATCAGATCCAAt                                                  >  2:6910536/1‑93 (MQ=255)

                                        gTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATg                                                 >  1:4851125/1‑93 (MQ=255)

                                         tATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATgg                                                <  2:6436338/93‑1 (MQ=255)

                                          aTCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATggc                                               >  1:3831969/1‑92 (MQ=255)

                                          aTCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGt                                               >  1:6371998/1‑93 (MQ=255)

                                           tCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTc                                              <  2:6338485/93‑1 (MQ=255)

                                            ccGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTcc                                             <  1:4581249/93‑1 (MQ=255)

                                                aTCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCgg                                         >  1:7020719/1‑93 (MQ=255)

                                                aTCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCgg                                         <  1:6264976/93‑1 (MQ=255)

                                                 tCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCggg                                        >  2:7460174/1‑93 (MQ=255)

                                                  cAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGc                                       >  1:8280001/1‑93 (MQ=255)

                                                   aGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCa                                      <  1:2608020/93‑1 (MQ=255)

                                                    gTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCAt                                     <  2:6177753/93‑1 (MQ=255)

                                                       gCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGa                                  >  1:3470170/1‑93 (MQ=255)

                                                        ctctATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCTATGGTCCGCGGGCATTGAg                                 <  1:4647504/93‑1 (MQ=255)

                                                        ctctATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAg                                 >  2:7578608/1‑93 (MQ=255)

                                                         tctATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGc                                >  2:7170332/1‑93 (MQ=255)

                                                         tctATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGc                                <  2:5684318/93‑1 (MQ=255)

                                                         tctATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGc                                <  2:9174908/93‑1 (MQ=255)

                                                           tATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCg                              <  1:2506845/93‑1 (MQ=255)

                                                            aTCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCgg                             >  1:7057758/1‑93 (MQ=255)

                                                             tCCGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCggg                            <  2:2271407/93‑1 (MQ=255)

                                                              cccgCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACCCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGc                           <  2:8568718/90‑1 (MQ=255)

                                                              ccGGCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGc                           <  1:7626689/93‑1 (MQ=255)

                                                                ggCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGa                         <  1:6439287/93‑1 (MQ=255)

                                                                 gCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACACCACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGAt                        <  1:5607140/93‑1 (MQ=255)

                                                                 gCAAACTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGAt                        <  1:4461081/93‑1 (MQ=255)

                                                                   aaaCTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATcc                      <  2:7765690/93‑1 (MQ=255)

                                                                    acttCGGTCGCAGCTGTGGC‑GGGGCGGACCAACTACACCTTCCGTCGTGCCCATCAGACCCAATGGTCCGCGGGCATTGAGCCGGGCGATCcg                     <  1:9683613/90‑1 (MQ=255)

                                                                    aaCTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCcg                     <  1:9049069/93‑1 (MQ=255)

                                                                    aaCTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCcg                     <  1:2910569/93‑1 (MQ=255)

                                                                    aaCTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCcg                     >  2:4868052/1‑93 (MQ=255)

                                                                    aaCTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCcg                     >  1:6814982/1‑93 (MQ=255)

                                                                      cTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCcgcg                   >  1:8922555/1‑93 (MQ=255)

                                                                      cTCGGTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCCGGCATTGAGCCGGGCGATCcgcg                   >  2:9222034/1‑93 (MQ=255)

                                                                         ggTGCCAGCTGTGGT‑GGGGCGGATCAACAACCCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTa                <  1:5496303/93‑1 (MQ=255)

                                                                         ggTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTa                <  1:7424678/93‑1 (MQ=255)

                                                                         ggTGCCAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTa                <  1:2260940/93‑1 (MQ=255)

                                                                              cAGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCg           >  1:883871/1‑93 (MQ=255)

                                                                               aGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGa          <  1:1547468/93‑1 (MQ=255)

                                                                               aGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGa          <  2:1055225/93‑1 (MQ=255)

                                                                               aGCTGTGGT‑GGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGa          >  1:8075461/1‑93 (MQ=255)

                                                                                   gtggt‑gGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTAc      <  1:5222805/93‑1 (MQ=255)

                                                                                     ggt‑ggGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACtt    <  2:1565430/93‑1 (MQ=255)

                                                                                      gt‑ggGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTc   <  2:5173026/93‑1 (MQ=255)

                                                                                       t‑ggGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTcc  >  1:2774188/1‑93 (MQ=255)

                                                                                       t‑ggGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTcc  <  1:1266443/93‑1 (MQ=255)

GATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGT‑GGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCC  >  NC_000913/4214962‑4215141

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