BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,597,313	(C)_7→8	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	93.7%	426.3 / 21.8	127	pseudogene (682/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand): ref base . (5/3); new base C (73/46); total (78/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

AGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597209‑1597400

agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa‑‑c‑c                                                                                              <  2:159048/100‑1 (MQ=255)

 gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑c‑cc                                                                                             >  1:2595161/1‑100 (MQ=255)

 gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑c‑cc                                                                                             >  2:1588001/1‑100 (MQ=255)

   cTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑c‑cccc                                                                                           >  1:67013/1‑100 (MQ=255)

   cTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑c‑cccc                                                                                           >  2:2734310/1‑100 (MQ=255)

    tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑c‑ccccc                                                                                          >  1:2705290/1‑100 (MQ=255)

    tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑c‑ccccc                                                                                          <  2:1814460/100‑1 (MQ=255)

     gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑c‑cccccc                                                                                         >  2:310987/1‑100 (MQ=255)

          ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑c‑ccccccctcgc                                                                                    >  2:2229758/1‑95 (MQ=255)

           gatcttgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTcgcg                                                                                    <  1:2740282/95‑1 (MQ=255)

           cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTcgcg                                                                                    <  2:1809972/100‑1 (MQ=255)

           cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTcgcg                                                                                    <  1:564264/100‑1 (MQ=255)

           cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTcgcg                                                                                    >  1:802657/1‑100 (MQ=255)

            gTGATGACATAGTCTTGATCCTGCACATTCAGGGTTCTGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑c‑ccccccctcgcaa                                                                                  <  2:2497785/100‑8 (MQ=255)

            gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGa                                                                                   >  1:290206/1‑100 (MQ=255)

              gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAAt                                                                                 >  2:2308/1‑100 (MQ=255)

                tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAg                                                                               <  1:337296/100‑1 (MQ=255)

                 gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGt                                                                              >  2:3263456/1‑100 (MQ=255)

                  aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGACACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATa                                                                                >  2:1157431/1‑97 (MQ=255)

                   cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCg                                                                            >  2:579768/1‑100 (MQ=255)

                   cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCg                                                                            >  1:545770/1‑100 (MQ=255)

                   cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCg                                                                            >  1:836840/1‑100 (MQ=255)

                   cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑C‑CCCCCCTCGCGAATAGTCGc                                                                           >  2:879161/1‑100 (MQ=255)

                    agagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑C‑CCCCCCTCGCGAATAGTCGcc                                                                          <  1:1401018/100‑1 (MQ=255)

                     gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGcc                                                                          >  1:46707/1‑100 (MQ=255)

                          ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTcccc                                                                     <  2:2628659/100‑1 (MQ=255)

                           tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCa                                                                    >  1:3263980/1‑100 (MQ=255)

                             aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGt                                                                  >  2:129435/1‑100 (MQ=255)

                               ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTc                                                                >  2:772130/1‑100 (MQ=255)

                               ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGCTc                                                                >  1:1224035/1‑100 (MQ=255)

                                cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                               >  2:1218635/1‑100 (MQ=255)

                                   cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            <  1:2308/100‑1 (MQ=255)

                                   cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            >  1:569546/1‑100 (MQ=255)

                                   cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                            <  1:3262430/100‑1 (MQ=255)

                                      aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         >  1:2806122/1‑100 (MQ=255)

                                      aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         <  2:2433804/100‑1 (MQ=255)

                                      aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CACCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         >  1:2701866/1‑100 (MQ=255)

                                       ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                        >  1:1600912/1‑100 (MQ=255)

                                        tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                                                       <  1:1597660/100‑1 (MQ=255)

                                          aGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCC‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                      <  1:476429/100‑1 (MQ=255)

                                           gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                                                    >  2:3276876/1‑100 (MQ=255)

                                           gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                                                    >  2:3134091/1‑100 (MQ=255)

                                             gTTCCGATACCGGTGGCATTGCCACCGATAAGCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTACCTTGATTACCAAtt                                                  >  2:690651/1‑100 (MQ=255)

                                              ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                 >  1:890223/1‑100 (MQ=255)

                                              ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                <  1:1539105/100‑1 (MQ=255)

                                              tacgattatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGACTCCCCAGTTCCTTGATTACCAAttt                                                 >  1:1670663/10‑100 (MQ=37)

                                              tacgattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                 >  2:2871910/10‑100 (MQ=38)

                                               acgattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                >  1:2216464/9‑100 (MQ=255)

                                                ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                                               >  2:864024/1‑100 (MQ=255)

                                                 cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                              <  1:1361590/100‑1 (MQ=255)

                                                  nttatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGata                                             >  2:963489/6‑98 (MQ=255)

                                                  gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTTGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             <  1:3182083/100‑1 (MQ=255)

                                                  gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             >  2:1442597/1‑100 (MQ=255)

                                                  gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             <  2:37387/100‑1 (MQ=255)

                                                  attatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             >  2:2417706/6‑100 (MQ=255)

                                                   ttatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                            >  2:2274737/5‑100 (MQ=255)

                                                   aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                            <  1:3175051/100‑1 (MQ=255)

                                                    tatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCCGTTCCTTGATTACCAATTTGAAAtt                                           >  1:1236669/4‑100 (MQ=255)

                                                    tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                           >  1:690741/4‑100 (MQ=255)

                                                    tACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                           >  2:3226347/1‑100 (MQ=255)

                                                     atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             >  1:521591/3‑97 (MQ=255)

                                                     atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                          >  2:3200039/3‑100 (MQ=255)

                                                     atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTCATTACCAATTTGAAAt                                            >  2:671939/3‑98 (MQ=255)

                                                      ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                         >  2:2910743/1‑100 (MQ=255)

                                                       cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa                                        <  2:2639379/100‑1 (MQ=255)

                                                       cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCACAGTCCCTTGATTACCAATTGGAAATTCa                                         >  1:1196991/1‑99 (MQ=255)

                                                        ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                       >  2:1536217/1‑100 (MQ=255)

                                                         gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTCGAAATTCAAtt                                      >  2:1133821/1‑100 (MQ=255)

                                                          tGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg                                     <  1:2424304/100‑1 (MQ=255)

                                                           ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGa                                    <  2:46706/100‑1 (MQ=255)

                                                            gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                   <  1:129436/100‑1 (MQ=255)

                                                             cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCCATTTGAAATTCAATTGAgg                                  >  1:299832/1‑100 (MQ=255)

                                                             cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTAATTACCAATTTGAAATTCAATTGa                                    >  2:1605955/1‑98 (MQ=255)

                                                               ttGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTTATTACCAATTTGAAATTCAATTGAGGaa                                >  1:1141632/1‑100 (MQ=255)

                                                                  ccACCGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCCCGAATAGTCCCCTCCCCAGTTCCTTGCTTACCAATTTGAAATTCAATTGAGAAATc                              >  1:893999/1‑99 (MQ=255)

                                                                     ccGATAATCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatta                          >  2:176739/1‑100 (MQ=255)

                                                                         tAATCGTATTTTCAGCCGTAACCAA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                     >  1:864810/1‑100 (MQ=255)

                                                                           aTCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    >  1:853087/1‑100 (MQ=255)

                                                                           aTCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    <  2:1236665/100‑1 (MQ=255)

                                                                           aTCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    >  2:2578278/1‑100 (MQ=255)

                                                                           aTCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    <  1:905847/100‑1 (MQ=255)

                                                                            tCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   >  1:1212987/1‑100 (MQ=255)

                                                                            tCGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   <  1:1071942/100‑1 (MQ=255)

                                                                             cGTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                  <  1:1975031/100‑1 (MQ=255)

                                                                              gTATTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                 >  2:538522/1‑100 (MQ=255)

                                                                                aTTTTCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg               <  2:266658/100‑1 (MQ=255)

                                                                                 ttttCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc              <  1:176740/100‑1 (MQ=255)

                                                                                  tttCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc             >  1:442979/1‑100 (MQ=255)

                                                                                   ttCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCa            >  1:993159/1‑100 (MQ=255)

                                                                                    tCAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCAc           <  2:971609/100‑1 (MQ=255)

                                                                                    tCAGCCGTAACCAAC‑C‑CCCCCCACGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCAc           >  2:1948913/1‑100 (MQ=255)

                                                                                     cAGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACt          >  1:2045283/1‑100 (MQ=255)

                                                                                      aGCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc         <  1:2229769/100‑1 (MQ=255)

                                                                                       gCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg        >  1:67746/1‑100 (MQ=255)

                                                                                       gCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg        >  1:2552922/1‑100 (MQ=255)

                                                                                       gCCGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg        >  1:2159378/1‑100 (MQ=255)

                                                                                        ccGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGc       >  1:1120055/1‑100 (MQ=255)

                                                                                        ccGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTCTTTTTTATTAGCCACTCGc       >  1:2583294/1‑100 (MQ=255)

                                                                                         cGTAACCAAC‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGcc      <  2:2709004/100‑1 (MQ=255)

                                                                                             accaac‑c‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc  >  1:299744/1‑100 (MQ=255)

AGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑C‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597209‑1597400

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